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Among the methods available to determine the complete 
state of stress at any point in a three-dimensional body, photo 
elasticity ranks very high. In photoelasticity, there are 
essentially two methods, (a) using transmitted light, and 
(b) using scattered light. The transmitted light technique 
presently needs a frozen stress model which has to be cut into 
slices and sub-slices to obtain the required data. Hence, the 
method is .destructive in nature. The scattered light technique 
however is nondestructive and has certain other advantages in 
the sense that no frozen st:^ess model is necessary. The models 
can be subjected to live loads at room temperatures and 
analysed. 


In scattered light techniqTj.e there are essentially 
two methods, one developed by Frocht and Srinath, and the other 
developed by Robert and Guillemet. The method developed by 
Frocht and Srinath is limited to cases where the effects of 
stress -induced optical activity is not high. This is equi- 
valent to saying that the rate of rotation of the secondary 
principal-stresses along the light iHth is small compared to 
the rate of retardation. The method developed by Robert and 
Guillemet does not involve any assumption regarding optical 
activity, but requires certain characteristic parameters which 
are indirectly related to the secondary principal-stresses and 
their directions. 

The present investigation is concerned with the 
following topics; 

1) Development of relevant optical equations using Stokes 
parameters for the analysis of problems involving optical 
activity, 

2) Experimental techniques to determine the characteristic 
parameters, that is the primary and secondary characteristic 
directions and the characteristic retardation; 

3) Reduction of a three-dimensional photoelastic model to an 
optically equivalent system; 

4-) New methods for scattered light analysis without involving 

* 

assumptions regarding stress-induced optical activity; 



5) Development of new graphical methods to determine the 
integrated effect in problems involving optical activity; 

6) Development of a large-field compensator using a few 
plates of arbitrary but equal retardation -an extension of 
Pancharatnam ' s andTuzi’s worh; 

7) Determination of integral fringe orders using different 
wavelengths of light; 

8) Experimental verification of the above methods using a 
model under combined torsion and tension; 

A short write-up on the above topics is giveb belowi 

The usual photoelastic analysis consists of trigono- 
metrical resolutions of the light vector at various stages, 
as it passes through the model* However, the study of light 
ellipse and its properties as it passes through the various 
optical elements provides a better and elegant method of 
analysis. This concept has been used throughout this inves- 
tigation. To start with, various previous systems of equation 
of three-dimensional photbelasticity are studied. Then these 
optical equations are kinematically derived in terms of Stokes 
parameters. Since the transformations of these parameters 
when the light ellipse passes through a retarder are linear 
in these parameters, the resulting optical equations derived 
in terms of Stokes parameters turn out to be a set of ordinary 
linear differential equatiohs which are very convenient for 



a.na,lysis and application. Using these equations, the analytical 
solution for the case where the ratio of the rate of rotation 
to the rate of retardation is a constant, is obtained (The rate 
of rotation itself can be any arbitrary function). This case is 
more general compared to the case investigated by Drucker and 
Hindlin where each of these rates is a constant. 

From these equations, the equivalence of any optical 
path in a three-dimensional photoelastic model to a system 
consisting of a retarder and a pure rotator follows directly. 

Two new direct methods of finding these characteristic para- 
meters experimentally in a transmission polariscope are presented. 
Extension of these methods to the scattered light technique is 
described. These methods of scattered light analysis do not 
involve any assumptions regarding optical activity. Also 
these do not require the rotation of the analyser at a parti- 
cular speed as is required in the method of Robert and 
Guillemet. A direct and whole field method which is derived 
from the optical equations in terms of Stokes parameters is 
presented. This method again is general and it does not in- 
volve any assumptions. 

Further some investigations are carried out on a few 
related topics. The graphical methods are very convenient to 
find the integrated optical effect when the secondary principal- 
stress difference and their directions are given. The existing 
graphical methods are essentially based on stereographic or 
parallel projection of the Poincare sphere or the quaternions. 


A brief review of these methods is presented. A more direct 
and simpler approach to develop the graphical methods which 
is based on the properties of the light ellipse, is then 
presented. Two simple graphical methods which have certain 
advantages over the previous methods are described. 

The compensator is an indispensable optical element 
in the photoelastic work. Construction of a large field com- 
pensator using a few plates of arbitary but same retardation 
is described. This is essentially an extension of the works 
of Pancharatnam and Tuzi. An extension of Pant’s method to 
find the integral fringe order using the fractional fringe 
order for different wavelengths of light, is presented. 

Some experimental results obtained from a rectangular 
photoelastic model under combined tension and torslcn in a 
transmission polariscope are presented as supporting evidence 
for the existence of an optically equivalent system, the two 
methods of determining these parameters experimentally, limita 
tions of the previous scattered light methods and the determi- 
nation of the integral fringe orders. 



CHAPTER 1 


im’RODUCTION 

Phot oelasticity and Electromagnetic Theory of Light 

Photoelasticity deals with the study of propaga- 
tion of light through elastic bodies under load in order 
to determine the stress distribution inside these bodies. 
Phot oplast icity, photoviscoelasticity, dynamic photoelas- 
ticity etc. are immediate extensions of the basic princi- 
ples of photoelasticity. The basic phenomenon of double 
refraction on which photoelasticity stands was discovered 
in 1816 by David Brewster, Neumann formulated the stress- 
-optic law in terms of strain in 18+1 . laber after ten 
years, the same law was formiiiated in terms of stress by 
Maxwell. The first engineering applications were tried 
by French engineer Mesnager and the British scientists 
Coker and Filon in the beginning of the present century. 

The electromagnetic theory of light explains all 

the photoelastic phenomena. According to this theory'', 

it is assumed that the transmission of all forms of radiant 
/ 

energy (including light) through space takes place in the 
form of electromagnetic disturbances. The radiant energy 
of light is equally divided between electric energy and 
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magnetic energy. At any instant, these electric and. mag- 
netic fields at a given point in a ray of light may he repre- 
sented hy two equal and mutually perpendicular vectors called 
electric and magnetic vectors in a plane perpendicular to 
the ray of light . Because of the interd^endence of these 
two vectors, we can for all practical purposes in the photo- 
elastic analysis, consider only cne vector say, the electric 
vector. The electromagnetic spectrum of radj-ant energy 
covers in addition to the visible light rays, the heat rays, 
X-rays, radio ■v>a,ves and the shorter gamma rays. In ordinary 
or unpolarized light, the light vector (electric vector) has 
no directional preference and its changes are chaotic. When 
the electric' vector changes its azimuth (orientation with 
respect to some fixed axes) and magnitude according to a 
certain law, we have polarized light. When the tip of the 
electric vector traces out an ellipse or a circle, it is 
elliptically polarized light or circularly polarized light 
respectively. When the azimuth of light vector does not 
change, that is, the electric vector stays in a particular 
plane, it is plane polarized light. This phenomenon of 
polarization is observed throughout the electromagnetic 
spectrum. 

Stress-Optic Law in Two and Three -Dimensions 

In a transparent , isotropic plate, in which the 
stresses are two-dimensional, and within the elastic limit, 
the phase difference or relative retardation p in wavelengths 



3 


between the rectangular wave components of light travelling 
through it and produced by temporary birefringence is given 
by 

p = h (p - q) (1 ,1 ) 

in which is a constant known as the stress- 

optic coefficient for the particular material for wavelength 
X , h is the thickness of the plate and p and q the prin- 
cipal-stresses. This temporary retardation or phase diffe- 
rence (this vanishes when the loads are removed) is intro- 
duced along the principal-stress axes. 

Stress-Qptic Law in Three -Dimens ions 

In a three-dimensional photoelastic model, a retar- 
dation (in wavelengths or fringes) equal to 

dn = dh (p' - q' ) (1 .2) ' 

will be added along the secondary principal-stress 
axes to the phase of the components of light vector along 
these directions. C-^ is the stress-optic coefficient as 
in two-dimensional case, dh is the actual length of light 
path (dh is so small that the rotation of the secondary 
principal-stress axes over this distance is negligible), p' , 
q' are the secondary principal-stresses for the direction 
of the given ray. The secondary principal-stresses for a 
given direction of light propagation are defined as the 
principal-stresses resulting from the stress components 
which lie in a plane normal to the direction of light ray. 

In both the above cases, the stress-optic coefficient may 




also depend upon the temperature of the model material. 

Thus, the photoelastic analysis essentially involves 
the measurement of retardation introduced in the model and 
finding the directions along which this retardation is intro- 
duced. A photoelastic model under an isotropic state of 
stress will not introduce any retardation and a general 
state of stress can he broken into two parts, the isotropic 
part and the deviatoric part. So, the photoelastic analysis 

will yield five stress components at a point out of the total 
2 

six components . The other unknown can be found by using 

either the shear difference method^, oblique incidence'^ , 

2 S 6 

electric analogy or by an Interferometric method 
The photoelastic technique in the two dimensional case is 
almost perfected. In this case, the model making and the 
simulation of the loading of the prototype take much of the 
time while the experiment is direct and simple. The bire- 
fringent coating method and photoelastic study with a ref- 
lection polarise ope eliminate model making without loosing 
the advantages of the transmission method"^ . 

Three-Dimensional Photoelastic Techniques and Rotational 
Effect s t 

In the general three-dimensional cases of stress 
distribution, the secondary principal-stresses vary conti- 
nuously along the wavenormal and this makes the analysis 
difficult. The three-dimensional photoelastic techniques 
in vogue are essentially two, (i) the stress freezing 
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technique where the stresses are frozen into the model 
and the model is carefully cut into very thin slices without 
disturbing the frozen stress pattern, so that in each slice 
the rotational effect is negligible (to freeze the stresses 
into it, the model is made to undergo a special thermal 
cycle) and (li) the scattered light methods where 

measurements on light scattered in a plane perpendicular 

to the light Path are made. 

17 

In 1939 Weller made use of the scattered light 

phenomenon in his analysis of three-dimensional state of 

stress in a photoelastic model. Using the ©ther theory of 

18 

light Drucker and Mindlin investigated the effect of 
rotation of secondary principal-stress axes along the wave 
normal, on the retardation in a photoelastic model. They 
solved the problem where the principal-stress magnitude 
remains constant and their directions change linearly along 
the wavenormal. They concluded that (i) the rotation incre- 
ased the relative retardation by a factor S = 
where R is the ratio of the rate of rotation to the rate of 
retardation and that (ii) at high stress the light vector 
rotated with the principal-stress directions, plane polarized 
light entering along a direction of princiiHl-stress remained 
practically pH^ne polarized in the direction of stress 
through the entire model. 

Mindlin and Goodman derived optical equations , 
of three-dimensional photoelasticity from Maxwell's electro- 
magnetic equations by neglecting the Inplane (plane perpendi- 
cular to the wavenormal) variations of the indeces of 
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refraction and the displacements and some other quantities 
of very low magnitude. They derived Maxwell -Neumann equa- 
tions^^ and showed the order of magnitudes of necessary 
approximations involved. They also showed that the Para- 
meter R was the essential one controlling optical wave pro- 
pagation in tliree -dimensional photoelasticity. 

PI 

Ginsburg derived yet another system of differen- 
tial equations governing the propagation of electromagnetic 

Op 

waves in an anisotropic, inhomogeneous medium. Aden 
simplified these equations by neglecting certain quantities 
and assuming a particular form for the solution. His equ- 
ations are finally reduced to two first order ordinary diffe- 
rential equations with variable coefficients. He proved 
the equivalence of any length of optical path in the model 
to a system containing three parameters called the charac- 
teristic parameters in analogy with the parameters for a 

23 

system of birefringent plates. He suggested that the 
experimental data of these parameters for various wavelengths 
of light would yield additional information regarding the 
stress distribution inside the photoelastic model. 

Frocht and Srinath developed a scattered 

light method in which they assumed that; when the amplitudes 
of light vibration along the secondary principal-stress 
axes at the entrance were equal, T*rith the addition of re- 
tardation along these axes at the entrance (i) the amplitudes 
along the secondary principal-stress axes at any point on 
the light Path did not vary and (2) whatever retardation 
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was added at the entrance, the same amoimt appeared as an 
additional retardation along the secondary principal-stress 
axes at any point on the optical p:ith. This is equivalent 
to ass-uming that the retardation vs. distance (along the 
wavenormal) curve shifts parallel to itself when a certain 
amount of retardation is added at the entrance. These as sum 
ptions however are valid for the cases where the rate of ro- 
tation is not high, and these are based on the findings of 
Drucker and Mindlin^® . 

Aderholdt, Mckinney, Ranson and Swinson^^ in their 
scattered light method treated the light components to he 
following the secondary principal -stress axes in cases where 
the rotation to retardation ratio is small. This was also 
based on the findings of Drucker and Mindlin^S . 

Robert and Guillemet^^»^^>^® followed an altogether 

different approach to overcome the rotational effect . They 

used the ellipticity measurement technique that was common 

23 

in crystallographic investigations , In this technique, 
the analyser is rotated at a particular speed and the ratio 
of the amplitudes of alternating and direct components of 
light from the analyser is measured with suitable electronic 
instrumentation. The analysis in this method is not direct 
but throiagh the characteristic parameters. This method 
however, does not involve any assumptions as in the above 
methods. 



8 


Present Investigation 

The objective of the present investigation is 

(i) to develop the necessary equations describing 
optical phenomena in a photoelastic medium in the 

-presence of rotation of the secondary principal-stress 

pq 

axes, in terms of Stokes parameters , 

(ii) to solve these equations analytically or numerically 
for various cases and study the validity of assump- 
tions of previous investigations, 

(iii) to use these equations for proving the existence 
of an optically equivalent system consisting of a 
retarder and a rotator, 

(iv) to devise methods (both with transmitted and 
scattered light) for the determination of the 
parameters of the optically equivalent system 
(characteristic parameters) experimentally, 

(v) to devise new graphical methods to describe 
passage of light through a photoelastic model in 
which the rotation of secondary principal-stress 
axes is present and to apply these graphical methods 
for the determination of the integrated optical 
effect and the characteristic parameters, 

(vi) to develop new techniques of compensation 

(vii) to develop methods for the determination of integral 
fringe orders (upto any value, when there is no 
stress free comer to count from), 
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(■viii) to devise and perform experiments to check the 
validity of (a) the governing optical equations 
(h) the existence of an optically eqiiivalent system 
and the techniques for the determination of the 
parameters of this equivalent system (c) the assump- 
tions of the previous scattered light methods and 
(d) to find the variation of characteristic para- 
meters with the wavelength of light. 

A brief outline of the contents of each chapter is 
presented below. 

In Chapter 2, a discussion on the available systems 
of optical equations in three-dimensional photoelasticity 
is presented. A new set of ordinary, linear simultaneous 
equations are derived kinematically in terms of Stokes para- 
meters . These .are very convenient for an, 3 ,lysis and application. 
Using these equations, the analytical solution for the case 

where the ratio R of the rate of rotation (0' ) to the rate 
of retardation p' is a constant (0’ can be any arbitrary 
function) is derived. Also presented is a discussion on the 
concept and properties of light ellipse. 

In Chapter 3j the existence of an optically equiva- 
lent model for any length of light path in the photoelastic 
model is proved. The importance of the C'Oncept of an opti- 
cally equivalent system is appraised. Some special cases 
in which there exists some relation between the parameters 
of the optically equivalent system (called characteristic 
parameters) are discussed. Two direct and simple methods of 
determining these parameters experimentally are described. 
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In Chapter 4, a discussion on the previous methods 
of scattered light analysis is presented. The validity of 
the assimptions involved in these methods is investigated 
using the analytical solution derived in Chapter 2, the 
nimerical solutions for some cases of loading ( of a rectan- 
gular har under combined tension and torsion) and the concept 
of 'optically equivalent system'. Modifications of the 
methods of determining the characteristic parameters with 
transmitted light (described in Chapter 3 ) Tor use with 
scattered light are suggested. 

Chapter 5 on Miscellaneous Topics 'deals with the 
following; . 

(a) 'Discussi'-n on various graphical methods in three - 
-dimensional photoelasticity' - In this section a 
brief discussion of various existing graphical 
methods developed by Jerrard ^ , Robert and Guille- 
liiot , Xuske , Schwieger26 ^ Cemosek27 ^ is pre- 
sented. A new, direct approach based on the concept 
of light ellipse with two graphical methods is des- 
cribed. Use of these graphical methods in finding 
the characteristic parameters is demonstrated. 

(B) 'Discussion of the existing methods of compensation 

and development of a large field compensator' - In 
this section use of (a) a system of four retarders 
each having the same retardation, with the pola- 
rizer or analyser and (b) a system of eight retar- 
ders each having the same retardation, as a 
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compensator is described. This is essentially an 

OQ 

extension of the works of Pancharatman and 
Tuzi ^2 

(C) ’Deteminaticn of integral fringe orders in . 

photoelasticity' - In this section the construction 
of a nomogram and a ready reckoner table which can 
be used to determine the integral fringe order 
(virtually upto any value) from the fractional 
fringe orders for two wavelengths of light, is des- 
cribed. This method is essentially an extension of 

30 

the method developed by Pant . Some experimental 
values (for a circular disc under diametral compre- 
ssion) are presented to demonstrate the value of 
these methods. 

In Chapter 6 are presented results from experiments 
conducted on a rectangular bar under combined tension and 
torsion. These form the supporting evidence for the validity 
of the governing optical equations and the existence of an 
optically equivalent system. The techniques described in 
Chapter 3 are successfully employed to determine the chara- 
cteristic parameters experimentally for various values of 
tension and torsion. These experimental values are compared 
with the values obtained by numerically integrating the 
governing optical equations derived in Chapter 2. The ' 
agreement between these two sets of values is found to be 
good. The characteristic parameters are also found experi- 
mentally for three wavelengths of light for the same loading. 
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The change in the characteristic directions is when the 
v^avelength is- changed from 500 *+ to 5950 °A (the correspond- 
ing change' in the material fringe constant is from 128 - 158 
Psi. in./fi^irge) . For the particular values of tensile and 
torsional loads considered, the assumptions of the previous 
scattered light methods are not found to be valid. 




CHAPTER 2 


THE uIGHT ELLIPSE MD OPTICAL EQU.^IONS 
IN TERMS OF STOKES PARAiMETSRS 

Introduction ; 

The photoelastic analysis consists of observing 
the changes in the properties of the polarized light in the 
photoelastic model experimentally and relating these to the 
(secondary) principal-stresses and their directions using the 
stress-optic law. This is usually done by resolving trigo- 
nometrically the light vector along different directions and 
examining the amplitudes and phase angles of each component. 

A simpler and more convenient way is to study the properties 
of the light ellipse and its transformations as the light 
passes through a model or a series of optical elements and 
then to use this information in the photoelastic analysis* 

Light Ellipse 

The most general type of polarized light is the 
elliptically polarized light which can be represented by 
an ellipse called the light ellipse. This can be represented 
by two simple harmonic vibrations along two orthogonal refe- 
rence axes OX and 01, Fig. 1., as; 
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Ex = ax Cosiirt 

and (2.1) 

Ey = SLj Cos(“ t +0) 

where 0 is the relative phase difference between the two 
vibrations, and ax, ay are the' amplitudes along these two 
axes. 


In the case of a stressed phot oel istic medium, it 
is conventicnal to make the reference axes coincide with the 
secondary principal-stress axes. 


Eliminating wt between the above two equations, the 
equation of the light ellipse is obtained as 




+ 


2 


a _ 2 % % ^ 


( 2 . 2 ) 


The major axis of the light ellipse makes, an angle 
with the X-axis such that 
2 av ay 

tan2=< = “2 ^ose 


= tan2 e Cos 0 


(2.3) 


where tane= = amplitude ratio 


(2.i+) 


o( is called the azimuth of the light ellipse with respect 
to the reference axis -CX. 

If a and b are the lengths of the semi -major and 
semi -minor axes of the light ellipse, their values are 
given by 





17 


pop 2 P 

a = af Cos °c + ay Sin c< + aY CoSo( Sino( Cos0 

P 2 2 2 2 (2*5^) 

b'^ = a^ Sin o( + ay Cos c< - 2 ax ay Cosoc Sino( CosG 
2 2 2 2 

so that a + b = ax + ay (2*6) 

The ratio ; 

tan a3=^=e (2.7) 


e is called the ellipticity. 

The light ellipse can also be referred to its own 
major and minor axes ?and n- giving 
Ejr = a Cos(“t +Co) 

^ ( 2 . 8 ) 


= b Cos (“it +£0+^11: ■^ 2 ) 


= + b Sin(wt + Gq) 

aY SinG Sines^ 

Where tanC^ = Oos=c + ay Sin=< Cos£ 


( 2 . 9 ) 


The plus or minus sign indicates left handedness 
or right handedness of the light ellipse. 

The following points are observed and used in the 
subsequent analysis. 

1 ) When the light ellipse is referred to its own axes, the 
relative phase difference is ± v/2* 

2) If E^i is ahead of E^ by ■b'/ 2, the light ellipse is right 

handed. It is ahead of by r/2, the ellipse is left 
handed, Fig. 2. ; 

3 ) If a phase difference of ±ir /2 is added along the major 
or minor axis of the light ellipse, it becomes a line ellipse 
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(linearly polarized), oriented at an angle o) such that 
tan w = + b/a. 

4) If the amplitudes and slj along the reference axes 
are equal, then tan2o( =oo , which means, that the axis of the 
ellipse will be at 45° to the X-axis. 

5) If a;^ = ay and G = Tr/2, the light ellipse becomes a 
circle and a = b (circularly polarized - from Eq. (2.5))- 

The state of polarization as represented by Eqs. 
(2.5) and (2.6) can be expressed in terms of the azimuth oc 
and the ellipticity b/a. These two quantities are sufficient 
to express completely the state of polarization if one dis- 
regards the absolute intensity (as is the case in most of 
photoelastic work). 

Stokes Vector 

The components of light vibration along the refe- 
rence axes OX and OX, Fig, 1, are given by Ej^ and E^ in 
Eq. (2.1). The components of this vibration along a set of 
new axes CX’ and OX’ , Fig . 3, inclined at an angle ofejf* to 
the OX, OX axes are given by 

Ev' = a.y Cos^^t Cos oc’ + av Sinef^ (CosoTt Cos0 - Sin“t Sing) 

^ X (2.10) 

Eyi = - Co Sot Since’ + ay Cosc(* (Cosot CosG - Sinot Sing) 

The amplitudes along the new axes are given by 

a^t = (a;j£ Cos3(» + a^ SiDof’ Cosg)^ + (ay Sirb(' Sing)^ 

a|.t = (- ax SiTtec’ + ay Cose(’ Cosg)^ + (ay Coso(* Sing)^ 


( 2 . 11 ) 
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. . a^! + a^i = + a^ since the total intensity will not 

change with the rotation of the reference axes. 

From Eqs. (2.5’) and using Eq. (2.3), 


a^ - = (a^ - a|-) Cos2cC + Sin2o( Cos£ 

= (a| - ay) Cos2oc + (a^ - a§-) tan2o( Sin2o( 

2 2 
a y - af 

~ Cos2=<: 

But, from Eqs. (2.4-) and (2.7) 

a^- b^ _ Cos2a) 

“ C»s2e 

. . Cos2wCos2^ = Cos2g (2.12) 

Multiplying this by Eq. (2.3), we obtain: 

Cos2 £0 Sin2o( = Sin20 CosG (2.13) 

Squaring and adding the above two Eqs., 

Cos^2 o = Cos^2e + Sin^2eCos^G 

or Sin^2 a ~ Sin^20 Sir^G 

Taking the positive square root on both sides, 

Sin2 w = Sin2e SinG (2.14-) 

The four quantities ^y)j ”* ^ 

32 33 34- 2 3. 

and 2a^ Sin£ are referred to as Stokes parameters'. ^ ’■ ’ 
which were introduced by the British fhysicist G.G. Stokes 
in 1852. The first quantity (a| a|) corresponds to the 

total intensity of light vibration. The variations of this 
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total intensity due to absorption in the medium are not of 
any consequence in most of the photo elastic investigations. 
For convenience, we set this quantity equal to unity. We 
shall represent the Stokes vector by V = (Sq S-j ^82^82). 


Throughout this discussion only the last three 

parameters will be used ignoring the first one as it will 

not affect the present analysis. The Stokes vector will 

;si'' 

either be -rritten in the form of a column vector 


s; 


or in 


the form of a row vector (8^ 8^ S_) depending on the convene 

' » ‘^1- i 

ience. 


^ Now, we have the following equivalent forms for 

the Stokes vector (with a^ + a|. = 1.)’, 




2 2 

- 4 


Cos2#’ 


Cos2w 

c\J 

CO 

0 

0 

Sq 

=: 

2a^ ay CosG 


Sin2e CosG 


3 

OJ 

CO 

0 

0 

Sin2o( 


I 

2a^ ay SinG 


Sin2e Sine 


Sin2w 



It may now be noted that these parameters are with 
reference to a set of perpendicular axes which make an angle 
ofb(with the major and- minor axes of the light ellipse". If 
these reference axes are rotated, the quantities £, o(^ 
a^ and aj change and so the Stokes parameters will also 
change cor re spending ly. 

Also sf + s| + sf = 1 (2.16) 

Since 8 -j = - a|., the first Stokes parameter S^ 
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depends only on the amplitudes of components along the 
reference axes. Also, = SinPo) and the third Stokes 
parameter depends only on the elliptic it y of the light 
ellipse . 

Transformations of Stokes Vector 

When a certain amount of retardation is added 
along one of the reference axes, say, OY-axis, Fig. 3, the 
phase differences , becomes (£ +£-^), and the amplitudes 
ay; and ay remain unaltered (i.e. 6 remains unaltered), i.e. 
the parameter $2 = Sin2 ©Cos £ becomes S 2 = Sin2 ©Cos( £+ £:i 
and S^ = Sin2 ©Sin € becomes S^ = Sin2 ©Sin( £+ £ 3 ) 



Si 




h 

i.e. the vector 

S2 

becomes 

^2 


S2 Cos £1 - S3 Sin £■_ 


S3 


% 


S^ Cos £1 + S2 Sin £;. 

I 





^ (0 . 


If the reference axes are rotated in the anti- 
-clockwise direction through an angle '^ 3 , Fig. 3, the azimuth 
<=< becomes ( =< - =< 1 ) and the ellipticity remains unaltered, 
so that S-j = Cos2 wCos2 a becomes S^ = Cos2“ Cos 2 (o( - oC 2 _) 
and S 2 = Cos2a) Sin2c< becomes S 2 = Cos 2 wSin2(oc - c< 2 _) • 
i.e. the vector 


-sf 


S^ Cos2°fi + S 2 Sin2oc5 


X 

^2 

becomes 

$2 Cos2=<3_ - Si Sin2®(iL 


^2 


■ i 

^3 

j 

"^3 


( 2 . 18 ) 
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Now, consider the addition of a retardation of £2_ 
along the set of axes OX', OY', Pig. 3, which make an angle 
of the original reference axes. For this process, 

the transformation is obtained by first decreasing the 
azimuth =< by an amount without altering the ellipticity, 
and then increasing the rh.aso difference by an amount 
without altering tile auiptitudes. This is equivalent to 
the successive application of t-he above two transformations. 
In this case, the vector (S-]jS2,S3 )becomes, 

"Si Cos2 "=<1+ S2 Sin2 fe'i~ 

(S2 cos2 Sin2°<x)Cos S3 Sin % = S2 (2.19) 

S3 Cos £]_+ Sin €3(82 Cos2'<i- S-] Sin2=<i) S3 

Here, the vector ("§-1^32^33 )is with respect to the 
new set of axes OX', OY'. If this vector is again referred 
to OX, OY axes, we will have 




(S'! Cos2<=(2+ 32 Sin2<=(2)Cos2oC2_ - Sin2oC]_^S2 Cos2oc ^ 



- S-, Sin2of^)Cos £]_ - S3 3in £]_ 

CM 

m 

=: 

|(S2 Cos2o(^ - S-| Sin2c<j^)Cos 0^ - S3 Sin0jcos2cc^ 



+ Sin2c5C2^(S-| Cos2oCj_ + S2 Sin2c(2_) 



S3 CosCj. + Sin 0^(82 Cos2o(^ - S-| Sin2cCq) 


( 2 . 20 ) 


Intensity Calibration Method To Determine Ellipticity 


Consider a system consisting of a x/4- plate and 
an analyser, the angle between their axes being 4-5’°. Let 
the Stokes vector for the light ellipse before this system 
with respect to the axes of the plate be A-^ ,= (31^82^33). 
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Let B.| be the Stokes vector for the light ellipse after 
the ^/+ plate referred to the (^A) plate axes and C-^ , the 
Stokes vector of the light ellipse after the (^A) plate 
referred to the axes of the analyser. Then the transfor- 
mation through this system is given by 




B- 


Si 

■ 

oi 




~Si' 

S2 


-S3 


-Si 

= 

S2 

.®3. 




S 2 


-^3 


( 2 , 21 ) 


Using the transformation laws given by Eqs. (2.17) 
and (2.18). If 




— 2 2 
and = a^ - ay 


sf 


Cos2t^ 0032*^ 

S 2 


Cos2u Sin2 =C 

.^3. 


Sin2w 


then, S-] = - Sin2w = 2a£ - 1 
or a| . 


Since 


a| + a^ = 1 


Also, the intensity along the OX -axis is given by (OX,OY 
now refer to the analyser axes) 

I=Ka^=-2(1 - Sin2io) (2.22) 

where K is the constant of proportionality. 

So, we have the following theorem; 

Theorem 1 ; Consider a system consisting of a ( >A) Plate 
and an analyser, the angle between their axes being 45°. 
The intensity of light I from the analyser is dependent 
only CO the ellipticity of the light incident on the 




5 “ Source of light 
P ” Polarizer 

(AM),,(V4)j- First and second quarter-wave plates. 

M ~ Model 
P “■ Fast axis 
S - Slow axis 

P,q~ Principal -stress directions 

A — Analyser 



OPTICAL SET-UP FOR TWO 
STRESS ANALYSIS. 


- DIMENSIONAL 
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(^A) plate and is given by Eq. (2.22). This theorem can 
be used with advantage in two-dimensional stress analysis. 
Ccnsider the optical set up shoxm in Fig. tj., The light 
after the (>./4-)'i plate is circularly polarized so that the 
amp'iLtudes along the p,q-axes are equal. Along p,q-axes 
the phase difference is Tr/2 before the model and (ir/2 + p) 
after the model, if the amount of retardation added in the 
model is p . From the property of light ellipse, the ellip- 
ticity of light after the model is tan (v/if + p/2). Using 
Theorem 1 , the intensity I from the analyser is given by 
Eq. (2.22), where tan 03 is the ellipticity. 

S o 

1=^(1- Sin (Tr/2 + p ) = ^ (1 - Cos p) (2.23) 

This equation relates the retardation introduced 
in the model which is proportional to the principal-stress 
difference (p - q) to the intensity I from the analyser. 

To find the constant of proportionality K, Tardy's method 
can be used. For this, remove the (^A )2 Pl^-^e and rotate 
the analyser for minimum intensity, thus obtaining the 
directions of the major and minor axes of the light ellipse 
after the model. Now set the axes of the (^A )2 Plate 
parallel to these major and minor axes so that the light 
after the (^A )2 Plate is plane polarized. Rotate the 
analyser for minimum intensity, thus finding the direction 
of plane polarization. The angle between this direction 
and the major axis will be w where tanw is ellipticity of the 
light ellipse after the model and to = ( p + ir/2)/2. Now 


2 ? 

arrange the (^/4)2 plate and analyser as shown in Fig. and 
note the intensity I. Since the quantities I and p are knom, 
the constant K can be evaluated using Sq. (2.23). Once the 
value of K is kno\Nni, the principal-stress difference can be 
computed using Eq. (2.23) from the value of + for any other 
point on the model. A calibration graph showing the relation 
between I and (p - q) can be drawn for convenience. 

Optical Equations of Three-Dimensional Photoelasticity 

Several systems of equations have been put forth 
to describe the optical phenomena in the presence of the 
rotation of secondary principal-stress axes. We shall give 
a brief account of some of these derivations. The first of 
these are the lyfeixwell -Neumann equations; 
da = - a Cos d0; da = - a Cos d0 (2.2^) 

and d = (-^ - •—) Sin d0 + dp (2.25) 

where d0 is the change in the orientation of the secondary 
principal -stresses and dp is the rdtardation introduced 

(due to stresses only) in an infinitesimal length dz of the 

20 

optical path in the model. Coker and Filon derived these 
equations kinematically. They assumed the components of the 
light vector in the trigonometrical form. Also, in this de- 
rivation, after a distance dz along the optical path in the 
model, the reference axes were first rotated to coincide with 
the directions of the new principal-stress axes and then the 
retardation caused due to stresses alone was added to obtain 
the components of the new light vector. 



28 


18 ‘ 

Drucker and Mindlin"f ormulated the problem in terms 
of the elastic ether theory of light. The ether theory is 
based upon the ideas of stress ^ strain and displacement 
and upon the equations of motion. The governing differen- 
tial equations are exactly the same as those obtained from 
the electromagnetic theory of light if the symbols repre- 
senting displ'O cement, density and shear modulus are inter- 
preted as magnetic force, magnetic permiability and reci- 
procal of dielectric coefficient respectively. They obtai- 
ned the analyi:ical solution for the case, where the secondary 
principal-stress difference remained constant and their 
orientation varied linearly along the light path. They 
concluded the following ; (1) Rotation increases the^rela- 

r T ^ 

tive phase retardation by a factor S = 1 + (-^) | 

which at times may be quite large, (^is the orientation 
of the secondary principal-stresses for the case under 
consideration, a is the relative retardation caused by 
the magnitude of stress alone). If the ratio - 7 -^ is not 


constant, the ratio of the increment in to the increment 
in A between fringes can be used in regions of high stress. 
(2) At high stress, the light vector rotates with the 
principal-stress directions, plane polarized light entering 
along a direction of principal-stress will remain practi- 
cally plane polarized in the direction of the stress through 
the entire model. 
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Mlndlin and Goodman derived a set of dynamical 
equations of three-dimensional photoelasticity using the.- 
Electromagnetic theory, with the purpose of revealing the 
approximations involved in the hinematical derivation. 

For this they made an analysis of the orders of magnitudes 
of the quint it ies involved in the equations. The orders 
of magnitude of these quantities are presented below. 


(1) 

"is 

ds 

= 0 

A^ 

OS 

(ii) 

3 

1 1 

d3_ 

= 0 

, M 

^ ds 


(ill) 


"b 


r bz 


= 0 

X z m 


Jco 

OS 

L 

cs 


G3 

6s 


(iv) 1 bz 

. to d3^ 



0 




( 2 - 26 ) 


(v) 




(vl) z 

w QS 6s 

S, n and b are representative lenghts in the 
X,Y plane (the plane of the wave front), index of. refrac- 
tion and wave velocity respectively, X is the wavelength 
of light in vacuum and 0 "^stands for”order of magnitude 
of”. To estimate the order of magnitude of these ratios, 
consider a wavelength of 6x10“^mm, a stress gradient of 


AZ 


OS 



30 


10 ° dyne s/cm 2 (s 1500 Ib/inT) per nm. and a stress-optical 


coefficient of 100 Brewsters. Under these extreme conditions 
' ^ and, for a model 5 thick, 2 = 0 jlO 

With regards to terms involving p, if the rotation of the 




0 


10 


principal axes were 2 ir radians in one millimeter, 


0 


oS. 


10 




Since the ratios involve squares or 


products of terms not exceeding 0 


selves donot exceed 0 


~io-6l . 


10 


-3 


the ratios th em- 


it is true that in the 


n eight) our hood of an isotropic piont, ^ becomes very large. 

'■0 

However, may be expressed in terms of ~ and the principal 
0'S 0- S 

stress difference through the use of the stress equations 
of equilibrium and the stress-optic law, following ^rhich 
it may be shown, by a limiting process, that the ratios 


involving dp are, again, small. Also, the order of magni- 

as 

tude of (n-5 - n2)/n and x M (^1,^2 refractive 

indeces) isjlO”^. It is not impossible for^ to become 
very large. An example is the case of combined tension 
and torsion of a cylindrical bar. For a wavenormal inter- 
secting the axis of the bar^ approaches infinity as 

dz 

the tension approaches zero and the product of the two has 
a zero limit in the neighbourhood of the axis of the bar. 
They simplified the governing equations by neglecting these 
quantities. They obtained the Maxwell-Heumann equations by 
assuming the solutions in the trigonometric form. By 
changing the independent variables, they obtained the 
following simplified form of the optical equations. 
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R u'= S-?' 

9 T 

-y* = — 

9 a 


Where 


R 




(2rr/x) (n 2 - ) 

a rroT (t - JK^"' dz) 


to(t - JK 2 dz) 



-1 -2 ? -1 -2 
VI = y n2 


(2.27) 

( 2 . 28 ) 

(2.29) 

( 2 . 30 ) 

(2.31) 

( 2 . 32 ) 


In the electrostatic system of nnitsy = 
where is the velocity of light in vacuum, u' and v* 
are the components of the light vector along the secondary 
principal-stress directions. 

Eqs. ( 2 . 27 ) and (2.28) show that R is the essential 


parameter controlling optical ’'rave propagation in three- 
-dimensional photoelasticity. Thus, we see that the kine- 
matical equations for all practical purposes are accurate 
enough as can he observed from the assumptions involved in 
the dynamical derivations. 


Aber?2 derived yet another system of differential 
equations governing the optical phenomena in three- 
dimensional photoelasticity. He used, the system of equa- 
tions for the propagaticn of electromagnetic waves in an 
anisotropic inhomogeneous medium obtained by Ginsbiirg 21, 
Using the fact that the |4iotoelastic medium is cnly 
weakly anisotropic, ' he obtained the following system of 
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ordinary linear differential equations. 
dBi - Aa 

^ C' - ag) ^ (2.33) 

and 

bI + i C’ ( cf ) Bp (2.3^) 

dz dz ' 2 X 2 a 

where c’ = GG^ (2.35) 

G is the velocity of light in vacuum and G^ is the photo- 
elastic constant, a 2 _ ^2 3-^® secondary principal- 

-stresses and b!| , B^ ai’e the components of light vector 
along the direct icsis of the secondary principal-stresses, 

Equaticns Using Stokes Parameters; 

We shcall now obtain a kinematical derivation of 
the governing optical equations in the presence of the 
rotation of the secondary principal-stress axes to faci- 
litate a better understandj-ng of the optical phenomena. 

An observation of the transformation Eq. (2.1 9)> 
shows that it is linear (i.e. no products, squares or 
trigonometric functions of the Stokes parameters are invol 
ved). This prompts us to use the Stokes parameters direct- 
ly as the variables in the kinematic derivation. 

Let the Stokes vector referred to the secondary 
principal-stress axes at points 0 and o' be (S.^ , Sp ? S^) 
and (S.^ , S 2 , S 3 ) respectively. The point 0 is at a 
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distance z from the origin idiich is at the point of entrance 
to the model, Fig. 5* Let the change in the orientation of 
the secondary principal-stress axes over the distance dz be 
d0 and the retardation introduced due to stresses alone be 
dp . Then using the transformation represented by Eqs. 
( 2 . 19 ), and noting that for sufficiently small distance dz, 
the quantities dp and dp will be very small so that 
SlJi2^^2dp, Sin dp ^ dp, Cos 2 dp % 1 ,Cos dp ^ 1 , djZ! dp 0, 

we obt ain 


~Si' 


S.| + dS-i" 


S.^ + 2S2 d^ 

S2 


S2 + dS 2 

= 

S2 - 2S.J d^ - S^ dp 

S3 


S3 + dS3 


S3 + S 2 dp 


(2.36) 


Dividing each of the above equations by dz and taking the 
limit as dz 0 , we obtain, 

1 d0 

®i = 'ST = ssg ■£ = 2 S 2 P (2.37) 

s' = - 2 S. - S. - 2S/' - S. p> (2.38) 

2 dz ‘ dz 3 dz ^ 3 

and 

S3 = -^ = S2 “ = Spp' ( 2 . 39 ) 

If the Stokes vector (S-^ , Sp , S 3 ) is referred to 
the fixed axes CK, OY, we have from Eqs. (2.20) (replace 
=( .j by p + d^ and £ .] by p ) 

S.[ + dS.^"” 

Sp + dS2 
+ <383 


5 1 

5 2 
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(S^ Cos2(^ + d^) + S2 Sin2(d^ +p)) Cos2(jZi + 60) 

- B±n2(0 + 60) J(S2 Cos202( + d^) - Sin2(^ + d^)) Cos dp 

- Sin dp j 

Cos2( 0+60) 1^(82 Cos2( 0+ 6l0) - Sln2( 0+ 60) Cos dP 
= - Sin dpj + Sin2( 0+ 60) (S^ Cos2( 0+ 60) | 

+ S2 Sin2( 0+60) ; 

S3 Cos dp + Sin dp (S2 Cos2( 0+60) - S., Sin2( 0+601 

- ( 2 .^ 0 ) 

As above, for sufficiently small distance dz, 6 0 
and dfP are very small so that Sin dp ^ d , Cos: 2d^ ^ 1 , 

Cos dp^ 1, di^^dp^ 0, Cos2( 0+600 Cos2 0 -2i2d 0 Sin2 0 ; 
Sin2(djZS ^2dj2Cos2jZi + Sin2;2f so that 


S1 


S^ + dS^ 


S-] + S3 dp Sin2^ 

S2 




S2 “ S3 dP Cos 2 jZl 

3 _ 


S3 4- dS3 


S3 4- S2 dP Cos 2 jZS - S^ dp Sin 2 jZ( 



— ^ 


WMM. 


(2*Vl ) 


Dividing each of the above equations by dz and taking the 
limit as dz — »- 0, we obtain 


dS^ 

^1 - dz 


and 

s' 

®3 ■ <lz 


S3 p’ S in 2^ 

- S^’ Go32;2i 

- S 2 p' Cos2^ - S-] p* Sin2p 


(2.42) 

(2.43) 

(2.44 ) 


These equations can also be obtained from Maxwell- 
-Neumann equations by converting the dependent variables 
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into Stokes parameters. 

Eqs. (2,42) - (2,44) will be particularly suitable 
for scattered light analysis. This will be discussed in 
a later chapter. The above sets of equations are, linear, 
ordinary simultaneous first order differential equations 
which are very convenient for analysis and application. 
Using the Eqs. (2.37) - (2,39), 'we shall obtain the analyti 
cal solution for the case where the rate of rotation ^can 
be any arbitrary function, but the ratio R is a cons- 

tant, (The case studied by Drucker and Mindlin where both 
and p' are separately constants is a particular case of 
this) . 


Let k2 = 

4r^ 


frcm Eqs. (2.37) and (2.39) 

s' 

3 2R. 


(2.45) 


Integrating both sides, 

S- 


S. 


2R 


+ Ai 


where A-j is an arbitrary const .ant. 

Substituting this expression for S^ into the equatim 

P P P 
+ S2 k 

(instead of using Eq. ( 2 . 38 )) 
we obtain an expression for S in terms of S^ . 
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Substituting this expression in Eq. (2.37) and integrating, 
we obtain 


A-1 1 A-j -S' 

= S3n(2K jZ!+ Ap) (2.4-6) 

2RK‘^ 


where Ap is an arbitrary constant. 

Using Eq. (2.37) again, we obtain 

S„ = % z: (1 - Cos (2K^+ Ap) (2.47) 

and (1 - 

^ ^ = ^ + 2 b,k ^ 2 } (2.48) 

The arbitrary constants A^ and Ap are to be evalua- 
ted from the initial conditions. 


Jones Vector 


It is sometimes more convenient to use Jones 


Op 

vector in photoelastic analysis. 

This vector has two 

e lement s 


KCy + oTt) 

P = ay e 

(2.49) 

and _ 

i(G + w t) 

Q = a e ^ 

(2.50) 


The quantities ay are the amplitudes along two 

reference axes OX, OY; and Cy phase angles of 

the instantaneous components along these axes. 

i(ey + (flt) iUt i^^ 

_ Since e = e e and the term. 

i“t 

e occurs in the expressions for both P and Q this can 
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be ignored in the photoelastic analysis where the varia- 
tion with time is of no consequence. So, we have 

P = a^ (2.51 ) 

i£ 

q = 3,J e (2.52) 

It can be noted that the transformations corres- 
ponding to the rotation of the reference axes will be 
linear in P and Q. So, if we derive the optical equations 
of three-dimensional photoelasticity in terms of the ele- 
ments of the Jones vector, -re should obtain ordinary linear 
differential equations in P and Q. 

Now, referring to Fig. 5? let (P,Q) be the Jones 
vector vdth respect to the axes OU, OV at the point 0.' If 
this vector is referred to the axes OU', OV^ its elements 
b ec ome 

P' = P Cos A jZl+ Q SinAjg 
= - P Sin A 0+ Q CosAj?) 

Now, if a retardation of Ap is added to the P' 

component, the elements of the new vector are given by 

iAp iAp i!ip 

p" = P' e = P e CosA jZJ + Q e Sin A 0 

and Q*'* = q’ 

For sufficiently small distance Az, Ap and L0 

are very small so that Sin Ap ^ Ap , Sin h 0 ^0^ Cos Ap ^ 1, 

iAp 

CosA^^I, Ap A4 )^ 0, e ^ 1 + lAp . We, then have 

p" r: P + AP = P(1 + iAp ) + q 
and Q*’ = Q + A, q = - P A j? + Q 
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Dividing these equations by Az and taking the limit 


CO 

0, ■ we 

1 obtain 

S = P'=L 

dp p 

+ Q 1^ = • 

dz 

dz 

^ dz 

and 



12 - o’ - 

dz “ ^ " 

dz 

P^ 

1 

II 


(2.53) 


(2.5+) 


Instead of adding a retardation of Ap to one com- 
ponent only, if we add retardations ^ and -Ap/2 to the two 

iAp/2 2 -iAp/2 

.ompcnents (i.e. P = P e , Q =: Q e ), we obtain 

P' = - | »'P + 0' Q (2.5J) 

and q' = - 0' P + I p’ Q (2.56) 

These equations are s.ame as those derived by 

Ab en . 


In this chapter, we have introduced the concept 
•f light ellipse as an elegant tool in photoelastic analy- 
sis. V.fe have given a brief review of some of the previous 
derivations of the optical equations of three-dimensional 
photoelasticity to show that, for all practical purposes, 
the kinematically derived equations are sufficiently accurate, 
i^e have also derived two new sets of equations for three- 
-dimensional photoelasticity, which are convenient for 
analysis and application.. Using these equations, we obtained 
the analytical solution for the case- where the rate of rota- 
tion 0’ can be any arbitrary function, but the ratio 
R = 0'/p’ is a constant. The contents of this chapter form 
the basis for subsequent chapters. 



CHAPTER 3 


CHARACTERISTIC PABAMETERS IN THREE-DIMENSIONAL PHOTO- 
ELASTICITY AND THEIR EXPERIMENTAL DEI ERM ELATION 

Poincare' Sphere Representation ; 

In this chapter, we shall discuss the existence of 
an optically equivalent system, its theoretical and experi- 
mental determination and its importance in three-dimensional 
photoelastic analysis. In crystallography, it is establi- 
shed that a series of retarders is equivalent to a system 
consisting of a single retarder (an optical element which 
adds a certain amount of retardation to the components of 
the incident light along a particular set of axes) and a 
Pure rotator (an optical element which rotates the incident 
light ellipse through a certain angle kno-'Am as the rotatory 
power). This has been proved geometrically by using the 
Poincare' sphere representation of polarized forms. We shall 
present this proof below. 

An elliptically polarized light can be characterised 
by two quantities (i) the azimuth and (ii) the ellipti- 
city tan o . This can be represented by a point P on the 

f 

surface of a sphere of unit radius called the Poincare 
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sphere with co-ordinates : x = Cos 2 w Cos 2 °fj y = Cos2co Sin 2 '=<, 
and z = Sin 2 w(these three correspond to the three Stokes 
parameters S-] , $2 ^< 3 . S^ of the light ellipse), as shown 
in Fig, 6. (Let us recall that the first Stokes parameter 
S^ = Cos2a) Cos 2 =< corresponds to the light intensity along 
one of the reference ■ axes. If the amplitude of light vib- 
raticn along X-axis is and the intensity is I, then 

I = K = ^ - ) . The angle 2 o( is measured clockwise 

from H. The angle 2 .^ is measured upward from the equator,. 
The upper and lo-wer poles and represent left and 
right circularly polarized light ( 2 u) = 90°, tan u = 1 ). 

Points on the equator indicate polarized light (plane, 20 = 0 
tan 0) = 0). The action of a linear retarder can be con- 
veniently represented on Poincare' sphere. Consider an 
elliptically polarized light represented by the point P, 

Fig. 8. with ellipticity tan w and azimuth o(, to be incident 
on a linear retarder with retardance A and azimuth cc^ . 

Its axes are represented by line RR-j in the figure. The 
X co-ordinate of the, point P corresponds to the intensity 
of this light vibration along one of the reference axes. 

The projected length of the radius vector OP along the 
axis ROR,^ corresponds to the light intensity along one of 
the retarder axes. (This can be proved as follows; X - Y 
plane is shown in Fig. 7 (a). Ox', OY' are the new reference 
axes with XOx' = 25 . "P is the projection of P on the 

X - Y plane. The projection of "P on OX 'is given by, 
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x' zz X Cos2oC 4 - y Sin2 5 = Cos2 ^+82 Sin2 ^ (say). But, 
S-j is also the first Stokes parameter obtained after rota- 
tion of the reference axes through an angle 5 . i.e. the 
projection of the point P (which again is the projection of 
P on X - Y plane) on any axis QX’ in the X - Y plane corres- 
ponds to the light intensity (of light vibration represented 
by P) along the same axis CX*) . 

Since this intensity remains constant during the 
addition of retardation, the process of addition of retar- 
dation in the retarder is represented by an arc from the 
point P to P' about the axis ROR| . This arc PP' on the sphere 
makes an angle of A ai the axis ROR-j . Rotation of the 
reference axes through an angle d=< corresponds to rotation 
about the Z-axis through an angle 2d°( cn the Poincare' 
sphere. During this rotation, the z--ic©ord.i.na'te andJ so the 
ellipticity remains unaltered. 

The process of addition of retardation along the 
reference axes can be explained as follows. For this, 
consider Fig. 7(b). Point P represents a certain ellipti- 
cally polarized light referred to axes OX, OY in the figure. 
For point P, x = OM, y = 0, z = MP. Now rotate MP about 
axis CX through an angle A so that its new position is MP' . 
Now, for point P', the co-ordinates are x* = Om' = OM, y' = 

- mm' = - MP' Sin A= - z Sin A , z = MP'Cos A = z Cos A . Now 
consider the three Stokes parameters S>i^S 2 ^S 2 of the pola- 
rized light represented by the point P. These are S^ = OM 
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= X, $2 = 0, S 3 = MP = z. Consider the polarized light P' 
obtained by the addition of a retardation of a along the 
reference axes. The Stokes parameters for this are given 
by s'-^ = X, S 2 = - z Sin A and S3 = z Cos A, The same values 
are obtained by rotation about the OX -axis. Thus, the effect 
of a single retarder of retardation A can be represented on 
the Poincare' sphere by a rotaticn through an angle A 
(upwards) about the corresponding axis in the equitorial 
plane. 

Combined Effect of a Series of Retarder s; 

The combined effect of a series of retarders can be 
represented by successive rotations of the Poincare' sphere 
about the corresponding axes in the equitorial plane. This 
in turn can be replaced by a single rotation about some 
general axis in Poincar4 sphere. Alternatively, this resu- 
ltant single rotation of the sphere can be resolved into 
two rotations about perpendicular axes - the first may be 
about the Z-axis and the other will then be about some ajcis 

in the equitorial plane which may be determined by the 
23 

construction of Fig. 8 . A rotation about AO through twice 
the internal angle at A followed by a rotation about BO 
through twice the internal angle at B is equivalent to a 
rotation about CO through twice the external angle at C. 

Thus the combination of a series of linear retar- 
ders is equivalent to a system containing two elements, one 
a rotator and the other a retardation plate. Since finite 
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rotations about nonparallel axes are noncommutative opera- 
tions, it is necessary to specify the exact sequence of the 
various eleraents. Interchanging any two of them would in 
general lead to a change in the final state of polarization. 

The above result has also been proved by the matrix 

qcr 

method by Hurwitz and Jones . 

Effect of a Series of Retarders Using Stokes Vector ; 

The above result will now be proved in detail using 
the Stokes vector ignoring the variation in the total inten- 
sity. The transformation thit the Stokes vector undergoes 
when a given elliptically polarized light traverses through 
a two-dimensional photoelastic model can be obtained from 
Eq. (2.20) as; 

V = M-, V 

where 





S1 

V = 

^2 

; V = 

”32 


% 




are the initial and final Stokes vectors and is given by 

(Cos^2°c^ + 6in^2otj CosGj ) ■Sin2°^i C6s2‘=^ Sin^^ Sln2'^-j 

(1 - CosC^ ) 

Sin2q Cos2=<-j (1 - Cos€^ ) Cos^2°f>[ CosC^ - Sine, Cos2®C, 

+ Sin^2q 

- Sin2«<^ Sine^ Cos2o(^ SinC^ Cose^ 
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In the above matrix and £| are the azimuth and retardance 
of the retarder respectively. 

This means, that the effect of a two-dimen si onal 
model at a particular point can be represented by a 3 x 3 
matrix. If the light represented by the vector V is inci- 
dent on n successive models, the out coming light vector is 
given by 


V = V = MV (3.1) 

where M = . . . .M 3 M 2 M-, 

This means that the effect of a series of two- 
-dimensional photoelastic models can also be represented 
by a 3 X 3 matrix. 

Now consider the relation (3.1) 


m-^ -| m-^ 2 3 


®11 ^21 ®31 


M = ; if = ”^2 “32 

®31 “^32 “^33 ®13 “^23 “33 


Now V = ( 1^,... M 3 M 2 M -1 )“*'(?) 
= (M^ ....) (V) 


-1 T 


One can show from the matrix that M^ = 
Similarly =: ]y^ ; M 3 ''=: etc. Hence 


V = (M^ .... ) (V) 


(3.2) 


= (Mj^ M^_^....M2M^)^ (V) 


l.e. V = f 


(3.3) 



FIG- 9. REPRESENTATION OF THE AXES OF 
THE BIREFRINGENT PLATES IN THE 
TWO EQUIVALENT SYSTEMS. 
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= S 3 = 0, S 2 = 1 ; mf'] + m22 + m |3 = 1 (3.16) 

_ __ _ 2 2 2 

S'] = S2 = 0 , S3 = 1 ; 1113-1 + 11132 + 11133 ~ (3.17) 


Thus, we have six relations connecting the nine 
quantities mij? i,0 = 1>2,3 so that the matrix M has only 
three independent quantities. So, a series of two-dimen- 
sional photoelastic models (for a particular light path) or 
birefringent plates, in general is equivalent to an optica- 
lly equivalent system with three independent parameters. A 
single birefringent plate has two parameters namely, its 
azimuth =< and its retardation p . A pure rotator which 
rotates the incident light ellipse has one parameter 
namely its rotatory power. The most simple optical systems 
consisting of only birefringent plates and pure rotators 
and having three parameters are (i) a system with a bire- 
fringent plate and a pure rotator in that order and (ii) a 
system with a pure rotator and a birefringent plate in that 
order. 

If the systems (i) and (ii) are to be equivalent, 
the pure rotatory power in both systems must be same, the 
angle between the axes of the birefringent plates in the 
systems (i) and (ii), Fig. 9 . must be equal to this rotatory 
power and the retardations of the plates in both systems 


must be equal. This can be proved as follows. For the 
first system, let the rotatory power be , the azimuth of 


the birefringent plate be p and its birefringence be 5. 

For the second sys4««»y-^Gt--*th‘e^"'"Mahtities be , B* amd d’ 

r^j. T,„ KAHPU'R, I ^ ' 

I CENl'BAt U BRARt I 

'TT23g"7i 
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respectively. Now consider polarized light to be incident 
along one of the axes of the birefringent plate in the 
first system, so that after this plate, it still remains 
plane polarized along the same direction. The pure rotator 
rotates it through an angle of so that at the exit of 
the first system ^^re have plane polarized light at an azimuth 
of ( P + ) . In the second system also, let the incident 

light be polarized at an azimuth of P . The pure rotator 
makes its azimuth ( P+ 4'^). This light should be polarized 
along one of the axes of the birefringent plate in the 
second system so that we have plane polarized light at the 
exit of the second system also with the same azimuth as in 
the first system. Therefore, p + = P + or = \jj' , 

That is, the rotatory powers in both the systems are same 
and the difference between the azimuths of the birefringent 
plates in the two systems is equal to the rotatory power of 
either system. Now consider circularly polarized light at 
the entrance of both the s/stems. For the ellipticities of 
the light ellipses at the exit of these systems to be equal, 
the birefringence of the plates in both the systems must be 
equal. That is, 5 = 

Optically Equivalent System for a Gintlnuous Body : 

Now, we shall discuss the existence of an optically 
equivalent system for continuous systems. Aben^^roved the 
existence of an optically equivalent system (consisting of 
a retarder and a rotator) for an optical path in a three- 



-dimensional photoelastic model alcog which the secondary 
principal-stresses and their directions vary continuously. 
Since the governing differential Eqs. (2.3'+) , (2,35) are 
linear, ordinary simultaneous differential equationsy the 
light propagation is governed again by equations of the 
form (3.^)-(3.6). The rest of the proof is similar to the 
one presented by Hurwitz and Jones Aben called the two 

axes of the retarder as primary characteristic directions 
and its retardation as the characteristic retardation. He 
also called the mutually perpendicular directions which 
made an angle equal to the pure rotatory power with these 
axes as the secondary characteristic directions. These 
three quantities namely, the primary characteristic direc- 
tions, the secondary characteristic directions and the 
characteristic retardation, are referred to as the charac- 
teristic parameters. His proof for the existence of these 
parameters is based on the theory of -unitary matrices. It 
may now be noted that, plane polarized light entering a 
photoelastic model along one of the primary characteristic 
directions at the entrance comes out as plane polarized 
light at the exit along one of the secondary characteristic 
directions. The corresponding primary and secondary charac- 
teristic directions are termed the conjugate characteristic 
directions. The existence of these characteristic para- 
meters in the case of continuous systems will now be proved 
in detail using the optical equations in terms of Stokes 
parameters. Consider Eqs. (2.37) and (2.38), we have; 
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r 1 




20' 


= S2 = - 23^0' - S3P’ 


The primes denote differentiation with respect to 
the distance co-ordinate z. 




,20'. 


+ 2S-10' 


■i' ' 

p 


" ^2 " 20 ' 


using Eqs. ( 2 * 39 ) ( 2 . 37 )* 

On rearranging various terms and simplifying, we obtain 


S-i - S" 


20" _p" 

0' p' 


+ S. 


^ - p'2 


0 


0 ' 


0 - 


+ 


4 S 


1 


0 '’ 0 ' - 

P' - 


= 0 


(3.18) 


This is a third order ordinary linear differential 
equation i.n S..J with variable coefficients. The initial 
conditions are; 

For a given initial set of (S-|, 82)33) i.e. for a known type 
of incident light, the following quantities can be deter- 
mined at z = 0 since the values 0 and p are also assumed 
to be kno'.m; 

rr 2820' 3.nd 8" = 2820" + 20,’ 82 = 2820" 

- 20' ( 28 -^ 0 ' - S3P’) 3-^6 also known. 

The solution for this initial value problem exists 
and is unique over the closed interval in vh.ich the co- 
efficients in Eq, ( 3 . 16 ) are continuous. The solution will 
be of the form 


(3.19) 


S-j = A-] f..| (z) + A 2 ^ 2 ^^^ "^3 ^3^^^ 


(3.20) 
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where ^ and are arbitrary constants which are to be 
determined from the initial conditions given by Eqs.(3.l9). 


i . e . S-] 


■ Z=:0 


= A-j f-j (z) 


z=0 


+ C 2 ) 


’ z=0 




= A-if^’Cz) 

Z=rO 


Z=0 


+ A2f2(2) + A3f3(2) =2820} 

2=0 Z=0 'Z=i 


— n u 

S-^ = A*^ f ( 2 ) 

z=0 . 


2=0 


+ A2 f 2 ^ ^ ^ 


2=0 




2=0 


vn 




20'(- 28.0 -SoP’)j 

J -7-1 


z=0 


From these relations, the constants A^ ^ A 2 and 
can be obtained as linear functions of the initial Stokes 
parameters* Aslo at exit, we have 



(28-10 + 82 ) 

and 80 = - frcmi Eqs. (2.81 ) and ( 2 . 82 ) 

p' 

where 0 ’, p' , are the values at exit. 

Hence, "Si , Sq and are linear functions of the initial 
Stokes parameters so that we can write 

( 3.21 ) 

where M is a 3x3 matrix. 

The elements of the matrix M are not functions 
of the initial Stokes parameters S-) , 82 and S 3 . They are 


8-1 

82 

% 


M 


^3 
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however functions of 0' and p’ and also of the wavelength 
of the light that is used. Since Eq. (3*21) is of the same 
form as Eqs. (3*^) - (3 •6), all the conclusions regarding 
the existence of an optically equivalent system and the 
characteristic parameters follow as derived in the previous 
section. Vie also note that these parameters are dependent 
on the wavelength of light that is used since the elements 
of the matrix M depend upon the wavelength of light. 

Some Special Cases; 

For some particular cases of stress distrihuticn 
in the photoelastic model, there exist some relations bet- 
ween the characteristic parameters. In such cases, the 
experimental determination of these parameters is direct 
and simple. 

Case (a)„ Pure Rotator; 

We shall now show that a system of two half-wave 
plates with the angle between their axes being consti- 
tutes a pure rotator. If V is the Stokes vector at the 
entrance of this system referred to the axes of the first 
half-wave plate, the transformations through this system 
are given by 


C os2 w Cos2 c< 


C os2w Cos2c( 


~Cos2to Cos(-2'<) 

Cos2a3 Sin2®f 

becomes 

- €os2u Sin2°f 


Cos2to S±n(“-2'<) 

Sin2 oj 


■- Sin2 to 

, -irt- 


- Sin2 (0 


V 
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*2 


3 

Gos2ojCos2 (-'=<') -=<) 


Gos2wCos2 ( of + ofl^ 

Cos2u Sin2 -°0 ; 


Gos2(‘JSin2 ( =< ) 

- Sin2w 


Sin2 CO 



Cos2<u Cos2 (oC + 2o(]_) 


Cos2to Sin2 (=< + 2o(j_) 


Sin2a) 



Where is the Stokes vector after the first 
half-wave plate with respect to its axes, tan w is the 
elliptic it y and ^ is the azimuth vdth respect to the axes 
of the first plate, of the incident light ellipse, Y 2 is 
the Stokes vector of light after the first half-wave plate 
with respect to the axes of the second half-wave plate, 
is the Stokes vector of light after the second half- 
-wave plate referred to its axes and is the Stokes vector 
after the second half-wave plate referred to the axes of the 
first half-wave plate. 

This shows that the effect of the combination of 
two half-wave plates is to rotate the light ellipse through 
twice the angle between the corresponding axes of the half- 
-wave plates and in the same direction as the angle through 
which the first half-wave plate is to be rotated so that 
its axes become parallel to the axes of the seccnd half- 
-wave plate . 
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Case (b). Reversal of Light Path; 


For studying the optical phenomena when the light 

32 

path is reversed, the Jones vector is more convenient 
since it contains only two elements. 


The Jones vector V is given by V 


'xT 

y 




The phase angle for the Y component is taken to 
be aero since we are only concerned ■'-.n.th the relative phase 
difference and' not the absolute phase angle in our present 
investigation. The transformation of the Jones vector 
through a retarder with retardation azimuth 

given by 


V 


Cos=< - Sinc< 
Sin=< Cos°C 

ie. 


iei 


e 

o 


0 

1 


Coso( Sino( 

- Sinoc Cos of 


H 


e ^ Cos. ^ + Sin^oC^ Sin2'^-] _ -j 


) 


Sin2=<i 


iG-] 

(e - 1) 


i£i 


fh£d_(e '’ + 1 ) 
2 


[v] = H,V 


(Both the Jones vectors V and V are referred to 
the same reference axes). 


T 

So, M-| is a symmetric matrix, (i.e. ) 

It is of the form; 
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Let the matrix corresponding to another retarder 
be given by 




3-2 b2 

b2 C2 


So, ^2 = (M-i 

Let be the matrix representing a third retarder. 

= M3 


T T T T 

Ncm, M2 = Rj 1 ^ M2 M, ) 


For n retarders, we have 


M-, M 2 M. 




Thus, the reversal of the light path is equivalent 
to transposing the matrix of the series of retarders. 

It is proved earlier, that a system of linear 
retarders is equivalent to a system containing a single 
retarder (of retardation € and azimuth =< ) and a pure rota- 
tor (of rotatory power The transformation matrix for 

such a system is given by 


Cos'^^ 

Sin'=f'i 


- Sin'<-| 
C os^f-j 


where c<^ = C5( + t|j^ 
iC 

e Gos®< Cos°(^ 
__ + SinpC Sinc<-j 

le 

: e Cospf SincC^ 


- Sino( C 03^*1 


ie 


r 

e 0 


Cosoc Since j 

0 1 


Since Cosce 1 



J( 


ie 

e Sin<^ Cosoc^ 
- Sinoc-| C oso( ' 

ie 

e Since Sino(^ 
+ Cog=( Cosof^ 
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The rotatory power for this system is given by 
= °C-] - °c 

For a system consisting of a pure rotator (rotatory power 
) and a linear retarder ( azimuth , retardation G ) in 
that order, the transformation matrix is again M-] * That 
is, the transformation matrix can represent either one of 
the above two systems. Now, the series of retarders M-j , 

M2 , M25....M^ is equivalent to an optically equivalent 
system with a retarder and a rotator. Let this be repre- 
sented by the matrix Kj * If the light path is reversed, 
the transformation matrix of the equivalent system x^ill be 
the transpose of the matrix i.e. m|'. Now, represents 
a system with a pure rotator of rotatory power - ^ = oc - oc^ 
and a retarder of retardation and azimuth of , 

Case (c). Reflection Polariscope; 

In a reflection polariscope, the characteristic 
retardation p is added along the primary characteristic 
directions and the light ellipse. is rotated through an angle 
for the forward ray of light. When this is reflected, 
it is rotated in the opposite direction through the same 
angle \|) and a retardation of p is again added along the same 
directions. Thus the combined optical path (forward and 
backward) is equivalent to a single retarder, (In this 
case, if the forward optical path is equivalent to a system 
represented by matrix M, the backward optical path will be 
equivalent to a system represented by the matrix M^.) 




FfG.10. DlSTRiBUTION OF SHEAR STRESS t AND 

ISOCLiNIC PARAMETER d) ALONG THE CHORD 
ABC IN AN AXI- SYMMETRIC CASE. 
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Case (d). Axi-syoimetric Case; 

In this case, the magnitude and directions of the 
secondary principal-stresses are symmetrical about the 
middle point of the optical path in a transverse plane* 
(example; the optical path along a chord in a plane perpen- 
dicular to the loading axis of a sphere under diametral 
compression.) If the total optical path is divided into 
two equal parts, and if the first part is equivalent to a 
retarder (azimuth =<, retardation p) and a pure rotator 
(rotatory power lii), then the second half is equivalent to a 
pure rotator (rotatory power -ip) and a retarder (retarda- 
tion p and azimuth «=() . 

Now obsernre Fig. 10, The optically equivalent 
systems for both the light paths AB and CB will be same. 

Let the corresponding matrix for these optically equivalent 
systems be M. Then for the optical path BC," the matrix 
representing the optically equivalent system will be . 

So, for the total optical path ABC the matrix representing 
the optically equivalent system will be M^M. So this total 
optical path ABC is equivalent to a single retarder with 
azimuth of and ret ardaticn 2p . The quantities of and 2p 
however depend upon the rate of rotation of the secondary 
principal-stresses and the rate of retardation along the 
optical path ABC> V 




FIG -11 . DISTRIBUTION OF THE DIRECTION <)> AND THE 

DIFFERENCE (p'-q') OF THE SECONDARY PRINCIPAL- 
STRESS ALONG THE OPTICAL PATH ABC • FIGURE 
SHOWS THE CROSS SECTION OF A RECTANGULAR 
BAR UNDER COMBINED TENSION AND TORSION 
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Case (e), Anti-syinT:ietrlc Case; 

In this case the secondary principal-stress direc 
tions are anti-symmetrical and the secondary principal- 
-stress difference is symmetrical about the middle point o. 
the optical path (example; the optical path in a plane perp* 
dicular to the direction of tensile loading in a rectangulc 
model under combined tension and torsion.) 

In Fig. 11, let the optical path AB be equivalent 
to a system with a retarder (of retardation p and azimuth 
«<) and a rotator (rotatory power if)). Then the optical path 
CB is equimlent to a system with a retarder (retardation p 
and azimuth -“C) and a rotator (rotatory power The : 

optical path BC is equivalent to a rotator (rotatory power 
i{) ) and a retarder (azimuth^- and retardation p). So the 
total optical path ABC is equivalent to a retarder (retar- I 
dation p and azimuth^), a rotator (rotatory power - 2’f') and^ 
another retarder (retardation p and azimuth -<=(). This 
system is again equivalent to a single retarder and a pure i 
rotator. The primary and secondary characteristic direc- 
tions will be inclined equally, but in opposite directions ! 
to the axes of the first and second retarders or to the 
secondary principal-stress axes at the points A and C. This 
can be proved as follows; 

From Eqs. (2.15) 

tan2ti3 = tanG Sin2p( i 

= tan (180 -£) Sin (- 2c<) and Cos2e = Cos2 w Cos2o( 




FIG.12. EQUIVALENT SYSTEM FOR THE OPTICAL PATH ABC. 
OA,OB ARE THE AXES OF THE RETARDER Ri , 2i)) 

IS THE ROTATORY POWER OF THE ROTATOR. OA'.Ob' 
ARE THE AXES OF THE RETARDER R2-0P|,0'P2 
THE PRIMARY AND SECONDARY CHARACTERISTIC 
DIRECTIONS OF THE SYSTEM RESPECTIVELY. 
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= Cos2u) Cos(- 2oC) i.e. along axes inclined equally but in 
opposite directions ( =<, - to the major axis of a par- 
ticular light ellipse of elliptic ity tan©, the phase diffe- 
rences will be C and 1 80 -£ respectively. Also, the ampti- 
tudes of light vibration are equal along these axes. In 
Fig. 12. OA, OB are the axes of the retarder R-j 2^ is the 
rotatory power of the rotator, Oa', OB' are the axes of the 
retarder R2 and OP^ , o' P 2 are the primary and secondary 
characteristic directions of the system respectively* Now, 
let the incident light on this system be plane polarized 
along the primary characteristic direction OP^ . Then at the 
exit of the system, the light will be plane polarized along 
the secondary characteristic direction o'P 2 . The changes 
in the properties of the light ellipse as it passes through 
the t .0 retarders will be as follows: 

The phase difference and the- amplitude ratio after 
the first retarder along its axes will be p and tan 6^ 
respectively. The light ellipse is then rotated through 
an angle 2i> by the pure rotator after which a retardation 
of p is added along the axes of the second retarder. 

Also, the light is plane polarized at the exit of 
the system along o' P2, one of the secondary characteristic 
directions. So, the phase difference and the amplitude 
ratio of the light vibration before the second retarder 
along its axes should be ( ■r- p) and tan(- © 2 ^ respectively. 
The same light ellipse has a phase difference of p and an 
amplitude ratio of tan®^ when referred to the axes of the 
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first ret?,rder. 

So, 

tane^ = tan( - 6 2 i 
or 9^ =-02 

Also, the angle between the axis of the first 
retarder and the major axis of this light ellipse is 

^ +°f). (Because along two sets of axes equ- 
ally inclined, but in opposite directions to the major 
axis of the light ellipse, the amplitude ratios are equal 
and the phase differences differ by ir . The effective 
(taking into account the rotatory power 2^) angle between 
the axes of the two retarders is C2iIj + 2=<). 

If if) = =C-| 

tan 2o(^ = tan2e^ Cosp - using Eqs. (2.15) 

or tan2( If) 4 ^) = tan2©^ Cos p 

This relation connects the angle e^i (the angle 
between the primary characteristic directions and the axes 
of the first retarder) and , =< and p . 

Thus, in the anti -symmetric cases, the primary 
and secondary characteristic directions are inclined equa- 
lly but in opposite directions to the secondary principal- 
-stress axes at the entrance and exit respectively. (Becau- 
se, the axes of the two retarders are inclined equally but 
in the opposite directions to the secondary principal- 
-stress directions at the entrance and the exit respectively). 
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Theoretical Determination of Characteristic Parameters; 

The matrix M in Eqs. ( 3 -S'! ) in the most general 

case will have three independent quantities namely the 

three characteristic parameters. In these equations, only 

2 2 2 _2 _2 

two are independent since S-| + S2 + = S-| + S2 + = 1 , 

These axe linear equations in S-j , 82^ S^, and • 

i 

To determine these three characteristic parameters we need 
the exit light ellipse parameters S-j , Sp and S3 for atleast 
two sets of independent initial Stokes parameters S-^, S2 
and S^ so that we have four equations in the three imiknowns* 
One of these three equations can be used for checking. The 
exit light ellipse parameters for two types of incident 
light ellipse can be found theoretically by integrating 
numerically the governing differential equations when the 
functions 0 * and p’are known. From these values, the three 
characteristic parameters can be calculated. 

Experimental Determination of the Characteristic Parameters: 

The exit light ellipse parameters for two types of 
incident light can be found experimentally and the characte- 
ristic parameters can be calculated from this data* If a 
reflection polariscope is used, the primary characteristic 
direction and the characteristic retardation can be found 
directly as the combination of the forward and (reflected) 
backward optical paths is eqiii valent to a single retarder. 
However, this can be used only for comparatively thin models. 
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These quantities can be foimd directly in a transmissicn 
polariscope for axi-symmetric bodies for, in this case also 
the optical path is equivalent to a single retarder,. 

It will be very usefui,if there exist direct and 

simple experimental methods to find these parameters in the 

general case without any calculations or iterations, Robert 
14,15-16 

and Guillemet ’ ’presented a method of finding these 
parameters experimentally in a direct way. The principle 
of this method is the same as the one suggested by Kent 
and Lawson for the ellipticity measurement. In this 
method, the ellipticity of light from the analyser is 
calibrated in terms of the ratio of the constant and alter- 
nating components of light intensity, when the analyser is 
rotated at a particular angular frequency a • In this case, 
the intensity of light from the analyser is given by I = 

(1 + Cos2ai Cos2a)t) where K 2 is a constant, Eqs. (2.15)» 
This method however requires additional instrumentation 
for measuring the alternating and constant components of 
voltage. To find the primary characteristic direction 
using this method, the polarizer at the entrance of the 
model is rotated until the ellipticity at the exit is gero. 

We shall now describe two direct and simple meth- 
ods to determine these characteristic parameters experi- 
mentally, which dbnot involve the rotation of the analyser 
at a particular speed or any additional instrumentation. 
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Method 1 : 


If the light at entrance to the model is circular- 
ly polarized, the light ellipse at exit will have its major 
axis at an angle of 4-5° to the secondary characteristic 
directions. This is so since the amplitudes along the pri- 
mary and secondary characteristic directions will be equal 
in this case (if the absorption of light along the optical 
path is neglected). If the light path is considered to be 
equivalent to a pure rotator and a retarder, the pure 
rotator will not have any effect on the incident circularly 
polarized light while the retarder increases the phase 
difference along its axes from 90® to (90 +S) where 5 is 
the characteristic retardation. Along these axes the amp- 
litudes will be equal even after the addition of retardation. 
These axes correspond to the secondary characteristic dire- 
ctions. This angle of ’(90 +5) can be measured by using 
Tardy's method. One of the primary characteristic direc- 
tions is found by rotating the polarizer until the inten- 
sity along one of the secondary characteristic directions 
is minimum. The final position of the polarizer axes are 
the primary characteristic directions. 

Experiment al Pr oc edure ; 

Keep the first quarter -wave plate axes at 45® to 
the polarizer axis so that we have “a circularly polarized 
light at entrance to the model. With the second quarter- 
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-wave plate out of the field, rotate the analyser for mini- 
mum intensity and note this direction . This direction 
corresponds to the minor axis of the light ellipse. Along 
the directions the amplitudes of the components 

of the light vector are equal and these directions corres- 
pond to the secondary characteristic directions. 

Now, introduce the second quarter -wave plate so 
that its axes are parallel to the and D-j + 90^ directions . 
Rotate the analyser for minimum intensity which will he 
absolute minimum. Let the minimum intensity position be 
The 2 (D 2 - ) = (90 +5). Now remove both the quar- 

ter-wave plates and keep the analyser axes along one of the 
secondary characteristic directions. Rotate the polarizer 
for absolute minimum intensity. The final orientation of 
the polarizer axis corresponds to one of the primary chara- 
cteristic directions. 

Method 2 : 

The principle of this method is the same as that 
of the intensity calibration method described in chapter 2. 

Procedure; 


Let the optical set up consist of the polarizer 
(at arbitrary orientation), the stressed model, the quarter- 
-wave plate (at arbitrary orientation) and the analyser. 
Rotate the analyser for minim-um intensity thus finding the 
orientation of the major and minor axes of the light ellipse. 
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Rotate the analyser throiigh 45° in either direction and note 
the intensity with a photometer. All ellipses ( of 
different orientations and ellipticities ) will have an in- 
tensity equal to along directions which make an angle 
of 45 ° with their major or minor axis. Now keep the angle 
between the axes of the second quarter-wave plate and analy- 
ser equal to 45*^ and rotate the polarizer until the intensi- 
ty from the analyser is equal to the previously found value 
of The final orientation of the axis of the polarizer 

corresponds to one of the primary characteristic directions. 
The reasoning for this is as follows; When the intensity 
of light from the analyser is equal to it means that 

the axes of the light ellipse are at 45° to the analyser 
axes. That is, they are along the axes of the second quar- 
ter-wave plate. Since the phase difference alaig the major 
and minor axes is 90 °} the phase difference along the same 
axes but before the second quarter-wave plate must have 
been zero 'v^ich corresponds to the case of plane polarized 
light. The orientations of the polarizer axes which pro- 
duces plane polarized light after the model correspond to 
those of the primary characteristic directions. 

Now remove the second quarter -wave plate and rotate 
the analyser for absolute minimum intensity. The final 
orientations of the analyser axes correspond to the secon- 
dary characteristic directions. The characteristic retar- 
dation can be found as in the first method or by any other 


method 



73 


VJith minor modifications, the first method can he 
used in the scattered light technique where the model acts' 
as an analyser and the second method in the technique where 
the model acts as a polarizer. These will be described in 
the next chapter. 

It may be noted now that in a transmission polari- 
se ope only the fractional part of the characteristic retarda- 
tion can be found directly while in general, it is difficult 
to find the integral part of the characteristic retardation. 
The usual method of using the Babinet-Soleil compensator with 
white light, as in two-dimensional problems can not be used 
here since the characteristic retardaticu in general, is a 
function of not only the wavelength of light, but also the 
orientation and magnitude of the secondary principal-stresses 
along the optical path. The importance of the concept of 
optically equivalent system and the characteristic parameters 
lies in the fact that in a three-dimensional model, in gene- 
ral, the preferential directions (directions which can be 
directly found by experiment) are the characteristic direc- 
tions and not the secondary principal-stress directicais* 

In this chapter, the proof for the existence of 
an optically equivalent model is presented in detail. The 
behaviour of light when the light path is reversed is dis- 
cussed. The relationship between the two optically eq'uiva- 
lent systems one ■'•dth the retarder followed by a pure rota- 
tor and the other with a pure rotator followed by a retar- 
der, is derived. Using this information, the optical 



74 


phenomenon for each of the following three cases are discussed. 

(a) Reflection Polariscope 

(b ) Model having axi-symmetric stress distribution 

(c) Model having anti -symmetric stress distribution 

Two direct and simple methods to determine experi- 
mentally the characteristic parameters are presented. These 
methods can be extended for use in a scattered light pola- 
riscope also. This will be discussed in detail in 
Chapter 



CHAHTER If 


sc ATT SUED LIGHT METHODS IN THREE-DIMENSIONAL 

PHOTOELASTICITY 

Unlike the stress -freezing 'technique, the scattered 
light method is non -destructive for analysis of stresses in 
three-dimensional photoelasticity. In the scattered light 
technique, the models can be subjected to live loads at 
roan temperatures and analysed, and no frozen stress model 
is necessary. In this chapter, we shall first make a brief 
survey of the existing techniques in scattered light photo- 
elasticity, investigate their limitations, and then suggest 
refined technioues. 

17 

It was first suggested by Weller that the three - 
-dimensional photoelastic analysis could be done by measur- 
ing the intensities of the scattered light components. His 
suggestion was based on the theory of absorption and re- 
-emission of light by vibrating particles. • . 

According to this theory, (i) the light that is 
scattered at right angles to the direction of propagation is 
always linearly polarized (ii) the intensity of this 
scattered light is proportional to the square of its appare- 
nt amplitude when viewed from that direction. By ’apparent 
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amplitude’ is meant,. the component of amplitude normal to 
the line of sight. His method in the three-dimensional case 
consists of computing the secondary principal-stress differe- 
nce from the minimum spacing of the scattered light fringes. 
That is, the model is rotated about the light path under 
consideration until the spacing of the fringes becomes mini- 
mum. If ’d' is the fringe spacing, the secondary principal- 
-stress difference (p' - q' ) is given by 

I I ^ 

(p - q ) = — where F is the material fringe constant 

( P si .in/fringe ) (^.1 ) 

The directions of the secondary principal-stresses 
are found by rotating the model until the fringe system 
vanishes. This method doesnot take into account the effect 

of rotation of the secondary principal-stress axes. 

18 

Drucker and Mindlin investigated the effect of 
rotation of the secondary principal-stress axes along the 
wavenormal on the retardation in a photoelastic model. They 
obtained an analytical solution for the case where the prin- 
cipal-stress magnitude remains constant and their directions 
vary linearly along the wavenormal. They concluded that 
(i) the rotation increased the relative retardation by a 
factor S = (1 + where R is the ratio of the rate of 

rotation to the rate of retardation and that (ii) at high 
stress, the light vector rotated with the secondary primcipal- 
-stress directions; plane polarized light entering along a 
direction of one of the secondary principal-stresses rema- 
ined practically plane polarized in the direction of the 
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stress through the entire model. They suggested that the 
factor 'S’ be incorporated in the scattered light analysis 
so that we have 

(p' - q' ) = . (^. 2 ) 

Sd 

37 

Saleme suggested a method of analysis which in- 
volves the measurement of absolute intensities in two sets 
of mutually perpendicular directions, the angle between 
which is His analysis is as follows; 

Let = (S'], S 2 , S 3 ) be the Stokes vector of light referred 
to the secondary principal-stress axes at a particular point 
in the model. 

s-i cosse 

= S^ = Sin2eCose (^.3)(a) 

S^ Sin2e Sine When the same light ellipse is 

referred to OXY making an angle <=( with Opq, we get 

S-| Cos20 Cos2c< + Sin20 CosC Sin2=C Cos2e^ 

¥ 2=^2 = Sin2e Cose Cos2=c - Cos2e Sin2=( = Sin2&i Cose-i 

S^ Sin20 Sine Sin20^ SinC-i 


(4.3)(b) 


Let the amplitudes of light vector components along X,Y-axes 


be a^jr and a^.* tanS-j = is phase . difference corres- 


2 2 . — 2 
ponding to Stokes vector Vx* a^ - a-^ = S^ = 2 a^ - 1 


1 - 2 a: 


Y 


■^,2 . .2 _ 1 


The intensities and along these axes are given by 

■ ? 1+ S'], 

Ix = KaJ = K( — ) 


2 


Q+»h-) 
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1 4- S 

= Ka| = K ( ■^) (4-.5) 

where K is a constant. 

Now consider the intensities and 1^' along the 
axes which are inclined at an angle of ^5° to the X,Y-axes, 
Fig. 13 . These are obtained hy substituting (=< +45°) for 
in the expression for in the above equations. 

So, = ■^ ('I - C os2 0 Sin2‘=< + Sin2© Cos0 Cos2c<) (4.6) 

V = I Cos2e Sin2o( - Sin2e Cose Cos2c<:) (4,7) 

Let z denote the distance along the wavenormal. 
Neglecting the rotation of the secondary principal-stress 
axes, i.e. 0-^ and ^ are independent of z, we have 

(4.8) 

(4.9) 

(4.10) 

(4.11 ) 

(4.12) 

from Eqs. (4.4) and (4,5) 

dix dl^ 

‘ 31 “' 

dz dz 


dl. 


X = 


dz 


- ^ Sin2e Sin2«( Sin^ ~ 
? dz 


~dz 


Sin2Q Sin2=< Sin^ ^ 
^ dz 


£1^'= - -I Sin2e Cos2=( Sin£ 
dz 2 dz 


^ ^ Sin20 Cos2p( SinC ~ 
^ dz 


Hz 


+ ly 


K 


tan2x 


(4.13) 
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from Eqs. (4.8) - (4.11 ) 


= KS^ = KCos2e^ • (4.14) 

iv' - ly'- Sinae^ cosCi (4.15) 


From the above Sqs. (4.12) - (4.15), the quantities 
K, ==< , 0-| and £'-| can be computed from the values of the 


intensities I,r Iv' and Iv' • 

A, X, X X 


Since (p 


' de 

q ; = F 

dz 


(4.16) 


The quantity (p - q ) can be calculated if the rotational 

effect is neglected. 
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Jessop attempted to arrive at a qualitative 
estimate of the magnitude of the rotational effect and 
deduced certain conditions under which this effect will be 
small enough to be neglected. He considered the Max^Arell- 
-Neumann equations 


M =z M 
dz dz 


+ 


(■ 


fX - fx 
ax ax 


) Sine 



(4,17) 


dp 


= a-^ Cos e ^ 

i do 


and 


day 

dp“ 


ay Cos ^ 

i dP 


(4.18) 


Mow, a relation connecting the quantity (to be 
determined) and the quantity £ (that can be observed experi- 
mentally) is to be arrived. 

The quantity is equal to at points where 
dz dz 

au, = a„ or SinC = 0. The first condition can be 


either 
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satisfied at the point of entry of the beam into the model, 

by orienting the model so that the secondary principal- 
-stresses at this point are inclined at to the direction 
of vibration of the incident light, but it is not in general 
satisfied elsewhere. The second condition is satisfied at 
points where whole or half-fringes appear but the value of 

at these points cannot be obtained from a smooth curve 

G 2 

plotted through the observed values. Also, in the absence 

of rotation, the quantity will be equal to To find 

dz dz 

the effect of neglecting the second term in Eq. (U-,17) which 
represents the rotational effect, we consider, 


dp 


^ay ax'' 


(2 + ^ + Cos€-|^ 


a. 


from Eq. (4.18) 


(4.19) 


^ 4Cos cf, if ay and ay are initially equal. 

If both 4^ and its derivatives are small 
dP 


4-C,osC 


Sine 

dedp 


I- 

4 — (Sine M) 
tip de 


so that cl|) _ fl, f:; 

y ay 

a^ a^ 


SinC 


Sin e-^ ~ 4%- 4? Sin^e 
dp d o dP 


where y 


d0 . d0 
d 0 dp 


frp, (4.19) 
^ . 2 

= y Sin £ 
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So, u is a small variable coefficient of the order 

estimate of the magnitude of ^ in any particular 

dp 

case can be made by observing the dimunition of the sharpness 
of the fringes In the model after rotating about the light 
path. If the total change in the orientation of the secon- 
dary principal-stress axes over the light path in the model 
is 30° and the relative retardation introduced (due to stre- 
sses only) is n wavelengths, the value of can not exceed 

dp 

(1/12n). If the fringes are sharplyly defined over a range 

of two or more wavelengths, the error introduced in by 

dz 

drawing the , z graph as a smooth curve through the points 

obtained by measurements of the positions of the fringes 

will not exceed something of the order of and this will 

generally be less than the possible error of measurement. 

In 11 IP 

Based on the above findings, Frocht and Srinath 
developed methods for measuring the secondary principal- 
-stress difference and their orientation. .An account of 
these methods will now be presented. 

Consider Eqs . (2.1) - (2.5) of chapter 2. Let 
X-and Y-axes represent the directions of the secondary 
principal-stresses. The intensity for a direction of 
observation making an angle with the X-axis is given by 
the expression 

= 3.^ Sin‘^=< + a.j Cos^^ - 2a.-yr_ ay Coso( Sino( CosG (2*5) 
if the proportionality constant is assumed to be unity. 

This light intensity will be maximum when the 
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direction of observation is along the minor axis of the light 
ellipse and is a minimum ^7hen the direction of observation 
is along the major axis. That is, the intensity is a maxi- 
mnm or minimum when the angle satisfies Sq. (2.3) 

tan2 (>< ■= -^ — “Cose = tan2 eCos e (2.3) 

From Eq. ( 2 . 5)5 note that for any fixed direc- 
tion of observation o(and for fixed values of and aYj the 
intensity of scattered light varies from a maximum to a 
minimum depending upon the value of 0 . Idien oC is less than 
'n /25 the intensity is a minimum (although not necessarily 
zero) when £= 2mv and becomes a maximum vrhen 0 = (2m - 
, When is greater than v/2, the condition for maximum 
and minimivn will be reversed. Except for special cases, £ 
will vary from point to point along the axis of propagation 
and in general, will assume the values necessary for a maxi- 
mum and a minimum light intensity. This gives rise to a 
sequence of bright and dark regions which forms the scattered 
light pattern. Also for approximately equal amplitudes 
a^ and a^, the axes of the light ellipse are at 45° to the 
directions of the secondary principal-stresses and therefore 
this direction gives a pattern of sharpest contrast. 

If the directions of the secondary principal-stress 
axes donot vary along the optical path and the amplitudes 

of the light components along these axes at the entrance to 

^ ' 

the model are equal, (i) these .remain equal all along the 
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optical path and. (ii) whatever retardation is added along 
these axes at the entrance, the same will appear as an addi- 
tional retardation at all points on the optical path in the 
model. We then have, 

(p' - q' ) = P ^ (4.20) 

dz 

A compensator can be placed at the entrance to the 
model with its axes parallel to the secondary principal-stress 
axes. If the axes of the polarizer, which is kept before the 
compensator make an angle of 45° ’-/ith the axes of the compen- 
sator, the amplitudes along the secondary principal-stress 
axes in the model will be equal. The compensator then acts 
as a continuation of the model. Along directions making an 
angle of 45° with the secondary principal-stress axes, the 
scattered light fringes will have the sharpest contrast. 

This is so since, the major and minor axes of the light 
ellipse at any point on the optical path in the model will 
be along these axes. From these fringes, the n-z curve can 

be plotted from which the slope — can be computed. 

dz 

Even in the presence of rotation of the secondary 
principal-stress axes, the above method of plotting the 
n-z curve and thus calculating the -secondary principal -stress 
can be used provided that the conditions mentioned by 
Jessop ^^are satisfied. 

Directions of the Secondary Principal -Stresses ; 

' ■ ■ , . ■ ^ ' 

Approximate Methods; 

The approximate locations of the secondary 
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principal -stress axes can be readily found by visual observa- 
tion vitliout a compensator or a photometer. 

(a) Displacement Pattern Method: 

When looking along either one of the two secondary 
principal-stress axes, the center of the local scattered 
light pattern is of uniform brightness, and the fringes on 
opposite sides of this region become displaced relative to 
each other by half a fringe. The discontinuity in the frin- 
ges arises because, for a direction of observation on one 
side of the uniform region =< is less than y/2 and on the 
other side it is greater than ir /2. Therefore, the same 
retardation 2mir will produce dark fringes on one side and 
bright fringes on the other as can be seen from Eq. (2.5)» 

The stress pattern resembles the picture of an X-ray of the 
spiral column with emanating ribs. 

(b) Sharpest Ccxitrast Method; 

The model is rotated until the scattered light 
stress pattern appears of sharpest contrast to the naked 
eye when viewed along the line of vision of the photometer. 

In this position, the secondary principal axes are approxi- 
mately at + ^5° to the line of vision of the photometer. 

(c) Minimum Intensity Method; 

If the light is circularly , polarized at the entrance, 
the axes of the light ellipse at any point are independent 
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of the birefringence and are always inclined at 4-5°, approxi- 
mately to the directions of the secondary principal-stresses. 
Hence, the determination of cne of the axes of the ellipse 
also determines the directions of the secondary principal- 
-stresses. The direction of the major axis of the ellipse 
can be readily and accurately located with a photometer 
since, the intensity of the light is a minimum for observa- 
tions along the major axis. 

(d) Exact Methods; 

(i) The Constant Intensity Methods 

l-ftaen ^ is small, the changes in the amplitudes 
dz 

a,^ and a^, when retardation in the compensator (whose axes 
are set parallel to the directions of the secondary princi- 
pal-stresses) is varied, are negligible. So, from Eq. (2.5)?- 
when = 0, the intensity is independent of the bire- 
fringence and it remains constant even thoi:igh the retarda- 
tion at that point is changed. 


( ii) The Zero Intensity or Line Ellipse Method; 


From Eq. (2.5 )j it can be observed that for a fixed 

value of°{ ( <ir /2), the intensity is a maximum when 6= (2m - 1 ) 

and a minimum when £ = 2mTr , In both these cases, the 

light ellipse degenerates into a straight line. I'flien 

„1 ay. 

6 = 2mrr , the line is inclined at an angle = tan 

and when' e = (2m - 1 the line is inclined at an angle 
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of -1 T^Y\ 

= tan X-axis. The intensities of scattered-' 

light for obse^ations along these lines will be zero* To 

determine the directions of the secondary principal-stresses, 

the retardation in the compensator is adjusted to make the 

intensity at a point minimum for an arbitrary direction of 

observation say OM, Fig. Ilf, Assuming =< to be less than 

'^/2, ^ will then have assumed a value of 2mir and the light 

ellipse will have become a straight line CO, Fig. 14. making 

an angle tan"*"^^--^) with the X-axis, which can be located 

ay 

with a photometer. Next, keeping approximately the same 
directicn of observation OM, the retardation in the compen- 
sator is adjusted to give the maximum light intensity. When 
this occurs, the light ellipse will have become once again 
a straight line, only this time it will be DO inclined at 

an angle = tan" C-^) to the X-axis. This line is also 
I ay 

located with a photometer. The directions CO and DO, are 
equally inclined with respect to the direction OX which is 
one of the secondary principal-stress axes. 

(iii) Absolute Intensity Method; 

This method which involves the measurement of the 

39 

absolute intensity has been suggested by Menges and 
Salame .. This has been described earlier. The main dis- 
advantage with this method is that it involves the measure- 
ment of absolute intensity which will depend upon the absorp- 
tion along the varying optical paths. 




FIG. 15. SCATTERED LIGHT TECHNIQUE USING THE 
MODEL AS A POLARIZER. 
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40 

Cheng described a new method which records simul- 
taneously the intensities of scattered light along two direc- 
tions of observation 45° apart in a plane normal to the beam 
of light. The theory of this method is the same as the one 

37 

described by Salame except that the effect of varying 
lengths of light path is eliminated by recording the inten- 
sities before the loading of the model. He made use of 
Druclcer and Mindlin's conclusion that for cases where the 
parameter R is small, plane polarised light entering the 
model along one of the secondary principal-stress axes at 
the entrance follows that axis throughout the model. The 
polarizer axis is set along this direction and both the pola- 
rizer and the model are rotated in unision until the light 
intensity is minimum at the particular point under conside- 
ration. l^fhen the minimum occurs, the direction of observa- 
tion is along one of the princijal-stress axes. 

Robert and Guillemet ^“^’^^a^pted an altogether 
different a]?proach towards scattered light analysis in the 
three-dimensional cases. Since any length of light path 
in a photoelastic model under load can be replaced by the 
combination of a linear retarder and a pure rotator, the 
characteristic parameters of this equivalent system can be 
found at each point on the optical path with respect to the 
point at the entrance to the model, Fig. 15* After knowing 
these parameters, the secondary principal-stress difference 
and orientations can be found by using the governing optical 
equations. To find these parameters, they adopted the method 



that Is popular in crystal physics, of rotating the analyser 
at a particular speed and then measuring the ratio of alter- 
nating and constant components of light from the analyser. 
This ratio is dependent only upon the ellipticity of light 
from the analyser (see Chapt-er 3)* The light incident on the 
model in this method is unpolarized and the model acts as a 
polarizer. The intensity from the analyser is given by 

1 = K 

2 

where is the first Stokes parameter of the light 
ellipse and 1 is a constant, 
or I = I (1 - Cos2o Cos2'=C) 

Tan 0 ) is the ellipticity and is the azimuth (with 
respect to the analyser axes) of the light ellipse incident 
on the analyser. If the analyser is rotated at an angular 
velocity of w , '=< = co t (t is the time.) The ratio of the 
alternating to the constant components of this intensity is 
obviously 

1 - e^ 

C os2 o = s 

1 + e*^ where e = tan o) 

Also Cos2w = _ Sin^2 ® Sin^C 

= Y ( say) 



then Cos 6 . 
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To find the primary characteristic directions the 
model is rotated with respect to the incident light -vibra- 
tion -until the value of Y observed is zero, which corres- 
ponds to plane polarized light (e = 0) before the analyser* 
The position of the incident light beam with respect to the 
model axes will be the primary characteristic direction. 

Now this beam is rotated through an angle of 1+^° and the 
val-ue of Y is noted. Then Y = Cos^ where G corresponds to 
the characteristic retardation. The direction of the major 
axis of the light ellipse before the analyser when the inci- 
dent light beam is along one of the primary characteristic 
directions corresponds to one of the secondary characteris- 
tic directions. For finding these parameters the model is 
rotated about the axis OO' , Fig. 15* The§e parameters are 

fo-und for the optical laths 0lyL| , OM 2 , .... 00*. 

13 

Aderholdt, Mckinney, Ranscxi and Swinscn in their 
scattered light method assumed that for problems where the 
rotation to retardation ratio is small, the light waves 
follow the secondary principal-stress axes. That is, the 
amplitudes of the canponents of light along these axes re- 
main constant along the light path. This assumption is 

1 o 

again based on the findings of Drucker and Mindlin . Keep- 
ing the Ught at the entrance to the model circularly pola- 
rized, the. directions at each point of the optical path in 
the model along which the fringes appear to be fading (or 
the intensity is constant) are noted which give the directions 
of secondary principal-stresses. The fringes are counted by 
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moving along the light path and at the same time, rotating 
the observation angle such that the angle of observation at 
any point is ^5° from one of the secondary principal-stress 


axes . 
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Srinath made an analysis of some of the previous 
scattered light techniques. He suggested a method called 
mini -max method, which is a combination of the zero and 
absolute intensity methods, to find the directions of the 
secondary principal-stresses. In this method, the orienta- 
tion of these directions is found from the formula 

’^1 +"^2 


tan c< 


V- 


‘3 + 


where 


2 2 

= K (a Sln°C + b Cos=<) = (!=<) 


= K (a Sin°c - b Cos°c)^. = (Ic() 
^ . ) • 

= K (a Coso( + b Sino()^ = (Ip) 


= K (a Cos°( - b Sin=<)^ = (Ip) 




and the retardatioi p is obtained from the equation 


max 


mtn 


max 


min 


Cos P = 


- 2 L + (y-i +^ 2 ) 


(yi - yg) 


where 


Jc( = K (a^ Sin^oC + b^ Cos^oC - 2ab Sinp( Coso( Cos p) 
- obtained from Eq. ( 2 . 5 ) 



95 


Critical Analysis ; 


A critical analysis of the above described scattered 
light methods will now be made. Except the method described 
by Robert and Guillemet, all the methods are. ccnnected with 
the findings of Drucker and Mindlin . These findings in 
turn are based on the analytical solution for the case where 
the secondary principal-stress difference is constant and 
their orientation varies linearly along the optical path. 
These results are confined to the cases where the parameter 
R is small. 

We shall now examine the validity of the assumptions 
involved in the previous scattered light methods for the case 
where the rate of rotation 0 ' can be any arbitrary function, 

but the ratio R = is a constant. 

p' 

The analytical solution for this case is given by Eqs. (2.^6) 

- (2.4-8). These are; 

•^1 1 ■ 


^1 


2RK‘ 


^ ^ (1 - Sin(2KiZ) -f A 2 ) 


^2 

$2 = (1 - Cos(2Ki3 + Ag) 


and 


_ Si 
S3 = ^ 


k2 


(-1 

2RK 


Sin (2K0 + A 2 ) 


( 4.21 ) 

(4-. 22) 
(4.23) 


where A-| A 2 are arbitrary constants which are to 
be determined from the initial conditions. 

Let us now examine this solution for two types of 
incident light, (i) plane polarized light alcog cne of the 
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secondary princixal-stress axes at the entrance and (ii) 
circularly polarized light. 


C ase (i) ; Incident Light is Plane Polarized Along one of 
the Secondary Principal-Stress axes At the Entrance 


At Z = 0 , 0 = 0 (say), S-^ =1,82=83^0 


From Eq. ( 4 , 23 ) 

A-i = - 1/2R 

and from Eq, ( 4 . 21 ) 

1 = -J + 1 (1 1 )i 

4e2 + 1 K ItRS + 1 


SinA2 


i.e SinA2 = 1 
or A^ = ir /2 , 

So, at any point on the optical path we have 

= — (COS2K0 - 1) 

2KK2 

The light is plane polarized at points where 

83 = 0 

i.e. when 2K0 = 2mr 
At these points 

g _ - _ ^ i ('I — ^_)‘2' 

4R^k2 K 4r2i? 

= 1 


f 


( 4 . 24 ) 
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l.e. the intensity attains its absolute m.aximum ’value. So, 
the light vector will be along the direction of one of the 


s ec on dary pr inc ipal -stresses. 


The change in the orientation of the secondary 


principal-stress axes between two fringes is given by 
3 q» (^.24-) as 

d0 = r/K = ddp 


C-^.25) 


i.e. dp = 

i.e. betv/een two fringes, the change in the retardation dp 
introduced due to the stresses only is given by Eq. (4-. 25). 

If the effect of the rotation of the secondary principal- 
-stress axes is neglected, this is given by 

dp = 2. TT 

So, the effect of rotation is to decrease the fringe 
1 1 

spacing by a factor — or 1 , 

2KR (4j^ 2 ^ ^ )-g- 


Case (ii) ; Incident Light is Circularly Polarized. 


At z = oj = S2 = o; S3 ^ I , = 0 

Ercm Eq. (4-, 23 ) 

A^ = 1 

and from Eq. (4-. 21) 

° ' S? 

i.e. SinA 2 = 1 or A 2 = tr /2 
At any point on the optical path, 

Cos2Kj^ 

2 " 4-r2 4 . 1 ^ 4-r2 +1 
== f 
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if Cos2K0 = - 

i.e. 2K0 = (2n + 1 ) -^ 
At these points, 


S, 


2RK^ 


+ ~ (K^-l)V^Cos2K0 


1 1 

■ o o 9 os2K0 

2RK"' 2HK^ 


1 


(1 - Cos2K0) 


2Rr' 

2R (1 + 4r2) 

(4r2 + 1 ) 




0 , 


So= --(1 


- Sin2IC0 


(4.26) 

(4.27) 


1 

(4r 2 + 1 )i 


Sin2K0 


- Sin ^ — 

2 


(4.28) 


The second Stokes parameter S^ gives us the inten- 
sity along axes making an angle of 45° with the reference 
axes (in this case, the secondary principal-stress axes). 

From Eqs. (4,2?) and (4.28) we observe that between 
two successive fringes. 


d0 = Rdp = ^ 


(4.29) 


So, in this case also, the effect of the rotation 
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of the secondary principal-stress axes is to decrease the 

fringe shacing by a factor — — i . The same conclusions 

(4r2 + 1 

1 © 

were dram by Drucker and Mindlin for the case where both 

of the fijnctions 0* and p'are constants. 

Now, we shall examine the above analytical solution 

in light of the assumptions of the previous scattered light 

methods. In the scattered light method developed by Frocht 

10 , 11,12 

and Srinath , the following assumptions are made: With 

the addition.' of a certain amount of retardation AGq at the 
entrance along the secondary principal-stress axes when the 
amplitudes of light vector components along these directions 
are equal, (1) a retardation of appears as an additional 
retardation at every point on the optical path along the 
secondary principal-stress axes and (2) the amplitudes (or 
the intensities) along the secondary principal-stress axes 
at any point on the optical path do not change with the 
change in the retardation added at the entrance. At 
Z = 0 let =0 (i.e. the amplitudes are equal along the 
seccodary principal-stress axes at the entrance). 

= 310 ^ 0 , 0 S2 = Cos€q 


The light intensity I = K a- 


X 

_ 1 + Si 
K(_-^ 


) 


where K is a constant . 


d^Q 2 d0Q 


(^. 30 ) 
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Let us eTaluate the constants and A 2 of the 
solution for these initial conditions. 

From Eq. (^.23) 

= SinC^ = A-^ (4 - .31 ) 

From Eq. (4-. 22) 


Sin^" 


0. 


S 2 = CoseQ = (1 - Cos(2K0j_ + A 2 ) 


or Cos ( 2 KjZij _ + A ^) = Cos£q 


(4.32) 


(1 - 


Sin^^ 


1/2 


and from Eq. (4.21 ) 


Sin(2K0j_ + Ag) 


Sine 


n. 


Sin2£o^l 

2HK(1 )i 

k2 


(4.33) 


Substituting these expressions in Eqs. (4,21) - (4.23) so 
that the constants A^ and A 2 are eliminated, we obtain the 
following expressions for S-j , S 2 and at any point. 


Sine 


Si = - 


0 1 


2Rr 


S 2 = C 0 SK 2 coseQ - 


~ + (SinK2 CoseQ + 
SinK2 SinGo 


CosK„ Sine 


0 


2RK 


) 


2RK 


and 


SinC, 


Q3 ^ ^2 ^ ^ (SinK2 CoseQ + 

where K 2 = 2K(i9 - 0j_) 


C0SK2 SincQ 
^RK “ 


(4.34) 

( 4 . 35 ) 

(4.36) 

(4 .37) 
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dl _ ^ 

^0 (1 + L 

- CosK^ Cose^^ 

= 0 

for all ■values of 0 


[coseQ + SinK2 Sine^ (1 + 4-R^)i 


0 


when 


= 2nir = 2K(0 


01 ) 


(4.38) 


Also 


dl 


dc 


0 


0 


for all values of 0 


0 


when R«1 (4.39) 

The phase difference 0 along the secondary principal-stress 
axes is given by 

Sin0Q SinK 2 CosGq C 0 SK 2 Sir£o 




2RK 


4r2k2 


tane 


h_ 

^2 SinKp Sin0o 

CosK^ CosCn - 

^ 0 2RK 

4r^ Sin0o + (4R^ + 1)^ SinKp CosCb + CosKp Sin^ 
(1 + 4R^) C 0 SK 2 CosCq - (1 4- 4R^)'5' SinK2 Sin€^ 


tanS, 


0 


if K 2 = 2nv 


(4 .40 ) 


tan€ tan(0Q + K^) 


(4.41) 

if 4R^«1 

o 

The quantity hR may be neglected when compared to 
unity if it is less than 0^01 i.e. for R<0.05» 

So, the assumptions in the above scattered ligljt 
method may be considered valid for the cases where the 
parameter R is less than 0.05« 
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13 

Aderholdt, Mckinney, Ranscn and Swinson , in their 
scattered light method asstune that the intensity along the 
secondary principal-stress axes remains almost constant (for 
certain class of problems where the rate of rotation is not 
high) along the optical path. 

Also from the above solution, i.e. Eq. (4.3^), 

b 

•0 


M = - ^ = 0’ 


dz 


2RK 


j - CosK^ CosCq) 


< 20 ’ 

i.e, da^ <2d0 

2 2 2 2 
(where ^ ) 

i.e. the intensity variation along the secondary principal- 
-stress axes depends upon the rate of rotaticai and not on 
the larameter R. The change in the square of the amplitude 
(and so the intensity) will be less than 0.1 if the rate of 
rotation is less than 0.05 radians per one tenth of a milli- 
meter. The above scattered light method may be taken to be 
valid for rates of rotation lesser than this value. Even if 
these assumptions are satisfied at each point, the errors 
involved may be cianiiLative and prove to be disastrous over 
long light paths. This sets. a serious restriction on these 
methods. Although the above remarks are based on the solution 
for a special case where the ratio R is a constant j- this 
solution holds over a length of light path evai in a very 
general case. This makes the above remarks relevant. 
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38 

Jessop from his analysis concludes that the error 
introduced by drawing theg , z grajti as a smooth curve through 
the points obtained by measixreaent s of the positions of the 
fringes will be p Sin^0 where y is a small variable coeffi- 
cient of the order (4^). In his derivation he assumes that 

dP 

both the quantities and its derivatives are small so 

dp 


that 


Cos 


dp 


A- (Sine-^) 

dp d£ 


Sine 


d^0 

d£dp 


^ (Sine 

dp dp 


Analysis as Applied to a Model Under Combined 
Tension and Torsion : 

Figs. 16-25 show the plots of various quantities 

for the case of a model under combined tension and torsion. 

The model material is paraplex. The material fringe constant 

is 139 Psi. in. /fringe. The values of shear stress are com- 

42 

puted from the series solution, Timoshenko . The integra- 
ted optical effects are computed by numerically integrating 
the governing optical equations in terms of Stokes parameters. 
Runge-Kutta-Gill method of numerical integration is used. 

It was foimd that the results didnot vary -f/ihen the nuaaber 
of steps is increased from 50 to lOQ. At each step the 

values of the three Stokes parameters S-j , S2 and S3 were 

2 ' 2 "' '' 2 '.' ' 

found to satisfy the equation + S 2 + S 3 = 1 , 
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Frcan Figs. 16 and 1 7 ve obserre that the quantities 

R and^-S are of the same order of magnitude along most of the 

length of the light path in the model. Hence Jessop’s 

assTomption of neglecting the term Sin£ compared to the 

ot 

term R.'jCX) 5€ is not justified in these cases. It can be said 
■that in general, the quantity need not be small and negli- 
gible. 


The conclusions of Drucker and "lindlin^^ are based 
on the solution for the case where 0' and p' are constants 
and their ratio is a small quantity. However, these con- 
clusions cannot be extended to the general case as can be 
seen from Figs. 19 -21, Fig. 19» shows that the phase di- 
fference alcng the secondary principal-stress axes calculat- 
ing without considering the rotation is greater than that cal- 
culated with rotational ef feet along some part of the op)ti- 
cal path. This is contrary to the prediction by Drucker 
and Mindliri^^ 

Now consider the assumptions in the methods descri- 
bed by Frocht and Srinath^^’^^’.^^These are not satisfied for 
the particular loadings considered as can be seen from Figs. 
22 and 23 . The variation of the amplitudes alcxig the seccn- 
dary principal-stress axes with -distance z along the light 

path is quite marked for these cases, Fig. 18. and the method 

13 ' 

described by Swins on et . al is not applicable here. So, 
before using either of these methods, one has to check for 
the problem under consideration whether the rate of rota- 
tion 0* or the parameter R is within the limits suggested in 



24. VARIATION OF THE SEC. PR. SIR. AXES ORIENTATDN 
4 ALONG THE OPTICAL PATH IN THE MODEL. 
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the previous section. The plot of JZ5 Vs. z and (p*- q*) 

Vs. z are presented in Figs. and 25* for reference. 

14- 15 16 

The method described by Robert and Guillemet- ’ 
doesnot involve any of these above mentiaied assumptions. 

But the analysis here is through the characteristic para- 
meters. This method requires the rotation of the analyser 
at a particular speed and some additional instrimientation . 


New Methods: 

It is highly desirable and useful to have methods 
which donot have any of these drawbacks. These methods should 
essentially be based cn the Maxwell -Neumann equations or 
their equivalent forms since these are the equations availa- 
ble at present which describe the optical iSienomena in the 
presence of the rotation of the secondary principal-stress 
axes, accurately as sho^-m by various dynamical derivations 
(Chapter 2). It is now clear that at any point on the optical 
path in a photoelastic model under three-dimensional state 
of stress in general, the preferential directions (the dire- 
ctions which can be determined directly by experiment) are 
the characteristic directions and not the secondary princi- 
pal-stress directions. So, in general, the secondary princi- 
pal-stress axes cannot be found directly from the experimen- 
tal at a single point unless it is coupled with the data 

at the neighbouring points. 

The scattered light technique involves the follow- 
ing steps; (1) measurement of scattered light intensity 
along different directions taking into account, the variation 
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of the lengths of the light path involved, (ii) finding the 
parameters of the light ellipse and the characteristic para- 
meters at the point under consideration and (iii) evaluating 
the functions 0 and p using these values and the governing 
optical equations. 

One of the methods of avoiding the error due to the 
variations of the length of optical path is to obtain absolu- 
te minimum (zero) intensity at the point imder consideration 

10 11 'IP 

along any particular direction (Frocht and Srinath ’ * ). 

Another method is to note the intensities along the required 
directions before and after the loading of the model (Cheng 
The former method has the drawback in that, it would not be 
possible to obtain the zero intensity in general, using the 
methods described previously because of the rotational 

14 1 cr 1 c 

effect. In the method described by Robert and Guillemet ’ ’ 

the variation in the length of light path doesnot enter into 
picture. 

Now, let us recall that (i) the amplitudes along 
directions at inclination to the major axis of a light 

ellipse on either side are equal and (ii) any part of the 
light path in the model is equivalent to a system consisting 
of a linear retarder and a rotator. Now consider the passage 
of light through a system consisting of a linear retarder 
(retar dance p-| , azimuth o(^) and a pure rotator of rotatory 
power . The transformations through this system are 
given by 



117 


p -t 




i S1 


Sl 


S.^Cos 24 ' + Sin2ii'(S2Cos 





- S2Sinp.j ) 

CVi 

m 


S 2 Co 3 f'.j - S^Sino^ ' 


(S^CosP.^ - S^Sinp^j ) Cos2!fJ 




" 

- S.|Sin2it' 

.S3 


S^Cosp^ + SpSin p.| 


^S^Cosp.^ +S^Sinp^ 

- ^ 

I 

r— -1 

- 

- " J 


b1 


=3 


(4.5+2) 


V = (S-j , S2 j is the Stokes vector of the incident 

light referred to the axes of the linear retarder V = (S^ 5 
S2 j ) is the Stokes vector of the light at the exit of 
the system referred to some arbitrary set of axes ^ = (of - ) 

I t 

whereof is the angle between the set of arbitrary axes X, Y 
chosen and the axes of the linear retarder. 

Also , 

- 2 2 

^1 ~ '^4e amplitudes of 

the components of light along the x' , y' axes. 

So, 

— p 2 

a^i - ay’ = S^Cos2i{' + (S2C0SP-1 - S^Sinp^ ) Sin2 4' 

2 


'1 


Consider now, the following cases; 


1 ) 

(a) 

3i 

= S2 - 0, S3 = 1, S'] = 

ax’ 

- ay’ 

= - Sin Sin2 ip 


(b) 


= S2 = 0, S^ = -1, 's-^ = 

2 

®X' 

2 

- ay! 

= Sinp.j Sin2 

2) 

(a ) 

S'] 

= Cos2e, S2 = Sin2p, S3 : 

= 0, S3 = 

(Gos2p Cos2 w 


+ Sin2P Cosp-j Sin2tij) 
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(b) = Cos2(p + 5), S2 = 3in2(P + |), = 0, 

3-j = - Cos2P Cos 2!}< - Sin2P Cos?-j Sin2-p 

3) (a) S-^ ■= 0, S2 = CosCj = SlnG, = Sin2iKCose-^ Cosp^ 

- Sin£ Sinp-j ) 

= Cos(€ + P^) Sin2'^' 

(b) = 0, = Cos(e + 7r), Sin(e + ir) 

S>j = - Sin2if; Cos( £+P-i) 

VJe note that by changing the incident Stokes para- 

meters from the case (a) to those of case (b), the sign of 
the first Stokes parameter of the light at exit is reversed. 
This means that, if the amplitudes along the x' axis is ax' 
in case (a), it 'v/ill become ay' in case (b). This is true 
for any direction of observation since the angle oC is taken 
to be arbitrary. We are thus able to obtain the component 
ayf along the axis X'. This can be achieved in many more 
ways. But the above mentioned ways can be said to be repre- 
sentative of the rest. Of these again, the first two cases 
(1) and (2) are of inport ance since these donot require the 
knowledge of the primary characteristic directions. In case 
(1), right circularly polarized light is converted into left 
circularly polarized light at the entrance. This can be 
achieved by introducing a half-wave plate or by rotating 
the quarter -wave Plate after the polarizer through 90°. In 
case (2 ), the polarizer is rotated through 90°. Using these 
r esplt s.,, we hf.ve.. the,. f.oll-ov|;ing , techniques,.^ 



S - Light source( laser) 

P - Polarizer 
{A/4)- Quarter -wave plate 
EM - Direction of observation 




FIG. 26. SCATTERED LIGHT TECHNIQUE WITH THE MODEL 
ACTING AS AN ANALYSER. 
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(1 ) The Model Acts As an Analyser 

A schematic diagram of the set up is shoxmi in 

Fig. 26. 

Method (a) 

Let the incident light be right circularly polari- 
sed. Note the intensity along the direction of observa- 
tion SM at the point under consideration M. Now make the 
incident light left circularly polarized and note the inten- 
sity I2. Now change the direction of observation until 
= I2. 

If I-| corresponds to an amplitude of then 

2 2 

corresponds "to where ax + ay = 1 . Ifnen, the condition 
= Sy is obtained, it means that the direction of obser- 
vation is making an angle of 4-^° with the major axis of the 
light ellipse at the point under consideration. So the 
final position of the direction of observation will then, 
correspond to one of the secondary characteristic directions 
(for the light path.OM). Now keep the direction of observa- 
tion in this position, rotate the polarizer (after removing 
the quarter -wave plate) until the intensity of light becomes 
zero. The final positions of the axes of the polarizer are . 
the primary characteristic directions. Now, keep the axes 
of the quarter -wave plate 4-^° to the primary characteristic 
directions and the direction of observation 4-5° to one of 
the secondary characterisiiic directicns. Rotate the 
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polarizer until the intensity observed is zero» Then the 
characteristic retardation will be equal to 2 < where is 
the angle betxreen the axes of the polarizer and the quarter- 

43 

-wave plate (inverse Tardy method, Cheng ), This is so 
since the zero intensity (corresponcSing to plane polarized 
light) will be obtained at the point under consideration 
only if the retardation in the linear retarder of the equi- 
valent system is nullified at the entrance. 

The total optical path in the model can be broken 
into small segments so that each segment can be represented 
by a single linear retarder. The axes of these retarders 
can be taken to be the directions of the secondary principal- 
-stress axes at the middle points of the corresponding seg- 
ments. The secondary principal-stress difference can be cal- 
culated from the retardance. of this linear retarder. Consi- 
der the segment AB of the optical path, Fig. 26. First find 
'the characteristic parameters for the optical path OA as 
described above. Keep the axes of the quarter-wave plate 
at to the primary characteristic directions. Keep the 
angle between the axes of the polarizer and quarter-wave 
plate as (v A -°f/2) so that we have circularly polarized 
light (say right handed) at the point A. Note the intensity 
along any direction at the point B. Now change the light 
at the point A from right handed to left handed circularly 
polarized light by rotating the polarizer through 90° a-nd 
then note the intensity Change the direction of obser- 
vation until The final direction of observation 



122 


corresponds to one of the secondary principal-stresses in 
the region AS. To find the retardation corresponding to the 
segment AS, nullify the effect of the equivalent system for 
the optical path OA. Let the characteristic parameters of 
this equivalent system be P-[ (primary characteristic direc- 
tion), P (characteristic retardation) and (rotatory power). 
Now keep a compensator at the entrance with its axes along 
these primary characteristic directions and with a retarda- 
tion of -p . Keep the axes of the quarter-wave plate along 
a direction making an angle of (v/lf - ) with secondary 

principal-stress directions of the segment AE. Fix the dire- 
ction of obsorvaticn along a direction making an angle of 
with these secondary principal-stress axes. Now rotate 
the polarizer until the intensity is minimum. If the angle 
between the axes of the polarizer (in the final position) 
and the quarter-wave plate is oc , then the retardation corr- 
esponding to the segment AB is - 2=C. 

The orientation and difference of the secondary 

principal-stresses can also be found from the characteristic 

parameters using the governing optical equatiors referred 

to some fixed axes, Eqs. (2.4-2) - (2-4-4-). For a particular 

incident light the three Stokes parameters 3-^, Sp and S3 

with respect to some fixed axes are found at each point 

using the characteristic parameters determined experimentally. 

The derivatives S| , 3^ and S3 (with respect to the distance 

z along the optical path) are also found. Then from Eqs. 

(2.4-2) - (2.4-4) Si 

tan2jZS = - ~ 

bp 




FIG. 27. SCATTEREQ LIGHT OBSERVATION ALONG 
TWO SETS OF AXES 45® APART. 
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and 

Si, 

SqSin20 

(Td ) Determination of Characteristic Paramsters From the 
Measurement of Intensities Along Two Directions 
Apart 

Consider any point £ on the optical path in the 
model. Choose two directions which are 45° apart and denote 
these by the numbers 1 and 2, Fig. 27 * Let the incident 
light be circularly polarized. Note the intensities I-] and 
I 2 along these directions. 

I'Je then have 

= iC| af 
Ij = K 2 a| 

Where K are constants which take into account 
1 » 2 

the absorption through the corresponding lengths of scattered 
light Paths. Now make the incident light left circularly 
polarised and note the intensities and I 2 along the dire- 
ctions 1 and 2. 

= It,(1 - af) 

= K, 

I 2 = K2(1 - a^) 

= K2a;^ 


and 



125 


So I-l -f = IC| 


and 



I'Jow conoider the Stokes vector transformation when the 
reference axes are rotated through 4-5°. 



i.e. in this transformation the first and second Stokes 
vectors are interchanged with a sign change. 

Also 

2 2 2 
S. + S 2 = Cos 2co 



Cos^2co 



tan 2k 


(l 2 - I2) (ll + I1 ) 
(I2 + 4 ) (ii - 4 ) 


The characteristic retardation p of the equivalent 
system (for the optical path from the point of entrance of 
the beam into the model to the point P) is given by 
p = ( 2w - Tr/2 ) 
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The orientation of the secondary 

= °C + ir/f . 

characteristic direction 

To find the primary characteristic directions rotate 
the polarizer (with the quarter-wave plate out of field) at 
the entrance to the model until the condition 


g 2 l£l 
s-j 4- S 2 =i 


ii 


.,2 


: ' 2 

|I2 - I2i 


I-l + I 


li 


I 2 l2j 


= 1 


is satisfied. 

The final positions of the axes of the polarizer correspond 

to the primary characteristic directions. Alternatively, 

keep the angle between the axes of the polarizer and the 

quarter -wave plate equal to (v/4 - p/ 2) and rotate them in 
II I2 

unision until — = (IC, and Ko are already known). When 

Ki K2 

this condition is satisfied, the light ellipse at the point 
P is circularly polarized and so the final directions of 
the quarter-wave plate axes will be at 45° to the primary 
characteristic directions. 

In Appendix A are described two more scattered 
light techniques, (i) determination of characteristic para- 
meters when the direction of observation is fixed and (ii) 
a whole field method. 

(2) Scattered Light Technique \^Jhere the Model Acts 
As a Polarizer: 


A schematic diagram of the set up is shOT-Ts in 


Fig . 28. 
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The technique to "be described is essentially an 
extension of the intensity calibration method for the deter- 
mination of characteristic parameters (Chapter 3)* Let the 
angle betv/een the axes of the quart er-uave plate and analyser 
•^e 45'^, Fig. 28. Note the intensity . Rotate the analyser 
through 90° and again note the intensity I 2 . Rotate the 
model or the incident light beam about the direction of 
observation until I-j = I 2 • The final posj.tion of the 
incident beam corresponds to one of the primary cnaracteris 
tic directions. The reasoning is as follo'/zs. 

The intensities I-^ and I 2 are the intensities of 
the components of light vector along two mutually perpendi- 
cular directions. If these are equa.1, this means that the 
light ellipse before the analyser has its major and minor 
axes along the axes of the quarter-wave plate so that before 
this plate the light is plane polarized. The corresponding 
position of the light beam is one of tlie primary characteris- 
tic directions. To find the secondary cho-ractej. ist ic dire 
ctions remove the quarter-wave plate and rotate the analyser 
for minimum intensity. In the final position, the axes of 
the analyser are along the secondary characteristic direc- 
tions. Now keep the light beam and the quarter-wave plate 
axes at 4-5° to the primary and secondary characteristic 
directions respectively. In this case the amplitudes of 
light coro.ponents along the primary characteristic directicns 
are equal and so the amplitudes of light components along the 
secondary characteristic directions are equal at the exit of 
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the scattered light path in the model. The light ellipse 
at the exit has its major and minor axes along the axes of 
the quart er -wave plate. Kow rotate the analyser for minimum 
iiitensity and let the axes of the analyser in this position 
make an angle >d.th the axes of the quart er ---^ave plate. 

The ellipticity of the light ellipse before the quarter-mve 
plate is tanc( , The phase angle along the secondary charac- 
teristic directions is 2 .< so that P = 2-< where p is the 
characteristic retardation of the optical, path Fig. 28. 

The characteristic parameters of the equivalent systems for 
the optical paths Mq M-j , M-j ^ , .... etc . are thus 

found. Using this data and the governing optical equations, 
the difference and the directions of tne secondary principal- 
-s tresses are found as described earlier. 

In this chapter, we have discussed the assumptions 
involved in so-.?e of the previous scattered light methods. 

This is done using the analytic.al solution obtained for the 
case where the rate of rotation is arbitrary, but the ratio 
R = jZS'/P* is a constant and numerical results for optical 
phenomena Tor some loadings of a rectangular model under 
combined tension and torsion. For the cases shown, the 
scattered light methods developed by Frocht and Srinath and 
Swinson et . al. are found to be inadequate. Five methods of 
scattered light experimentation and analysis are described 
which neither involve any assumptions nor the rotation of 
analyser at a particular speed as in the method described 
by Robert and Guillemet . Of these, method 1(b) which 
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involves the measurement of intensities along axes ^5 
apart is a very fast method of determining the characueristic 
parameters at various points, if the computation of one pro- 
ducts and quotients of the measured inrensities is made auto 
matic. As nondestructive techniques these metnods hold high 
potentiality for three-dimensional stress analysis. 



chapter 5 


MISCELLMSOUS TOPICS 


Introduction; 

In the previous chapters methods have been descri- 
bed to determine the characteristic parameters for an opti- 
cally equivalent system either for a series of retarders 
or for any optical path in a photoelastic model under a 
three-dimensional state of stress. Using the concept of 
light ellipse and Stokes parameters, we shall investigate 
three additional topics in ihot oelasticity, namely; A. The 
Graphical Methods, B. The Theory of Compensation and C . 
Determination of Integral Fringe Orders In Photoelasticity. 
These topics are useful in the general three -^mensional 
photoelastic analysis. 

A. Graphical Methods in Photoelasticit y 

Graphical methods in photoelasticity are essentia- 
lly pictorial representations of the changes in the para- 
meters of the light ellipse as it passes through a photo- 
elastic model, or the polarise ope elements. These are 
useful (1) in finding the changes in the light ellipse 
parameters as the polarized light passes through various 
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polariscope elements, without cumbersome calculations, 

(ii) in finding the characteristic parameters of a parti- 
cular light path in a three-dimensional photoelastic model 
or for a series of linear retarders and (iii) in estimating 
the effects of various errors in the elements of the polari 
scope (ex'i error in the orientation of the axes of the pola 
riser, deviation of the retardation of the quarter -wave 

plate from -rr /2 etc . ) • Most of the existing graphical 

23 

methods are derived from the Poincare' sphere (see 
Chapter 3) representation of polarized forms. After a 

brief description of the graphical methods described by 

94 14,15,16 25 26 

Jerrard ^ , Robert and Guillemet , Cuske , Schwieger 

27 

and Cernosek , we shall describe a new and simpler appro- 
ach to graphical methods. 

Jerrard^'^ has given a detailed account of Poincare 
sphere representation, its derivation and plane representa- 
tion of the polarized light by using comiorrnal transforma- 

14,15,16 

tion. Robert and Guillemet also have discussed these. 

In the derivation of the Poincare' sphere representation 
the light vibration is first represented in a complex 
plane. The most general type of polarized light can be 
represented by a combination of two simple harmonic motions 
along the axes OX and OY as? 

X = a^ Cos (St = Real Ae^^ (5»'i) 

with A = 




FIG .29. REPRESENTATION OF POLARIZED LIGHT FORMS IN 
A COMPLEX PLANE. 
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Y — S-Y Cos ( 6jt; 4- = R03.1 Bs ' ' ( 5" * 2 ) 

• -t-i " iGy 

witn 15 = a,^e ^ 

5;here w is the angialar frequency, t is the time and 
Gy are the phase angles. 

The ratio B/A is in general a complex quantity and 
can be -/.rritten as; ' 


VJ = B/A = u + iT 


ay ” ^ 

3-X ^ 


where 


t an e 



= tan 0 (Cos 0 + i Sine) 

(5.3) 


tan 6 is the ratio of the amplitudes of the components of 
vibrations along the reference axes OX and OY, Fig. 29, and 
0 is the difference between the phase angles of these 
components. So, any general type of polarized light can • 
be represented by a complex quantity B/A = u + iv as des- 
cribed above. This is represented by point M in the com- 
plex plane u - v, Fig. 29(a). 

OM = -^ Amplitude ratio 
^X 

The inclination G of OM with the u axis is the 
relative phase difference. If this phase difference is 
increased bye , the point M moves to- m' where OM = OM' 
and MOM' = £ . In this mapping, all points on the real 
axis represent linear vibrations since in this case, the 
phase angle G (inclination of OM with the u axis when M is 
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on the u axis) is zero. The length of OM again represents 
the ratio of the amplitudes along the reference axes CK and 
OY. If B/A is purely imaginary, the representative point M 
lies on the imaginary axis, the phase difference £ is F/2 


and so the major and minor axes coincide uith the reference 
axes OY and OY. (Because, with respect to the major and minor 
axes the phase difference is ir/2,) Right hand vibrations 
in which the OY vibration leads the OX vibration, are repre- 
sented by points above and left handed vibrations by points 
below the real axis. Thus the amplitude ratio and the phase 
difference of a particular light ellipse referred to a par- 
ticular set of axes, can be represented in the complex plane. 
It will however be useful if the ellipticity and the azimuth 
of the light ellipse (vm.th respect to a particular set of 
axes) are also represented in this mapping. For this consi- 
der an elliptic vibration ’-dth ellipticity b/a and azimuth cc 
(with reference to the axes-X,Y), Fig. 29(b). This can be 
defined in terms of the vibration along the axes u,v ( i.e. 
the major and minor axes are along the u,v axes), so 
u = a Cos ot 
V = b Sin oTt 
b/a = e = ellipticity 

Resolving along the axes OX and OY, we have, component alcng 
the OX axis = x 


(a CoS Mt Cosc< - b Sina3 t Sinoc) 


= Real 


(a Coscf + i b Sino() e 


icut^ 


(5 ) 
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C omponenu alcfng^ _ j Cos cot Sino( + b Sin co t Cos=<) 
uhe Ul 3.yjLS j 


= Real I (a Sin=( - i b Cog=c) e 


icot 


(^.5) 


)i:r, paring these equations ■•dth Ens. (5*1) - (5«3), we obtain, 


a Sin- - i b Cos=( 
a Cos3( + i b Sin=< 


tan-~< - i e 

1 -r i e tan'a' 

This equation is a special form of the bilinear transformation 

( 5 . 6 ) 

c + dt 

with ad ^ be 

(condition for conformality) 

In this conformal transformation circles and stra- 
ights lines are mapped into circles and straight lines. 

(a) Mapping of Light Ellipses of Constant Orientation 
(’^= Constant) and Varying Ellipticity (- “^ < e < 


w = - 1 e 

1 + i e tanc< 

= u + 1 V 

2 2 
- iMfL (1 - e ) - ie (1 + tan cQ 

(1 + tan^rf) 

so that u = (1 - 

1 + tan-^^f 

u 0 as o( -» 0 or 'ir/2 (5.7) 



:v 



FIG. 30. CONSTANT ELLIPTICITY (e) AND CONSTANT 

AMPLITUDE RATIO CIRCLES IN THE CONFORMAL 
MAPPING. 
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2 

and V = - Q 

(1 + tan^ 

T ■* - e or-1/e as <=( -* 0 or ■Tr/2 

ilov/ consider two values of e i.e. e^ = e and 
Qj then for these values (with same ), the values 
of u are the sa:ae ’-here as the values of v are equal in mag- 
nitude out opposite in sign. So, the raapping is synsnetrical 
about the u - axis. A.lso u varies with e so that it can 
take nonzero values also (if°f ^ 0, v/2). Therefore, the 

mappings of constant azimuth =( light vibrations are circles 
with centers on the u - axis, (-.hen ’^=0 or v/2, this dege- 
nerates into a straight lines., the v - axis), llso u = 0 
and v = + 1 ife='t.1 irrespective of the values of oc . 
Therefore, all the circles pass through the points P(0,1 ) 
and P' (0, - 1 ), Pig. 30. 

(b) 'Mapping of light Vibrations of Constant Ellipticity e 
and Varying azimuth o( 

Consider two values for tan ^ in Sqs. (5»7) i*e., 
for tan=<-i = f and tan =<2 = - f (e is same for both). For 
mappings of both these vibrations the values of v are the 
same, but tbe values of u are equal in magnitude but oppo- 
site in sign. Hence these mappings are syCTetrical about 
the V - axis. Also v varies with if e ^ 0 so that it 
takes nonzero values also. Therefore, these mappings are 
circles viith their centers on the v - axis, then e = 0, 
this degenerates into a straight line, the u axis. 




These two families of circles are 'orthogonal 
since OF. OP' = OH. Oil' = 1. Thus, if a vibration is represen 
ted 'ey a point M the oriontation of the ellipse is obtained 
from the circle throiogh P, ?' and M with OQ = tan=c, and the 
ellipticitj is obtained from the orthogonal circle through 

H v/ith 8 = OH. 

(ii) ivasks’s^'^ J Circle Methods 


Consider a plane polarized light of unit amplitude 
represented by the straight line AOA as shoc/n in Fig. 31* 

The components along the principal-stress axes p, q are Cos 6* 
and Sine’. If a phase difference oft' is added between 
these two components in the model, 

Vibration component along) 
p axis ) 


Cos e' Cos oj t 


Vibration component along 3 _ 

j) = Sin 9 CosCoit +0 ) 
q axis ) 


Now '".'ith respect to the axes OE and OC which are 
inclined at an angle Y to the p, q axes, the vibration com- 
ponents are given by 

Component along OB = Cos^t (Cos© Cos y+ Sin© Cos 0 Sin y) 
- Sin9* Sin^t Sin0 Sin y 

= a ’ Cos (“t + 0^ ) 

Component along OC = Cos“t (Sin©* Cos £Cos Y 

- Cos9 SinV) - Sin w t Sin©^ Sin0 Cos y 
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and 


Sinks’ Sin^Y + Cos^e' Cos^v + Srn2e' lose Cosy Siny 

(5.8) 

Sin'^©' Cos^Y 1 - Cos'^e' Sin^Y - Sin2e' CosC CosY Sin y 


(5.9) 


o find tie maximun and 'Tiinimiiin values of 


and b' (corres- 


ponding to the lengths of seoii -major rod minor axes a and 
b), ore differentiate the above expressions vith respect to 
Y and eqeiate to zero 


,2 ,2 
da db 
dY - dY - 


giving tan2Y = Cosetan2§’ (5»'10) 

Substituting this expression in terms of y in Sqs. (5*8) 
and ( 5 * 9)5 obtain 


a2 _ ^ j 1 + (Costae’ + Cos^e Sin22e* )2j (5.'>'l ) 

=: 2" (1 J ) 

where <1 = (Cos22e’ + Cos^e Sin^2e' )^ 

and b^ = i (1 - J) (5-12) 

The quantity (003^2 0* + Sin^ e’ Cos^e is called 

J. The ellipticity and the orientation of the major axis 

can be represented by an arrow (vector) having length equal 

to J and oriented at a chosen angle. The angle between 

this direction and a chosen reference direction are doubled 

and so are all the other angles involved. 

Now, <1 = 2a^' - 1 =1 - 2b^ = 

1 + e‘ 

• vrtiere e = § = ellipticity 

St 


(5.-13) 
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So, the magnitude of J depends on the ellipticity 
only. Further, 


e2 = ., 1 . I - 4 (5-1^) 

1 + J 

So, J = 1 represents plane polarized light and 
J = 0 represents circularly polarized light. That is in 
Pig. 32, 3- tip of the 1 Tector on the circumference of the 
circle necins plane polarized light, a tip in the center 
corresponds to circularly polarized light while all other 
positions represent ellipticaHy polarized light. Since 
all the angles are doubled, the mutually perpendicular 
principal-stress axes are represented by the straight line 
qOp, since qOp = ir . If plane polarized light is at an angle 
6* to one of the principal-stress axes, the components along 
the principal-stress axes are Oos e' and Sin 0* . Consider 
the plane polarized light represented by the point A on the 
c ircumf er one e of the circle. The line AC is perpendicular 
to line qCp. 

Now, 

q B = 1 - CosZe' = 2Sin^e 


and BP = 2Cos^e' 

i.e. q B and BP are proportional to the squares of the amp- 
litudes of components along the principal-stress axes. In 
order to find the intensities of light components along two 
mutually perpendicular directions', a perpendicular is to 
be drawn from the tip of the J vector on the straight line 
which corresponds to these mutually perpendicular directions. 
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The intercepts correspond to the required intensities. Since, 
during the addition of retardation along the principal-stress 
axes, the intensities of components of light along these 
axes remain constant, the representative point moves along 
the straight line AC which is perpendicular to the line qOp- 
At the point A, 

j = (Costae' + sin^se' cos^e)^ = 1 

or Cos^^ = 1 

i.e. C = 0, ir, . . . . , nir 
At the point B, 

J = OB 

= - Cos2©‘ 

= (Cos^a©' + Cos^ Sin^2©* )r 
or Cos^£ = 0 

, , (2n -f 1 V 

i.e. £ = ir/2, 3*^/2 , 2 

This can be interpreted as follows; At the point 
A, the phase difference is zero. As the phase difference 
is increased toir/q, the representative point moves to B 
along line AB. As the phase difference is further increased 
t B TT it moves to the point C along line BC . When it is 
3^/2, the representative point comes back to B and so on. 
(lit) Schwieger's Analysis; 

Schwieger presented some rotation rules to be used 
•vfith J circle method. He introduced the Wulff ’s grid method 
of crystaUc^raphy into photoelasticity and described methods 
to find the characteristic parameters for a combination of 




parallel pro JECT fON' 

POLARiZ ATION (7 C(ROgMaifti|fei^ 
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linear retarders- A brief suimary of his analysis is pre- 
sented below; 

The four parameters of light ellipse, namely, <=(. , 
the azimuth, the angle G (or tan 9, tne amplitude ratio), 
the angle w (or tan o , the ellipticity) and G the phase 
difference. From the equations of Chapter 2 we have the 
following relations connecting these four pa,rameters. 

SinC Sin2G = Sin2o (5-1?) 

Cos€ tanpe = tan2c< (5-1^5) 

C os2°< Cos2 o = Cos2 9 (5 *"17) 

The angles of a spherical triangle satisfy these 
relations and so these parameters can be represented on 
Poincare’ sphere. 

(a) Parallel Projection: 

In Fig. 33 . all points of the upper hemi-sphere 

represent left elliptically polarized light while points in 
the lower hemi-sphere represent right elliptically polarized 
light. The area of the J circle can be obtained by parallel 
projection on a plane having a normal which is parallel to 
the axis connecting the poles of the sphere. 

The chords drawn in the J circle, Fig. 3^* are 
lines of constant amplitude ratio, the ellipses are loci of 
points of constant phase difference. The radial lines are 
lines of constant azimuth. The actual azimuth will be half 
of the angular inclination of these radial lines with the 
OX axis. The refer ehee axes are the lines C3K and OY. The 
scale on the right hand side of this plane representation 
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FIG. 3 5. STEREOGRAPHIC 
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is calibrated in terns of ellipticity. This plane represen- 
tation can be readily used to find the resulting photoelastic 
effect after the Passage through a set of linear retarders. 
Ad'dition of phase difference corresponds to moving along the 
chords of constant amplitiiide ratio, the length of inovement 
corresponding to the value of retardation added. Rotation 
of reference axes corresponds to moving a.long a constant 
ellipticity circle until the angular inclination to OX is 
increased by double the amount of rotation of reference axes, 
(iv) ’i-.uff's Grid! 

This is the stereographic projection of Poincare' 
sphere of polarization. 

The stereographic map is the projection of a sphere 
on the plane of a great circle, Fig. 35, in such a way that 
the projection of a poout A cn the sphere is the point in 
which she line through A and a pole of that great circle 

intersects tlie plane of the great circle, iliis projection 
44 

is conforrrral 

The chords of constant amplitude I’atio in the para- 
llel projection correspond to the curved lines surrounding 
the YOX axis in the stereographic projection. The radial 
lines here again represent lines of constant azimuth. The 
loci of points of constant phase difference are arcs of 
circles in this projection. 
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Determination of Parameters of the Eqiii valent Optical 
System For A Combination of Two Linear Retarders-: 

= Retardation introduced in the first plate 

= Retardation introd’-iced in the second plate 
a 

The optically equivalent system consists of a linear 
retarder and a pure rotator. If the incident light is plane 
polarized along the axes of this linear retarder it remains 
plane polarised after the retarder but is rotated by the 
pure rotator. That is, it is plane polarized at the exit. 

To find the rotatory power of this pure rotator and the orie- 
ntation of this linear retarder draw two lines QX^ and CK 2 
corresponding to the axes of the two plane models. Cons- 
truct the circular arcs on these lines corresponding to 

= constant and 180 - A 2 - = constant. From this point of 
intersection P of these two arcs, proceed to the point Q and 
R on the circumference of the circle along the constant 
amplitude ratio ciurves corresponding to the axes QX-] and 
OX^* The angle ROQ will be twice the rotatory power of the 
system. Also, X-| OQ = 2pj,gp]i_, where is angle bet- 

ween the axes of the first plate and the axes of the linear 
retarder in the optically equivalent system. The procedure 
can be easily understood from Fig, 37. Ir this figure 
Ai P 2 A 2 represent s the passage of light through the two 
retarders.: 

P^^P^ = 20 , 

The actual problem is to locate the points P^ and 




FIG. 36- DETERMINATION OF THE ANGLE 9 ACTIVE 
OF THE AG TIVE 

SVSTEM FOR TVV0’':^L|litMSiEL^^^^^^^^^ . ^ 
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P_ which lie on the constant phase difference curves & zz 
and A = 180 - respectively such that 

P^OP^ . 26 

This is done by rotating one curve '.ath respect to 
the other b;/ an angle 20 as in Fig. 36. so that their 
point of intersection P corresponds to the points P^ and P^ 
of Fig. 37. Fig. 33. shows the orientation of the axes of 
the re'carders. 

In Fig . 36 . , 

R^Q = 2 (e^ - ) 

A 

ROQ f as mentioned ^y Schwieger. 

To find the retardance of the optically 

equivalent system consider plane polarized light along the 
axes of the first plate represented by the point P in Fig. 
39 . Retardation A^ in this plate brings it to P^ . Rotation 
of the reference axes to make them coincide with the axes 
of the second plate brings it to the point P^- Addition of 
a retardation A^ in the second plate is represented by the 
path P 2 P 3 so that the point P^ represents the state of pola- 
rization at the exit of the tx^fo plates when referred to the 
axes of the second plate. We should arrive at the same 

point P_ even if we use the 'optically equivalent system. 

■ ■ ' ' , ■ ^ 

low proceed back to the point where ^ 36 ^ 14 . = 

and Oli, = OP^:. Find the point P^ on. the circumference such 

that pfe^ = 26^ . That is, -the point, P^ represents the stats 
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of polarization at the entrance vhen referred to the axes 

of the eq‘'jiiTalent retarder. The points and P, lie cn 

5 M- 

the same carve of constant amplitude ratio. (This fact can 


be used to check the values of and The path 

P^P^_ represents the equivalent retardation In the 

case of systems with many number of retarders, this method 

can be used successively. Schmeger, in his example with 

ten plates, mentions an iterative procedure where the azimuth 

of the incident plane polarized light is varied until the 

light at the exit is plane polarized. 

(vi) Graphical Method Using Quaternions 
* 27 

Cernosek devised a geometrical method based on 
the properties of quaternions to find the characteristic 
parameters of a system of linear retarders. He made use of 
Jones vector in this derivation. The principle of this 
method is described below. 

Let the Jones vector for a certain elliptically 
polarized light referred to the X,Y axes be given by 


i(£^ + w t) 


a e 
T 


i(eY +(i: t) 


where phase angles and are 


the amplitudes of the cemponents along the X, Y axes 


respect iyely. ^ 

,1 0 

Let , I = 

10 1 


i 0 
0 -i 


; 3 = 


0 1 

-1 0 


k = 


0 i 

1 o 


( 5 . 18 ) 




FIG- 40 - CO-ORDINATE SET USED FOR 

THE DESCRIPTION OF THE GENERAL 
QUATERNION Mi, M2, M3, M4 ARE 
THE COEFFICIENTS OF THE QUATERNDN. 
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The genera.! unitary 2x2 matrix, which represents 
the tr'Uisf orraation of polarised light from its entrance th- 
rough a .s:v’'stem of linear retarders to its emergence can be 
iiniquely expressed as a quaternion 

M = M.^ I + i + 3 + k (5*'!9) 

The coefficients of the quaternion M are connected 
by the unitary nor.m relation 

lif + M2 + M3. = 1 (5 *20) 

The quaternion can be described graphically as 
shOv'.na in Fig. 4-0. The co-ordinates of the four points 
K-j K 2 K 3 and in Fig. 4-0. are the coefficients of the 
quaternion M defined by Eq. (5»'19)» 

■and M 2 % respectively in four different co-ordinate systems 
(m m ), (m m ), (m m ) and (m m ) which are in the same 
plane but have different orientations. The values of ^^ 2 , 
M 3 and for a rotator and a retarder placed before or a 
after the system can be easily found from the transf ormticn 
relations. Bach of those transformations corresponds to a 
particular type of rotation of some of the four radius 
vectors OK^ OK.^ OK^ and OKj^. Thus, the final positions 
of the points K.j ^2 ^^3 can be found graphically very 

easily for either a retarder or a rotator when it is placed 
either before the system or .after the system by simple rota- 
tion of the radius vectors and-S^ ^-5 and K^. Two 
other points in each step hwill then be found by drawing 
lines parallel to the coordinate systems. 
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The Characteristic Quantities of a General System of 
tietarders? 

Before studying a system of retarders, ve shall 
preso't the quaternions for simple operations. 

Quaternion of rotation: 

S (^p) = Cos=<2- I - Sirf^g* ^ 

is the rotation angle. 

Quaternion 'Uc<„ = S M when S (cCp) is placed after M. 

2 a ^ 

then 

U = A-j^l '^ 32 ^ ■'^ 2 ^ (5*21 ) 

where 



= 

CoS<^ 

+ I'L 

Sinpe^ 

(5*22) 

12 

1 

2 

3 

2 


OJ 

= Mg 

Co^2 


Sinp<2 

(5*23) 

A 

= M 

Cos°C 

- M 

SinoC 

(5*24) 

32 

3 

2 

1 

2 


\2 


Cosoc^ 

"«2 

SinoC^ 

(5*25) 


Quaternion for a linear retarder of retardance gi'ven 

by 

If this retarder is placed after M, the quaternion 
for the combined system is given by 

U ( A 2 ) = G (^2^ ^ ^ B'jpi + ■^22^'*'®32^ ^4-2^ (5*26) 

B = CosA^ 

12 1 2 


where 


- M SinA^ 
2 2 


(5*27) 
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C os oinA^ 


CosA^ - SinA^ 

3j_^2 = \ CosA^ + SinA^ 


(5 * 28 ) 

(5.29) 

(5.30) 


If the rotator and. the retarder are placed before the system 
M, the corresponding quaternions for the combined systems 
will be as given bolow? 


U (°^ Y ^ = M S (°^q ) = A-] -j I + A 2 '{ i + -^31 ^ 
where 


and 

U^^(= M G (A^ ) = B^-jI + B|^i + j + Bj_^^k 
where 


(5.31 ) 


"Si 


\ 

0 

0 

.M 3 

Since. 

1 

(5.32) 

■'^21 

= “2 

003“=^ 

1 

. 

Since 

1 

(5.33) 

A31 

= M3 

CoSoC| 

r M, 

Sirb(-| 

(5.3^) 


= \ 

cos^f 

1 

-“2 

Sirff., 

1 

(5.35) 


(5.36) 


hi 

^1 

COSA^ 

1 

-«2 

SinA_, 

1 

(5.37) 

®21 

II 

C os A 

1 

+ M., 

1 

Sin A_, 

1 

(5.38) 

B 

31 

II 

Cos A 

1 


Sin A 

1 

(5.39) 

^1 

= \ 

COSA.^ 

-M 3 

SinA| 

(5.^0) 


The coefficients of the resulting quaternion TJi 
after a rotation through an angle found by 
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FIG. 43. 


THE CHARACTERISTIC OJANTITIES OF A 


GENERAL SYSTEM OF RETARDERS Wj ARE 
THE ORIENTATIONS OF THE PRIMARY AND 
SECONDARY CHARACTERISTC DIRECTIONS 


2p IS THE CHARACTERISTIC PHASE 
RETA^W-,. ,, , 

•" tan, 

^ s . 

tan fW- Q> = J_!^ . . 
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rotating the radius vectors of the points and in the 
clockwise direction, Fig. 41. as can be seen from Eqs.(5*22)- 
(5*25). In the same way the coefficients of the quaternion 
for a retarder whose axes are inclined Zit an angle oc, with 
the reference axes can be found by rotating the radius vec- 
tors of and (multiplication by the quaternion of 
rotation) through an angle in the opposite directions 
and then rotating the radius vectors of and by an 
angle A-| , both in the same direction (for the addition of 
a retardance of A-j in the linear retarder), Fig. 42. 

Similarly, the geometrical procedures correspond- 
ing to the equations (5*32) - (5*35^) a-nd (5*37) - (?.40) 
can be established. 

Using .the above relations, ve have the following 
representation of the quaternion for the combination of a 
linear retarder and a pure rotator (of the optically equiva- 
lent system). Fig. 43. Following the procedures described 
above, the quaternion for a system of linear retarders can 
be found graphically. Comparing this with the quaternion 
for the optically equivalent system. Fig. 43. the characte- 
ristic parameters of the system of linear retarders can be 
found. In Fig. 43. the characteristic quantities of a gene- 
ral system of retarders are given by 

■ : .ei + p 

'^ 2 ' . 
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So, the Tariation of any two of these parameters with respect 
to one another can be plotted keeping one of the other two 
parameters constant. For example, the azimuth cc and the 
phase difference 0 can be plotted along X and Y axes so that 
we have constant amplitude ratio curves and constant ellip- 
ticity curves. We will have other systems of graphs if 
other combinations of the parameters are taken with different 


coordinate systems. 

From the above equations, we have 


and 


tan€ 


_ tan2t>J 

- Sin2P( ~ (1 _ e 


Say- ay 

tan2=< = — 2 

^ - a^ 


2e 

S3n2^ 


(1 -e|) 


Cose 


where 

e = b/a = elliptic it y = tan tu 
ej= = amplitude ratio = t^e 


(SM) 


From these relations, it is clear that for given 
values of e and , the values of C is fixed and for given 
values of e^^ and € , the value of is fixed. The values of 
these x)arameters can be tabulated and this table can be used 
for computing the effect of a series of retarders*' 

From Eqs. (5.^^) and (5 •^5), note that for a 

given value of e the eULipticity (not equal to unity), the 
phase difference € has a minimum value which occurs at 
oC — For a given value of ey the amplitude ratio (not 

equal to unity), the azimuth «< has a maximum (or minimu m) 




FIG. 44. THE FIRST GRAPHICAL METHOD 
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■valu 0 whicli occurs at^ = 0 . wTion the aupli'nude ratio 
equals unity and. the phase difference £ 7 : 90'"'^ > th® azimuth 
equals hlien e^ = 1 and£ = 90 ^^ it is circularly 

polarized light and takes all values. 

One can plot Eqs. (5*^*+) and (5*^5) as shown in 
Pj_g. 44 .. In this figure, the horizontal axis represents 
the azimuth and the vertical axes the phase differencee . 
Use of Sq. (5*^) gives ’’Constant Elliptic it y" curves. To 
get these, one plots vs C ; for fixed values of e getting 
curves like APB, AEB, etc. Use of Eq. (5-^5) gives ’'Constant 
Amplitude Ratio” curves. These are obtained by plotting 
=< vs ^ for fixed value of At points A and B the liiase 

difference € is v /2 and the azimuths =<are 0 and v /2 respec- 
tively. So points A and B represent states of polarization 
when the reference axes coincide with the major and minor 
axes of the light ellipse. The line represents circu- 
larly polarized light and the lines AO, OC and CB represent 
plane polarized light. This form of graphical representa- 
tion of light ellipse is very caivenient, since in photo- - 
elastic analysis, we add a certain amount of retardation 
alcng a particular set of axes (during this addition, the 
amplitude ratio along these axes remains constant and so 
this process is represented by movement alcng the constant 
amplitude ratio curves) and rotatim of the reference axes 
(during this process, the ellipticity remains constant and 

so it is represented by movement along the constant exlip- 

ticity curves). The same graph can be used for ranges 




OR- Axis of the first retarder 
OR2--Axis of the second retarder 
OP -Reference axis 


FIG- 45- ORIENTATION OF THE AXES OF 
THE TWO RETARDERS. 
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beyond Tf /2 radians for the azimuth and the phase difference 
by taking into accoijnt the sign variations and using Eqs. 

and This is demonstrated at the and of this 

section. 

Procedure to Use the Gr a phi cal Method; 

Consider two retarders, and whose axes are 
inclined with respect to each other by an angle ©2 1®"*^ 

the axes of *make angle 0^ and ©-[ + 7/2 with respect to a 
given fixed axis, Fig. . Let the initial light be plane 
polarized along the refering axis. This initial light elli- 
pse is represented by point 0 (with € = 0 and°f = 0), Fig. 

^6. When this is referred to the axes of , the initial 
light ellipse ’/ill be given by point (with e = 0 and amp- 
litude ratio given by point ). Addition of retardation 
by is represented by movement along constant amplitude 
ratio curve P-j P 2 up to the point P 2 , such that the ordinate 
(phase difference G ) of P 2 is greater than that of the point 
P.^ by an amount C^. Now, we refer this light ellipse to the 
axes of the second retarder by moving to the point P^ frcm 
the point P 2 along the constant ellipticity curve such that 
the azimuth of the point P^ is greater than that of the 
point P by an amount ©-, where © is the angle between 

2 €L. 

the axes of the two retarders. Addition of a phase differe- 
nce of €2 in the second retarder is represented by movement 
along the ccostant amplitude ratio curve P^ Ilj. from the point 



F(G.47. VARIATiaN OF LIGHT ELLIPSE PARAMETERS 
IN A PLANE POLARiSGOPE. 


I 

I 
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to the point where the ordinate £ of the point 
is greater than that of the point hj an amount £2? "where 
^2 is the retardation introduced by the second retarder. 
Finally, the point Pj^ in the chart represents the light 
x'^ibration at the exit of the second retarder with respect 


to its axes. So, during the representation of the v.aria- 
ticn of light ellipse as it passes through a system of li- 
near retarders on this graph, the point under consideration 
always represents light vibration with respect to the axes 
along which the retardation has just been added or about to 
be added. 


Passage of Light through Elements of a Polariscope: 


(i) Plane Polariscope: 

In Fig. 47* point 0 represents the plane polarized 
light along the axis of the polarizer 5 point P^ , the plane 
polarized light referred to the principal-stress axes of the 
model at "its entrance; point P^j the light at- its exit after 
an addition of £-j retardation referred to the axes of the 
model, and P^, the same light vibration but referred to the 
axes of the analyser which is kept cro ssed with respect to 
the polarizer. If one of the axes of the retarder (i.e. 
model) coincides with the direction of the initially pola- 
rized light, then we move from 0 to P2', where OP^ = ♦ 

When the analyzer is at ir/2 to polarizer axis, we move along 
the constant ellipticity line .P2OLC (e =0) until we come 
to the orientatlcn of the analyzer axes, that isjPOiat G, 
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The amplit-ude ratio is zero and no light emerges 
analyzer . 

(ii) Circular Polariscope: 

Fig . 4-8 shous the elements of the circular pola- 

riscope. The- rcints P. and in Fis". 4-9. represent light 

'I ^ ^ 

vibrations before .rad after the first quarter-wave plate 
referred to its .axes, she points P 2 , P^ ,the light vibra- 
tions before and after the model referred to its principal- 
-stress axes at the particulxr point under consideration. 

The points Ilj., represent the light vibrations before and 
after the second quarter-wave plate with respect to its axes 
and the point P , , the light vibration before the analyser 

D 

referred to its axes. Now, we have to prove that ©x-i » 

amplitude ratio corresponding to the point P depends oily 

6 

on the retardation introduced in the model and not on the 
angle 0, the orientation of the principal-stress axes with 
respect to the axes of the quarter-wave plate. In a graph 
of Fig. 49 ,. which is drawn to scale, point P^ will always 
lie on the line, irrespective of the value of 0 

JL 

used. This can also be shown by Eqs. (5»4^) (5*^5) 

that are used to construct Fig. 49- To show this we proceed 
as follows: 

Frcan Eqs. (5»44) and (5»4-5), we obtain 

tane = (5.^6) 

Sin2c( 

and tan2=< = tan2 0CosC (5*4-7) 
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where tan© is the amplit’ide ratio. 

Let e 3 ,£k, € 5 , 06 hG the phase differences at the 
points ?3 respectively. 

Prom Fig . 4-9 

£3 = ^/2 - e* 


PC _ TT/l^. 

3 “ ' 


tan-, o = _ tenpoi 

Sin2 c 3 


2^)3 = = ^/2 - e 


For the point P, 


4 


tan^i^ = 


t an 2 0 ). 


3 


TT 


Sin 2U + 0) 
4 


= Coig 

Cos20 


e ^ = 0, - /2 

5 4 

0 . ^ C os20 

. . tan 0;- = - Cott), = ^ 

5 ^ Cot0* 

from Eq. ( 5 - 5 ^ ) 

Also, for points rlj. and P^ 


tan 2 (jl + 0) = - Cot20 = tan2ei^ Coaei^ 
and tan 2°C^ = tan 26^ Cos 0^ = tan 2e^_ 3in0 4 . 

. . tan 2 ^^ = - Cot 20 tan0 ^ 

_ _ Cot 20 Cot0* 

Cos 20 

from Eq. (5 .49 ) 

' _ _ Cot £* 

Sin 20 


(5.^8) 


(5.^9) 


(5.50) 


(5.52) 

(5.53) 


( 5 - 54 ) 
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Nov/, for tile points and we have 


and 


tan 2iu. = tanc ^ Sin 2=c_ 

555 

tan 2“^ = tan 2^^ 


i5-5y) 


= tane ^ Sin 2=c^ 


tanG.^ Cos 


(5.56) 


van S 

D 


tan G^ tarpC2 ^ 

5 5 


_ IT /1+ 

6 - 5 


= - Cot 20 


from Eqs. (5*50) and (5*5^) 


= ' n "/2 + 20 


(5.57) 


[■low for point ? 


tan 20> 


tan Cotp<2 


5 


* CosG^ Sin 20 

from Eq. (5 *57) 

Sin 20 


Cot G* Sin 20 

from Eq. (5 *5^ ) 
= tan (ir-G*) 


or 


2Q^ = (it -G ) 


% 


1 I 

2 "2 


(5-58) 



179 


So, the amplitude ratio tan© ^ = 
retardation € introduced in the model and not on the orien- 
tation 0 of the principal -stress axes -^ith respect to the 
axes of the first quarter -wave plate. In the case of light 
field (the quarter-wave plates are crossed but the polarizer 
and analpsoi- arcj: parallel) polarise ope, 


depends only on the 


6. = nv 

whenever C* ^ = (2 n + 1 ) v 
then 


tan © 


6 




0 


i.G. the extincticn of light occurs. 

In the dark field set up, the analyser and polari- 
zer are crossed and the two quarter-wave plates are crossed. 
In this case 

Here Eq. (5 *56) becomes 

tan = tan Cos 

and tan = Cot 20 

= tan (I - 20) 

2 

or ^6 = - 20 (5.57) 

becomes 

®6 = 




FIG. 50. OPTICAL SET-UP IN TARDY'S METHOD. 
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so that 

tan = tan£'*/2 
6 

0 


i.e. extinction of light occurs. 

•vhen - 2ffir 

Those results cm bo proved in exactly the same way 
for other al.ternative set ups of the dark and light field 
arrangements. 


(iii) Tai'dy's Method; 

The set up is shoim in Fig. ^ 0 . Point 0 in Fig. 51 


represents plane polarized light referred to the polarizer 
axis. Points 1 and 2 represent light before and after the 
first quart er -wave plate referred to its axes. Points 3 
and 4 (path 3-2-4-) represent light before and after the 
model ivLth respect to the principal-stress axes. Points 5 


and 6 represent light before and after the second quart er- 
-wave plate. In order to cut off light, the amplitude ratio 
along the analyser axis must be oo , This can be achieved 
by making point 6 come to point 7 along the constant elli- 
pticity line. The orientation of the analyser axis is then 
_ o(^) vrhose value can be found as follows: 

Fran Fig. 51 we note that 

= 7r /2 - e* 




tan 2 




tane Sin 2 ^^ 
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= tan 

h 

= tan CL - C* ) 

3|e 

or 203^^ = Tr/2 - C 

Als Oj 

= Tr/2 

From Sq. (5.^3) 

Sin 2© SinG = Sin 2 a) 
when G = Tr/2 

6=0) 

. . 6^ = 0)^ = 

tan 2=<^ = tan 2©^ Cos = tan 2>i^ 

3fC 

i.e. c<6=o>^=u)^=|-| 

The angle between the axes of the second quarte'. 
-wave plate and analyser is given by 

=<7 -«6 = W 2 -«6 

. . The orientation of analyser axes is given by 
5 = (c(^ - o(^) + ttA 
= 3 ir A 



' or €*-,^ - -r 

From this graiSiical representation the essential 
features of light ellipse can be cfosei*ved clearly. For any 
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light ellipse, for' directions nearly at 45*^ to its major or 
minor azis, the -variation of phase difference vith direction 
is very small. (Because, the constant ellipticity curves 
are alm ost flat for values of =< around 45° ) . This variation 
is high for directions close to the major or minor axis. 

Also, 'TO see that as t'he ellipticity nears unity, the varia- 
tion of amplitude with direction hecomes small. This chart 
will be very useful in finding the effects of small errors 
in the retardances of quarter-wave plates readily. 

Determination of Characteristic Parameters: 

To find the characteristic parameters of a system 
of linear retarders Schwiegsr suggested a geometrical cons- 
truction and an iterative method, and CJernosek presented 
a method based on quaternions. A much simpler and direct 
method to find these characteristic parameters will now be 
described. If the incident light on the equivalent optical 
system is circularly polarized, the amplitudes of the compo- 
nents of light vibration before and after the linear retar- 
der along its axes (primary characteristic directions) are 
equal, and these components are just rotated in space by the 
pure rotator to coincide vith the secondary characteristic 
directions (Because, the angle between the primary and sec- 
ondary characteristic directions is equal to the rotatory 
power). So, the mut-ually perpendicular directions at the 
exit of the system along which the light components have 
equal amplitudes are the secondary characteristic direc-ticns. 
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This mc-thod can be used with any of the above described 
graphical riethods (except the one using quaternions) and 
with the second graphical method to be described now. Once 
the secondary characteristic directions are fomd, the pri- 
mary characteristic directions and the characteristic retar 
dation can be found directly by using a different type of 
incident polarized light other than circularly polarized 
light. An example of finding the cliara-ct eristic parameters 
for a combination of two retarders using this chart will be 
presented at the end of this section. 

The Second Graphical Method; 

Squaring and adding Eqs . (5.^1 )and (5.^2), 

we obtain 

Sin^ 2eCos% + Cos^ 20 = Cos^ 2a) 

Cos^ 2esin2€ + Cos^C = Cos^ 2a) 

Cos^ 2a! = X and Cos^ 2© = y 

y ^ - - - Cot^ (5.59 ) 

Sin^0 

For a given value ofG , this equation represents 
a straight line. 

Since the set of values x = 1 , y = 1 satisfies 
Eq. (5*59) irrespective of the value of€ , Eq. ( 5 * 59 ) repre 
sents a family of straight lines all passing through the 
point X =1, y = 1. Also, 

if 0, X =1; and if € ::ir/2, then X = y 


or 

Let 

then 
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if 7 = 0 

X = Cos^e^ = Cos^ 2w 
or £= 2u 

Squaring and adc31ng Eqs. (5»^2) and (5*^3), we obtain 

Sin^ 26 = Cos^ 2tD Sin^ 2 <=( + Sin^ 2o 
or j = X Cos^ 2n( (5 *60) 

If X = 1 , y = Cos^ 2=<= Cos^ 2© 

or c( = e 

The families of straight lines represented by Eqs. 
( 5 . 59 ) and (5.60) can be plotted as shovm in Fig, 52. 

Eq. ( 5 . 59 ) represents straight lines of constant 
phase difference AC, AG, AF, AB,.., . etc. Eq. (5.60) re- 
presents straight lines of constant azimuth BA, BE, BD, BC, 
etc. The ordinates are lines of constant ellipticity 
and the abscissae are lines of constant amplitude ratio. 

Geometrical Construction of x = Cos 2w And 

y = Gos^ 2 6 Scales (Figs, 53 And 5^) 

A o 

In' Fig. 53 . BC = AC = 1 and BCA = 90°. Draw a ^ : 
semi circle on BC with 0, the mid-point Of ^BG as the cehter* 
To construct the constant phase difference line 
t'ahe)''Pbint' ;D);on).the ‘heBhcihClev d^D = 2e^ . Draw 

the perpendicular to BC* Dj^A. is the required constant 

''The proof for thisf is, 
as follows: ' '5' ^ 
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In the graph point B is the origin and line BC is 
along the x-axis. Also x = Cos^ 2a3 and j = Cos^e . 

On x-axis, 

2t! 

as sho’-ni above 

Now in Fig. 53, 

ED^ = BO + OD^ 

= i + i Cos 20^ 

= Cos^^ 

If the ellipticity at the point D-| is tan <d, the 
phase difference € is given by the relation 
C= Zm 

Also, the value of x at is given by 

X = Cos^ 

= Cos2ei 
or2u = 



We know that all constant phase difference lines 
pass through, the point A (1,1). Therefore the line AD^ 
is the constant phase difference line \>^ith e=e^» 

The procedure for cc^istructicn of constant azimuth 
lines is exactly similar. ■ In Fig. 5^. IlPC is the circle 
on AC with Oj. as center (A0-| = O^C ). To eon struct the 

constant azimuth line <=( =c( ^ , locate the point E on this 

A 

circle such that A 0^ B = • Draw the perpendicular EEl^ 

on AC from the point E to meet AC in E^ . 
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For point S , 


X = 1 


Gos"^ 2 0= C 0^ -f- 0-jE 
Cos 


Cos^ 


. « 9 = «=( ^ 

Also ^=B when x = 1 (shoivTi above) 

» . c( = o( ^ 

So, the azimuth at the point K| isc<^. 

We know that all the constant azimuth lines pass 
through the origin B. So ^ is the constant azimuth line 
with ^ . 

To increase the phase difference from to 0^ 
we move alaig the constant amplitude ratio (i.e. larallel 
to X-axis) and for rotation of the reference axes from azi- 
muth to 2 move along the- constant elliptic ity line 
(i.e. parallel to If-axis). 

The same chart can be used for other ranges also 
i.e. beyond ir /2 but the proper signs are to be taken into 
account. This is demonstrated at the end of this sectian. 

The. 'Striking advantages of this method are that ■ 
(i) it is very simple and easy to draw since it involves 
only straight lines and (ii) the nonlinear scales can be 
easily constructed by using the method described above. 
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The analysis for plane and circular polarise ope us- 
ing this c.iart is exactly similar as that doscrjlbed for the 
first graphical method. 

Tardy slethod ; 

In Fig. 55- the notations are same as those used 
in Fig. 50 . Point 0 is plane polarized light referred to 
the axes of polarizer. Point 1 represents light before the 
first quart or -wave plate referred to its o\m axes. Addition 
of a retardation of ir/2 is represented by movement along the 
constant aniplitude ratio line (parallel to X-axis) from 0=0 
line (A C ) to £= ir/P line (A B), that is to the point 2, 
which represents circularly polarized light. The sanB vib- 
ration, but referred to the principal-stress (inclined at 
V5° to the first quarter -wave plate axes) is represented by 
the point 3 (which coincides with point 2). This rotation 
of the reference axes through ^5° is represented on the 
chart by movement along the constant ellipticity line (i.e. 
Parallel to Y-axis) from «< = line (B C ) to<=c = 90 ° line 
(A B ). Addition of a retardation of £* in the model is 
represented by movement from point 3 to point 4- along crxis- 
tant amplitude ratio line till the representative point is 
on the censt ant relative phase line with £= 'f/ZTr+ G . 

Light vibration after the model but referred to the axes of 
the' second quarter-wave plate is obtained by moving the . re- 
presentativG point from 4 to 5 along the constant ellipticity 
line. SiJDtraet icai of a retardation of v /2 in the quart er- 
-wave plate is represented by movement along the constant 
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amplitude ratio line from point 5 to point 6 . In order to 
cut off light, the amplitude ratio along the analyser axis 
must be ^ . This can be achieved by mPhcing point 6 come to 
point 7 along the constant ellipticity line. The orienta- 
tion of the analyser axis is then whose value 

can be fo-und as follows: 

From Fig. 55 . 

x^ - = X6 - X5 = y^ - y^ 

But x^ = y^ = 1 

\ = ^6 

i*B* C 03 ^ ~ C os^ ^^6 ~ ^6 

Also €q - 

( • . when y = 0 , 0 = w ) 

And 00 = = 0 


Also 


But 


ei, = 

^7 "“^6 = ® 7."®6 

( 

= 0 

o( _ 0 

6 6 

■'ll. = 2 ( 01 ^ = 2e^*2=c^ 


when X = 1 , © =oc ) 


If the retardaticsi added in the model is P then p = 
£ 4 “ IT /2, If the angle between the analyser axis and one 
of the principal-stress axes is p , then 




09 f. 



FIG. 57. FIRST GRAPHICAL METHOD ( 0<o« 27r, o :<■ €< 2ir ). 
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= v/k- - P 
e =e ^ - -rr/a 
= 2=<6 - 7 r /2 

= - 2 P 

Derivation of J Circle Method Mithoijt Referring To The 
Poincare' Sphere Representation of Polarized Forms 

Consider the polar coordinate system (r, 5), Fig. 56. 
Let ^ = J = (C 03^20 + Sin^a© Cos^C)"^ 

then ^ = Cos 2a> (5.62) 

- using Eqs. (5-^'l) and (5*^2) 

Therefore, the magnitude of the radius vector r is 
a measure of the ellipticity. So, the curves of constant 
elliptic it y are concentric circles with center at 0, the 
origin. 

Let = 2=< 

where °C is the azimuth 

then r Cos°^ = J Cos 2c( (5*63) 

But from Eqs. (5*^'1) and (5*^3) 

Cos = Cos 2^J (5.6^) 

s o that . 

rcos^ =• Cos 20 (5.65) 

i.e. the x component (= rCos5) depends on e only. So, the 
loci of points of constant amplitude ratio are ,str;aight 
lines perpendicular to the X-axis. 
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Mow, 


X = r Cos = Cos 2© ■ 
y = r Sin OC = Cos 2a) Sin 2t< 


• . • 


x2 + 

Cos2£ 


Cos2 2 

1 


Cos^2a> Sin^2o^ 

C os2£ 

- -using Sq. (5-^2) 


( 5 * 66 ) 


Eq, (5.66) represents a family of ellipses with 
their centers coinciding with the origin and with the 
lengths of semi -minor and semi -major axes as Cos6 and imity 
respecti-vely. So the loc-us of points of constant phase 
difference is an ellipse belonging to this family of 
ellipses. The Wulff’s grid (stereographic projection) 
method can also be derived in a similar way mthout refer- 
ring to Poincare' sphere by choosing a proper expression 
for r j the radius vector. 

The graphical methods described can be used for 
values of phase difference and azimuth beyond the range 
0 - /2, The first graphical method is plotted in Eig. 57, 

for tho ranges t«C<2r, 0<£<2tr 

Some Examples of Using the Graphical Method , " 

We shall now take up the two examples of Schwieger 
for finding, the characteristic parameters. 
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Example (a) 

Consider the conibination of two mcdels (linear 
retarders). Model (1) produces a retardation of 
hodex (2) produces a retardation of = 60 dog. The angle 
l.-i'^/'oe.i -^ne axes of these two models is i^iven hT 


^^2 - 

ihe characteristic parameters of this system are 
foimd by one of the methods described earlier, i.e. the 
secondary characteristic direction and the characteristic 
retardation 5 are found from the exit light ellipse parame- 
ters when the incident light ellipse is circularly polarized. 
The primary characteristic direction is found by deter- 
mining the azimuth of light at the entrance which produces 
plane polarized light along one of the secondary principal- 
-stress axes at the exit of the system. This is found by 
tracing back from the exit point to entrance point of the 
system. This is demonstrated in Figs. 58(a) and 58(b). 

From these graphical methods, 

S = 69 •k- deg. 

'^'-1 (rotatory power) = 13 .S deg. 

Orientation of primary ) 

characteristic directicns with ) 

respect to the axes of the ) " deg. 

) ■ . ■ ■ 

first model ) 
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Results from the second gralSiical method; 

(The procedure is demonstrated in Pig. 59* > 

= - 10«8 deg. 

= 26.7 deg. 

5 =84- deg. 

iheoretical Values using the analytical solution 
(Chapter 2) ; 

= - 9.8 deg. 

P-l = 27.4- deg, 

5 = 78.6 deg. 

In this section ’re have developed two new graphical 
methods to trace the passage of polarized light through a 
series of optical elements. These two methods are compared 
with a few of the existing graihical methods. Procedures 
for the use of the new graphical methods are described with 
application to various polarise ope arrangements. The first 
method traces the transformation of polarized light through 
the same constant ellipticity and constant amplitude ratio 
curves. The second method also makes use of the same 
principles, but involves only straight lines in the graih 
which are convenient from the c cn st rue tioi point of view. 
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B. THEORY OF COMPENSAHON 

In this section, we shall first present a brief 
review of the existing methods of compensation and the methods 
of meas-oring the integral fringe orders. Then, we shall des- 
cribe, the development of a compensator which has certain 
advantages over the previous ones, and a general method of 
measuring the i.ntegral fringe orders. The theory for these 
methods malces use of the concept of optically equivalent sys- 
tem and various methods of determining this system for a 
series of retarders which were described in Chapter 3* 

The compensator can be defined as an optical device 
which introduces a required amount of retardation along a 
particular set of axes. This is equivalent to a two-dimensio- 
nal photoelastic model in which the principal-stress differe- 

45 

nee can be varied at will. Jerrard has given an account of 
different types of optical compensators in vogue. 

One of the commonly used compensators in photo- 
elastic investigatiens is the Babinet-Soleil compensator. 

This consists of two quartz wedges of same inclination with 
their fast axes parallel and a rectangular block of quartz 
in series with these wedges with its fast axis perpendicular 
to that of the wedges. The retardation is varied by the 
lateral movemait of the wedges and thus their combined thick- 
ness in the optical path. Other types of compensators inclu- 
20 

de Coker compensator which is a tensile strip (preferably 
made of the same material as that of the model) in which the . 
tension can be varied by a loading frame, a model beam 
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compensator, (Jellyman and ), in which the fringes 

resulting from lending are frozen into the material; mica 
plate compensator, (Phenister in which the angle of in- 

cidence of light is varied to affect a variation of the re- 
uardaticn introduced; etc. The compensation methods which 
donot involve the use of any additional optical elements are 
the Tardy compensation, De Senarmont compensation and the 
inverse Tardy method (Chapter 2). The underlying principle 
in aH these methods of compensation is to Introduce a reta- 
rdation equal and opposite to that introduced along the prin- 
cipal-stress direction (or the characteristic direction in 
a three-dimensional case), so that total extinction of light 
is obtained when observed in the polarise ope. The retarda- 
tion introduced can be related to the (secondary) principal- 
-stress difference in the model. 

Some of the desirable characteristics of a compen- 
sator can be listed as follows: 

(1 ) the principle of operation should be simple and the re- 
lation between the movement (angular or linear) in the com- 
pensator and the retardation that it introduces should be 
linear, (2) it should be inexpensive, (3) be accurate, (4) 
be useful for all wavelengths of light in the visible region 
(5) have a large field so that it can be used directly in 
the diffused light polariscope and (6) its retardation should 
not vary with time. 

The Babinet-Soleil compensator involves a very 
careful preparation of the wedges and so is expensive. 
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But, it can be used ad-Gli white light to find the integral 
fringe order upto about six fringes. But it does not give 
correct results when the variation of the material fringe ■ 
constant of the model with ’.-avelength is not linear (assu- 
ming that it is so for the quartz wedges and the rectangular 
quartz block in the compensator). Its field of view is very- 
small and is about 1 cm. in .iiameter. Coker's compensator 

If 7 

is cumbersome because of the loading frame. Phemister's ' 
mica plate compensator has a nonlinear relationship between 
the angle of light incidence and the retardation introduced 
in the model. It cannot be used to find the integral fringe 
order directly. The compensators which involve freezing of 
stresses into a piece of photoelastic material usually exhi- 
bit time dependence of the retardations. The Tardy, De 
! 

Senarmont and inverse Tardy methods can not be used with 
white light to find the total fringe order directly. They 
involve some errors if the wavelength of light is changed 
since the quarter-wave plates that are used shov; some depen- 
dence of their retardation ’;ith the ’-/avelength of light. 

It is the objective of this section to develop methods of 
compensation which involve minimum of the above drawbacks 
using the theory of characteristic parameters discussed in 
Chapter 3* 

28 

Pancharatnam , in his analysis on achromatic cir- 
cular polarizer and achromatic quarter-wave plate suggested 
that two quarter -wave plates with a half-wave plate in bet- 
ween will act as a compensator. Tuzi^9 suggested that a 
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combination of tvo retarder plates of equal retardance can 
be used as a quarter-v;ave plate. But this canbination will 
have a rotatory power. Childs U3t«i the coEbination of 
3 retarder plates to act as a quarter-wave plate with no 
rotatorjT power. The present coinpensator to be described is 
essentially based on those principles. The analysis for 
this will be made using the Stoke's vector. 

epabinaticin of Two Linear Retarders of Boual R etard ance p 
C onsider two linear retarders each of retardance p 
(with the angle between their axes being ©]. The transfor- 
mation of the Stokes vector through this system is given by.* 
^t Vq and 'j' be the Stokes vectors of light at 
the entrance and exit of the first retarder, with respect to 
its axes), ■ ' , / 









S-ii 


3li 


; S^ ! 

^0 - 

Sgi 


S 2 C os P - S^'^ Sin p 


^2 




^ C os p + S 2 |j^ Sin p 


_^3_ 


Let V.| and V 2 be the Stokes vectors at the exit 
of the first and second retarders referred to the axes of 
the second retarder. We then have. 


Y 


Ir- 


S.J Cos 20 + 82 


S.|Cos 2e + S 2 Sin 2© 

S 2 Cos 20 - S., Sin 20 


(S 2 Cos 2e ■ 


- 

- Sq Sin BqX^os P — S^Sin p 



C os p+ Sin p 

= Yo 


(S^ Cos 20 - S^ Sin 2e) 
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(3^ C'js 26+3 


+ S2 Sin 26) Cos 29 


>in £2© 


00s p 


(S2 Cos 26 - 3i Sin 26) - S3 Sin p) 
Es 2 Cos 2^ - Sin 2e ) Cos p - 3. 5 


- o:^ Sin p iCos 2© 


+ Sin 2e (s^ Cos 20 + S2 Sin 29) 

S3 Cos p + Sin p (S2 Cos 20 - S-] Sin ^) 


S2 = V 


(5.67) 


where 7 is the Stokes vector of light at the exit of the 

second retarder with respect to the axes of the first 
retarder. 

If Vq and are the Stokes vectors of incident light 

referred to the axes of the first retarder (reference axes) 
and the retarder of the optically equivalent system, x;e have 
Si S-j C os 2 + S2 Sin ap 1 ■ 

^0 = ^2 S2 Cos 2P-, - S^ Sin 2P-, ' = 


S^ Cos 2 + S 2 Sin 2B^ 

(S^ Cos 2P^ - S^ Sin 2f^ ) Cos 5 - S^ Sind 
I S3 Cos® + Sind (S2 Cos 2 p^ - S-, Sin 2 p^ ) 
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wh-are 

the retarder 
its axes. 


is the 
'-■f the 


Stakes vector of light at the exit of 
optically eqxiivalent system referred to 


sj^'s'': 


r 


sreci. 


CP 


t Tx 


system is equivalent to an optically equivalent 
single retarder (retardations, aximuthp^ with 
a:ce£ of the first plate) and a rotator vrith 


a r 00 at cry i:.-''.'-sr of . The 
vector through the optically 


transformation of the Stokes 
equivalent system is given by 


Also, 



Cos 2^ (S^ Cos 2p^., + Sin 2^'^ ) - Sin 2i> ^ Cos d 
(S^ Cos Sin 2Bi) - S^ Sin 3 j 

Cos 24- [ (S^ Cos 2$^ - Sin 2e^ ) Cos 5 - Sin3 ] 

+ Sin 24) (S^ Cos 2p^ + Sg Sin 26^ ) 

S'^ Cos 5 + Sin 6 (S^ Cos 26^ - Sin 2$-^ ) 



( 5 . 68 ) 



where ip = B-] + Cy and 7 is the Stokes vector of the light 

at the exit of the optically equivalent system (after the 
pure rotator) with respect 'to the reference axes. 

We shall now find the parameters p-j, 5 and 4 --^ of 
the optically equivalent system in terms of p and 6 using 
Eqs. (5*67) and (5.68). 
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iiet the incident light he circularly polarized for 
the original system and the equivalent system. The Stokes 
vector of this incident light vill be (0,0,1 ) .Equating the 
St oKes vectors at the exit of both the systems and substi- 
tuting for the incident Stokes vector, the vector (0,0,1 ), 
we obtain 

S^ = Cos^P - Sin^pCos 2 0 = C os 5 (5*69 ) 

S^ = - Sin pSin 20 Cos 28 + Sin o Cos p Cos 2 8 Sin 20 

+ C os p Sin p Sin 2 0 


= Sinews in 2i> 


(5.70) 


How, let the incident light be plene polarized along 
one of the reference axes (i.e. the axes of the first retar- 
der), sucli that Yq is (1,0,0) for both the original and 
optically equivalent systems. Then, by equating the Stokes 
parameters of light at the exit of the two systems, we obtain 


S^ = - Sin 2© SinP = - Sin 2p^ Sin 3 (5-71 ) 


If (5iSTr/2 and the value of p is given, the value 
of© can be calculated, i.e., the value of e f or \^ich the 
system acts as a quarter -wave plate. Fr'Om Eq. (5.69), for 
5 =ir/2, we have 

tan^p = Sec 20 (5.72) 

Since Sec 2 > 1 > Tr/4 

Tan^ p 1.1 or P 1. (5«73) 


Also, frem Eq. (5*71 ), the value of 
Sin 2© Sin p = Sin 


can be foundas 
( 5 .? 4 ) 



J’roffi Sq. (5«70) 


Sin 2 P= Sin p Sin 2 e ( - Cos 2 e + Cos P C os 2e + Cos p) 

(5.75) 

from which the values of -p and C:m be calculated. 

G^M nation of Two Quart or -wave Plat.^g 

• ie shall now investigate the behaviour of a combi- 
nation jf two quarter-wave plates and its use as a compensa- 
tor. Substituting p =ir /2 in Sqs. ( 5 . 69 ) and (5*71), we 
obtain 

- Cos 2 0 = Cos5 

or e =Tr - 20 (5.76) 

i.e. the characteristic retardation is linearly related to 
0 , the angle bet'ween the axes of the quarter-wave plates. 

From Sq. (5.71 ) 

2 p ^ = p = ir/a 

or = TT/h ( 5 - 77 ) 

From Sq. (5*70) 

- Cos 20 = Sin 2^ 
or Ti) = "e - ^ /k- 

= + "^1 
= ir/k + 

or ^-1 =0 - Tr/2 ( 5 - 78 ) 

i.e. the rotatory power is also related to© , the angle bet- 
ween the axes of the tv/o quarter -wave plates. 

Hence, two quarter -wave plates with angle e between : 
their axes are equivalent to a retarder (aximuth p = tt / if 
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with raspoct- to the axes of the first ^iiarter-wave phate, 
retardation 5 = tr - 2 0 ) and a pure rotator of rotatory power 
4-'^ = (0 - 'ir/2). 

So this combination can be used as a compensator 
if this rotatory power is nullified. This can be done in the 
following way. Consider the system of a light source (unpo- 
larized), a polarizer and two quarter-wave plates, Fig. 60. 

In Fig, 60(a) - 

OP - Axis of the polarizer 

OF^ - Fast axis of the first quarter -wave plate CR - 
one of the axes of the retarder of the optically equivalent 
system (one of the primary characteristic directions) 

OF 2 - Fast axis of the second quarter -wave plate 
OM - One of the secondary characteristic directions 
Also, 

— IT — ■jr — IT ir 

T= 

i.e. the angle between the secondary characteristic direc- 
tions and the axes of the second quarter-wave plate isir/f 
which is a constant. So, with this combination of polarizer 
and quarter -wave plate , we can produce the required ampli- 
tude ratio and phase difference along. the secondary charact- 
eristic directions (say OM and ON). Along these directions, 
the amplitude ratio will be tan 6 and the phase difference 
will be ( TT - 2 e ). If the amplitude ratio is to remain cons-' 
t ant, but the phase difference is to vary, we rotate the po- 
larizer and the first quarter -wave plate together by the 



V F-, 

'■ C-'t 

\ Cu / 



FIG. 61(a) 


Orientations of the axes of the four 
retarders 

FiiF 2»F3,F4 are the fast axes of the 
firstjsecond, third and fourth retarders 
respectively. 

R-),R2 are the primary characteristic 
directions of the first and second systems' 
Mi,M 2 are the corresponding secondary 
characteristic directions- 




FIG'SKb) 


FIG. 61 (c) 
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requirod amount. If the- phase difference is to be increased 
by p wo rotate the polarizer first quarter-wave plate combi- 
nation through an angle -p /2, keeping the second quart er- 
-vrave plate fixed. 

Similarly, the combination of two quarter-wave 
pla.tes and analyser can be used as a compensator, Fig. 60(b). 
In this case, the second quarter -wave plate and analyser are 
to be rotated together so that retardation is varied along 
the fixed direction OR (the primary characteristic directions 
remain fixed in space if the first quart er -v/ave plate .axes 
are fixed). The retardation added along the primary charac- 
teristic directions will again be (ir - 2e). By altering the 
value ofe , v/e can obtain the required amount of retardation. 
Rotation of the analyser is necessary to compensate for the 
rotation of the secondary characteristic directions with the 


second quarter -wave plate. A quarter -wave plate will have 
the same effect -vjhen it is rotated through -r radians. Its 
effect will ccaitinuously vary as it is rotated through 0 - ir 
radians. So, the angle e can take values from 0 to ^ i.e. 

We can introduce retardations frcaa tr to - ir (Eq. (5-76)) 
radians by these coffipensation techniques. 

Nw, , c on sider. ■ two' pair s-' of ' 'retarders,' each "p-ir . 
consisting of two retarders of equal retardation i» . The 
axes of these retarders and the primary and secondary chara- 
cteristic directions of the corrbspoftding optically equiva- 
lent systems are shown in Figs. 61(a) -and 61(b). 
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’ ^2 ’ ^3 ’ first, second, 

third ina fourth retarders respectively. 

, ^2 the primary characteristic directions of the 
1 ir st and seC'-^nd systems 

are the corresponding secondary characteristic 
direct ions . 

How examine Eqs. (5-72), (5.73) and (o-Tf). In 
these equations if is replaced by -0^, we obtain 
for P -j , -4* for ’i' and - “I'-j for i|-‘^ (*.* = 4' - P ^ from Eq, 

(5.68). i.e. if the angle' between F-^ and ^2 is©^ and the 
angle boto/een F^ and Fj^ is - © , the rotatory powers of the 
equivalent systems for the pairs F^ , F^ and F^, F^ will be 
equal but in opposite directions. Also, the optically equi- 
valent system can be either a retarder and a rotator or a 
rotator and a retarder. In both cases, the rotatory powers 
and characteristic directions will be equal. But the retar- 
der axes will be along the primary characteristic directions 
in the first case and along the secondary characteristic 
directions in the second case (Chapter 3). Let the axes of 
the four retarders be oriented as shown in Fig. 61(c), 

Lot the combinatioi of the first two retarders be 
equivalent to a retarder (azimuth with respect to the 
axes of the first retarder of retardaticsi V2) and a pure 
rotator of rotatory :power ■. . d?hon the confcinatiai of the 

third and fourth retarders will be equivalent to a system 
with a retarder (azimuth - with respect to the axes of the 
third retarder of retardation v/2) and pure rotator of 
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rotatory power - i|;^,or a system with a pure rotator of 
rotatory power - ^ and a retarder (azimuth -"pi - with 
respect to the axes of the third retarder of retardation 
'^/2). So, when the four retarders are arranged in sequence 
as shO'/n in Fig. 61 . the rotatory power of the first two 
retarders will he cancelled hy the rotatory po-rer of the 
combination of the third and fourth retarders. And the com- 
bination of the four retarders acts lire two quarter-wave 
plates. The orientation of the axes of these quart or -wave 
plates will be alcsig the primary characteristic directions 
for the first two retarders and the secondary characteristic 
directions of the combination of third and fourth retarders 
respectively. The orientations of these directions are 
found from Eq. and these are fixed relative to the 

axes of the retarders once the values of 6 and p are fixed 
(pis the retardation of each of the retarders). Thus the 
combination of four retarders each of retardation (F > 
from Eq. (5 •73)) can be used as a combination of two quart er- 
-wava plates. This in turn can be used as a compensator as 
described abcve- By varying the value of e ^ , the angle 
betwTeen the retarders, this combination can be used as a 
compensator for any wavelength of light. The range of wave- 
length' is however limit ed by the value of F and the dispers- 
i'ln characteristics, {i.e,^ the variation of the retardation 
with wavelength) of .the ■ material ' of the retarders. : 

In the above discussico, we have shorn that the 
ccmibination of two qusirter-wave plates can be used as a 
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compensator provided the rotatory power introduced by the - 
system is nullified by rotating the polarizer or analyser 
as the case may be. To develop a compensator which does not 
introduce any rotatory power i.e. to make it act as a single 
retarder, let us consider the case of a half-wave plate placed 
between two quarter -wave plates- 

A Half-Wave Plate Betv/e e n Tw o i3uarter-Wave Plates 


It will now be shovm that a system consisting of a 
quarter-wave plate, a half-wave plate and a quarter-wave 
plate in that order, can be used as a compensator. The axes 
of the two quarter-wave, plates are parallel. The angle bet- 
ween these axes and the axes of the half-wave plate is e • 
The transformations of Stokes vector through such a system 
are given below 


Let and be the Stokes vectors of light at the 

entrance and exit of the first quarter -wave plate referred 
to its axes. Then, we have, 





Si 

^0 ^ 







Sa 


If and are the Stokes vectors of light at the 

entrance and exit of ihe half -wave plate referred to its 


axes 
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Cos 20 - Sin 20 

1 3 


S Cos 20 - S_ Sin 2e 
' 3 


-S 3 Cos 20 - S-j Sin 2 8 S 3 Cos 28 + Si Sin 2 0 


Also, 

Cos 2© (S^ Cos 20 - S^ Sin 20 ) - Sin 20 (S^ Cos 2© 

+ S-j Sin 2© 

% = Cos 20 (S^ Cos 20 + S^ Sin 2 ©) + Sin 2 © <S^ Cos 2 © 

- S 3 Sin 2©) 


Cos 2 © (S^ Cos 26-8, Sin 26) - Sin 2© (S, Cos 2© 
1 3 3 

+ S^ Sin 2© ) 


Cos 2 © (S^ Cos 2 © + S^ Sin ^) + Sin 2i(S^ Cos 2 © 


- S^ Sin 2 ©) 



( 5 « 79 ) 
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inhere, V^; Stokss vector of light Incident on the second 
quarter -wave plate with respect to its axes. 

Y^; Stokes vector of light emerging from the second quart er- 
-wave plate with respect to its axes. ' 

The Stokes vector at the exit of the optically equivalent 
system is given by i^q. (5.68). Comparing this equation with 
(5.79) for = S 2 = Oj = 1 , we obtain 

S2=Cos5=Co3>+eor5 =lfe (5,86) 

= Sin 5 Sin 2 ^ = - Sin 4-e 


or Sin 2 ^ = - 1 


6 = ^ e 



II 

• 

CD 

• 

•H 


( 5 . 81 ) 

for s 

t 

= 1 , s = s 

2 3 

= 0 


= - 

Sin 2p^ Sing 

= Sin 4-0 

( . . s = 4 ©) 

• 



0 

II 

1 


(5.82) 


* • The r ot at ory power of the ) 

) = t. = f - P. = 0, 

combination ) * * 

bo, this system acts as a single retarder with retardaticn 

4-0 and azimuth - with respect to the axes of the first 

quarter-wave plate. The retardation introduced caa be 

varied by varying the angle 6 . Thus, this system can be 

used as a compensator. 

The \diole system, i.e, a quarter-wave plate, a 
half-wave plate and a quarter-wave plate, can be constructed 
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With eight retarders each of retardation p, tan p>l 

(for t.ivi Part j.cular i»ravelength of light wider ccnsideration )« 
This can be done as follows. 

Consider two types of quart er-wave plates (a) 
and (b) each of which is a combination of two retarders. 

Tae orientations of the axes of these retarders are shown 
in Fig. 62. In the (a)-type combination, the rotatory pcwer 
13 clockwise where as in the (b) type it is anticlockwise. 

The dilference between these two types of combinations is 
in that the angles between the axes of the two retarders in 
types- (a) and -(b) is e and - erespectively. If the type-(b) 
combinatico is placed after the type-(a) combination, the 
rotatory powers of the combinations cancel each other so 
that the two combinations behave as two quarter-wave plates. 
So, to construct a compensator we can use four such combi- 
nations (combination of two retarders each of retardation p ) 
with tne first and third combinations of the type (a) and the 
second and fourth combinations of the type (b). Then the 
rotatory powers of the first and third combinatit^s are 
nullified by the rotatory powers of the second and fourth 
combinations respectively. Thus, the whole system of these 
eight retarders acts as a series of four quarter -wave plates 
which in turn can be used as a compensator as described 
abdve. By varying the angle e, the angle between the axes 
of the two retarders in a two retarder combinaticxi, accord- 
ing to Iq, (5.72) i.e. tan^ p = Sec 2 e , we can use the com- 
pensator for any wavelength of light for which the condition 
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tan P >1 is satisfied. (This, in practice is not a severe 
restriction since the retardation p does not vary appreciably 
within the visible region of the electromagnetic spectrum 
for most of the iSiotoelastic materials). Using this method 
large field compensators can be made using large sheet type 
retarders (ex* frozen-stress sheets). The angle 6 can be 
found experimentally for different wavelengths of light with- 
out any calculatiais. It is advisable to cut all the eight 
retarders frbm a single large retarder of uniform retarda- 
tion. The use of a large field compensator in a lens polari- 
se ope eliminates refocussing of the lenses when the- compensa- 
tor is introduced into or removed out of the field of view 
of polariscope, which is necessary with a Babinet-Soleil 
compensator because of its small field. The compensator 
described in this secticxi thus has some additional advantages 
like large field, cheapness etc. However,, the integral fringe 
order cannot be determined directly with this comp^sator 
as can be done with Babin et-Boleil compensator (upto a 
certain nraaber of fringes, usually about 6). However, the 
integral fringe orders upto vlrtmlly any nunier can be 
determined by an altc^ether diiferent method. He shall now 
develop this method in the following sectiem. 
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C • DgriRMDiATION OF IICTB^HAL FRINGE OHD^S III PHOTOELASTICIT Y 

In analysing a Eiiotoelastic snocimen having no 
stress-free corner or, in which the formation of fringes 
during the leading process cannot he observed, it is desi- 
rable to have some means of determining the Integral fringe 
order at a point when the fractional fringe order at that 
point has been determined by the Tardy or other method. A 
Babinet-Soleil or a similar compensator can be used to dete- 
rmine the integral fringe order upto a certain value. How- 
ever, there is a serious .point which needs attention when a 
Babinet-Soleil compensator is used. In practice, to deter- 
mine the integral fringe order at a point by the help of a 
Babinet-Soleil compensator, a white light source is used and 
retardations are added or subtracted tectil a zero-order 
fringe (a black fringe) appears cn the screens This process 
assiiaes that the compensator acts as a continiatioo of the 
model and its behaviour is similar to that of the model. 

This implies that the model and the compensator introdxjce 
retardations which are equal in magnitude but opposite in 
sign for all colours cemtained in white light. This is true 
only if the material of the model and that of the cempenaator 
exhibit same tjrpe of dispersicn characteristics (variatiOT 
of retardaticn with wavelength of light). This is not the 
case in general. Most of the compensators are supplied 
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v;ith a calibrati-cn chart which sho’^s the dependence of the 
compensator on the wavelength of light. Fiirther, experience 
has shovm that using a Babinet-Soleil compensator and a 
photometer, instances occur where the reading on the photo- 
meter is almost zero for two consecutive fringe orders, mak- 
ing it difficult to identify which one corresponds to the 
zero fringe order. This experience is quite noticeable with 
a urethane rubber model. 

In 1963, Pant ^ showed that by determining the 
fractional fringe orders ■''ith two •''avelengths and 
one can determine the fractional fringe order obtained by 
using a third wavelength A , His analysis ignores the de- 
pendence of the stress -optic coefficient A and also is res- 
tricted to the case where the difference between the fringe 
orders associated with A^ and .A^ is cxily a fraction, i.e., 
the integral fringe orders associated with A^ and ^2 
equal. 

We shall now shew, how by using two wavelengths, 
one can determine the integral values (and their multiples) 
for both wavelengths, associated with the fractional fringe 
orders observed at any point. The solution can be obtained 
by using a nomegram or constructing a ready redeemer for any 
two wavelengths A^^^and ^2 shown in this analysis. Fur- 
ther, the usual assumption that the stress -optic coeffioient 
is independent of wavelength, is not used. Its dependence 
on A is explicitly shown by experimental results. 



Stress -Op tic Law 
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■ihe totress-optic law can be expressed a 

(P - q) = ^ (p - q) 

wnere p j retardation 


( 5 . 83 ) 


h 


(P - q) 


thic m-jss of tne model or length of light path 
inside the model 

principal or secondary (for frozen stress slices: 
principal-stress difference 


A J wavelength of light used 
C : Stress -optic coefficient 
f : Model fringe constant in tension 


In practice, the stress -optic coefficient C is 
considered to be Independent of wavelength. However more 
careful observations have shown 20,49 g is 

to some extent by the wavelength. For our present analysis, 
ve shall assume that 0 depends cn 1 and hence, f is a func- 
tion of 1 . If, u, any experiment, C is found to be inde- 
pendent of j , the analysis remains unchanged. 

Let and be the wavelengths of the light used 
p-l , Pg the corresponding retardatiens. In general, 
these retardations will consist of an integral value and a 
fractional value. Let + 8^ ) and ^ 8^)^ 

«liere , Og are positive Integers and , 8g are positive 

t . 


fractlcns less than unity. Without loss of generality, 
let f^ be greater than f^, so that their ratio z - fl 

’ fa 


>1 



From Sq. (5.83) 

(ni ( p - q)' 

(^2 + fla) = 4 ( P - q) ’ 

"^2 

^2 ^2 _ £l _ 

+ ai ^2 

or 

n2 - 2 =2 5 ^ -52 =€ (5,9+) 

If the stress-optic coefficient happens to be inde- 
pendent of X , then, 

2 = Z*= 

The experimental investigations will show that using 
2 instead of 2 will lead to erroneous results. In practice 
the fractional fringe orders 5 ^ and ^2 determined and 
hence £ becomes an experimentally determined factor. Eq. 
(5.8^), will- then have to be solved for n^ and n 2 subject 
to the conditicfi that these are whole numbers. 

General Observaticns : 

a) Since z>1 (as experimentally cSsserved), P 2 > P| . This 
is obvious when C is independent of X , 

Hence, 

^2 ^ ^2 ^ +3-1 

and n 2 > (5-85) 

b) Since 0 <- 5 ^ <1 and 0<5 a given z. 
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and 


e 

max 


) 

max 

^ min 

€ . 

min 

It 

O) 

.-A 

^mln 

^^2^max 

i*e. 

2 > e 

> - 1 



( 5 . 86 ) 


c) In Sq, (5.^), if we add N to r \2 and N/z to n^ where 
N and N/z are integers, the equation is still satisfied. 
Let N be the smallest integer such that N/z is also an 
integer. Then the solutions for ng and n^ are 


Cn 2 + SN) and (n^ + SN/z), where S = 0,1,2, •.... 

( 5 . 87 ) 

In other words, if we find two basic solutions ng 
and n-j , satisfying Eq. (5.8^), then the multiple soluticns 
are given by the above equation. 


Constructicn of Hcmc^ram ; 

The wavelengtb in visible spectrum TOries 7700°4 

(extreme red) to 3900 °A (extreme violet) and for this range 
z = < 2. In practice, the usual range is fraa red (about 

7500®A) to blue green (about 5000 °a) giving for z a mlue 
less than 1 •5» 

In order to solve Eq. (5.8^), wie can construct a 
nomogram as shown in Fig, 6!3* BLne AF represents z, line 
AD, the integral values of n^ and line BE, the integral 
values of n^ . The scales for the values of z, Ug and 
are the same. Frcra A, rays A^ , Ag , A^ »•••• etc are 
drawn passing through the integral values of n-| = 1 ,2,3, • . -etc* 
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liie slopes of these rays are 1,2,3, . -.etc., hy construction. 

Let us say that in a particular example, the ratio 
of the two nodel fringe constants corresponding to A-j , A2 > 
is O'! 7^2 = 2^. In Fig, 63. A G = • A horizontal line G H 

is .ii'a;*m cutting the rays AA^ , at g,II,G.., 

etc. as sho^m. As the slopes of these rays are 1 , 2 , 3 , . . . .etc 

G K = z^; G N =2 z^ , G 0 = 3 ,etc. Line G H cuts 

the vertical lines representing the integral values of n2 
at L, N, P, 0,...., etc. Distances L K, M fv, P N, T N, etc 
represent the values of C that one determines experimentally. 
The procedure for using the nomogram is as follows; 

For a given problem, let the measured values of fra- 
ctional fringe orders be ^ so that z 5-| - §2 = €1 • 

Identify a point which is the intersection between the line 
G H and one of the vertical lines representing the values of 
11.2 (like L, M, P, T,...) such that the distance between this 
point and the slant ray immediately to the left of it is 
equal to £-j» Let us say in our case, M is such a point 
where M K =C -j where K is a point on the ray iimediately to 
the left of M, i.e. ray A-j » The vertical line on which pcdLut 
M lies gives the value of n2 (in this case ng = 3 ) and the 
slant ray immediately to the left gives the value of n-| 
associated with it (here has = 1 ), The proof for this 
is evident frc® the fig^ure, since M G - K G = M K, , i.e., , 
n^ “ = €-} , Sq. where 02 = M G = 3 and z-j n-{ 

= K G = with n-j = 1 . 
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Generally for a given value of 0^, there will be 
more than one point (lying on one side of the vertical lines), 
such that the distance between this point and the slant ray 
immediately to its left is equal to e^. These give multi- 
ple values according to Eq. (5.87). This aspect will be dis- 
cussed later, » 

One observes from Fig, 6)8i. the values of € can cnly 
be discr'ite. This is evident by the nature of Eq. (5.^). 

For, we can have a problem where n 2 » and z are at our 
command, but n-, and «-{ have to satisfy Eq. (5.8f) and also 
that n^ should be an integer. 

In the above procedure, we have assumed that 0 is 

positive. If 0 turns out to be negative say,0 = - 02» then 

we measure the distance frcm the point sixjh as M to the slant 

ray immediately to the right ; i.e.; if M is the point such 

that M lif (where^ N lies ce the ray iaiaediately to the right) 

is equal to | - 02 1 , then n 2 = 3 and n-j =: 2, since the ray 

on which H lies is A 2 with associated n^ = 2, Fac negative 

values of 0 , the point can also lie on the segm^it 0 G* 

Negative values cannot be numerically greater than i as ^own 

by our previous analysis. < 

To make the above procedure clearer, let us take 

a specific example. Let z = 1.25- A horizontal line pr 

correspcnding to this value is drawn as shown in Fig. 63, 

cutting the vertical lines at a , a^ , a etc. Let 

• 2 3 

us assume that we find experimentally z - ^2 = 0 » 25 * 
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Now we look for a point lying on one of the vertical lines 
(such as a^ or a^ or a^ or..,) such that the distance bet- 
ween this point and the slant ray immediately to its left is 
0.25* In our case is such a point since ^ b^ = 0*25. 
Hence n 2 = ^ (point ^ lies on this) and n^ = 3 (since b^ 
lies on ray A 3 ). We can see that these values of n^ and 
sr.tisfy Fr. (5.Sf) with z = 1 .25 and J = O. 25 . Qn the other 
hand i_ c, happens to - 0 . 5 , we select a point such that 
the disuance between this point and t)i.e slant ray immedia- 
tely to the right is equal to 0*5- Hence, we get a point 
3-2 since a 2 b 2 = - O .5 and b 2 lies to the right of a 2 » 

Then ng = 2 and n^ = 2. 

A value of €= - 0.25 will give us points a^ and 
b-j with n 2 = 1 and n^ =1; a value of e = - 0.75 gives us 
points a^ and b 3 with n 2 = 3 and Uj = 3 ; a value of p = 0 

gives us point a^ (b^^ coincides with it) with n 2 = 5 and 
n^ = 

Miiltiple Values for n^ and n 2 : 

From Fig. 63. we notice that for e = - 0.25, we get 
points a-j and b^ as noted above with n 2 = n^ = 1 , and also 
points a^^ and b^ such that a^ b^ = j - 0 . 25 | with associated 
values n 2 = 6 and n^ = 5* These also satisfy Eq. (5*S+), 
since, 6 - x 5 = -0.25 • This is in agreement witfe 
(5-87) which gives all multiple values. For our case, 
the lowest integer N such that N/z is an integer is 5 . The 
solution for n 2 and n^ according to Eq* (5.8?) are: 
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+ S H = 1 -r 5 S = 1,6,11,16, ,etc. 

and n-i + -§-1" = 1 + 4 S = 1 ,5,9,13, • etc. 

2 

all satisfying Eq. (5 .84). This should not cause much of a 
problem in .practice, since the difference between cne set 
of values for ni,n2 and the subsequent set is so large that 
an analjrsis of the problem c-an give a clue as to the appro- 
priate value. Further, as will be shown later, one can 
always choose two appropriate wavelengths such that the 
difference between two multiple values is so large that 
there will be little doubt as to the correct value. If 
ambiguity still exists, case can use the francticnal fringe 
order ^ for a third wavelength of light That is, one 

set of multiple solutions are found using the values of 5^ 
and ^2 ^ second set is found using the values of ^2 

and The Value of n^ common to both these sets gives 

us the correct solution. For example let these two sets of 
multiple solutions be given by; 

(n^, n^) = (1,2), (6,6), (11,10), (16,14), (21 ,18),.. .etc. 
and 

(n^, n 3 ) = (0,1), (3,3), (6,5), (9,7), (12,9), ....etc. 
then the correct values of the integral fringe orders are 
given by; 

(n^ , n^, n^) = (6,6,5 )» 

All Possible Values for g And Associated n 2 , Values: 

In practice it is not necessary to use a nomogram. 
For a given ratio z of -the wavelengths, one can write down 
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all L-he possible values of 0 that one can determine experi- 
mentally and write down the corresponding values of n^ and 
Tuis procedure has several advantages, since one has 
a sort of a ready -reckoner using which, the associated 
values of n^ and n^ with any given 0 can be read out. Also, 
If duG to some experimental error, the value of 0 is not de- 
termined exactly, the ready reckoner can immediately tell 
what should be the correct value of 0 for any z and a given 
range of n^ , n^ . This is so because, for any z, the values 
of 0 can only be discrete as was pointed out earlier and as 


seen in Fig, 63 . and also for any range of n2 , 
between any two successive rays A. A or A A. 
etc.) the values of 0 are fixed. 


n-i ( lying 

■^ * * * * 


In Fig. 64 -., a portion of Pig. 63 ., is shotm i-dth 
G H, the horizontal line cutting the rays A-j , A2 , A^,... 
etc. at distance of z, 2 z, 3 z,..., etc. from G. Let 
be the integer next to z, i.e. if z = a-, + where is 
the inttiger part and p-j the fractional part less than 1, 
then + 1 . If 0 , . Let Y2 be the integer 

n<-At to 2z; i.e. if 22 = b,^ where is the fractional 

part less than 1 , thovi a2 ^ «2 / otherwise, 

Y2 =: a2. This is continued until Y^j, tAere Y^jj is the 
integer such that Yju = mz. The points corre spending to 
T-j , Yg, . . ., are sho*/m in Fig, 6§ . Note that the 
slant ray in Pig. 63. imssing through the vertical Une " 
n2 = m will pass through the point y . 

Ill ' ' 
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From Fig, 6 3. the values of € and the ass>ociated 
■values of n2 and can be vrittan dO’v.Ti. This is shown in 
Table 5 . 1 . The procedure for constructing this table is as 
follows; 

^ = n2 - n-| z 
and z > € >_ •] 

Case (i) c. is Positive: 

Write down the values of £ , n2 for = 1 , then 
for r:^ = 2 etc. For n^^r: 1, the first possible value for £ 
would /■-] - z with Ug = Y-j * To obtain the other values, 
the value of n2 should be increased in steps of 1 so that 
we obtain 1 + - z, 2 + - z, 3 + r-j - z, ....etc. for 

€ and 1 + , 2 + , 3 + Y^,....etc. for n 2 . This is 

continued until the value of e exceeds z. Let the maximum 
value of e that does not exceed the value of z with ni= 1 
bee . 

i .e. 0 < z . 

=; 02 n^/' z : ' 

= 02 - z 

. (% ., i-j are the c,orrespondin.g mlues of 
02 and n^ respectively so that f = ^2 - z) 
i.e. 02- z < Z'' 

02 <'2 z 

The maximum value of 03 that satisfies this inequ- 
ality is (Y2 - 1 )• Only the maximum value is to be taken 
since S is to be maximum. 


( 5 . 8 ^-) 

( 5 . 86 ) 



2^ 

• 

• • - 1 - Z 

~ ^ n-j = 1 

Sinil:irly, 

for n-j =2 
^ - 2a 

TI2 ~ 2z "i a 
or ^2 C 3z 

"2 = ^3 - 1 

for = 1 

^ — H 2 f ^ a 
or n 2 < (1 + 1 ) a 

" 2 = ’' 1+1 - 1 

Case (il) G is Negative 

^ = ^2 - n-j a, a < 0 < - 1 

n-j can not be aero (If it is zero should be be 
gative to make C negative) 

ff = 1 a positive integer 
= n^ - n^ a = n^ - i z > - 1 
0 > n2 - 1 2 > -1 
i-e. 112 <12 
and n 2 > 1 2 - 1 
i.e.la>n2>l2-1 
^ , n2 = Y2_ - 1 

The table for negative values of 0 is ec^structed 
using , this equatieti# 
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T.iELB 5.1 



Pos 

sible 

Values r 

of 6 

and 

Ass --x: dated 

n 2 and n. 



Pos 

itivo C : 








e 

ni 

^2 


e 

ni 

^2 


^'1 

- Z 

1 

■^1 


^2 ~ 2z 

2 

^2 

1 + 

^1 

- z 

1 1 

+ Y^ 

1 ■-+ 

- 2z 

2 

1 

2 + 


- z 

1 2 

+ ""l 

2 + 

^2 - 2a 

2 

2 4- -^2 

3 + 

■^1 

• 

- z 

1 3 

9 

3 + 

^2 " 

# 

2 

3 + 


• 

« 


• 

• 

» 

• 


« 

• 

• 

• 

• 

• 


• 


• 

• 


* 

* 

• 

1 

C\J 

- 1 

- z 

1 

- 1 

^3 " 

1 - 2z 

2 

'^3 






e 


^1 


>2 


continue upto 

Yuj^l-Cm - 1) z 


m - 1 



n^ = m 

1 + Ygj.-j -Cm - 1) 

z 

(m - 1) 

1 


where m z is 

2 + -£m - 1 ) 

z 

(a - 1 ) 

2 

* 

an integer 

3 + Yia-i~(m - 1) 

z 

(m - 1 ) 

3 

+ lk-1 

= Yr« 

‘IE 

* 






(m -1 ) - 1 


- Z 




WAVELENGTH Vs. ABSORPTION 
FILTER NO. 2938 

MAXIMUM TRANSMISSION AT 5950 A 



WAVELENGTH 


WAVELENGTH Vs. ABSORPTION 

FILTER NO. 3214 

MAXIMUM TRANSMISSION AT 5004 


o< 






lie gat lire £ : 


e 

ni 


«2 

1 

1 

1 

^1 

- 1 

Y.-. - 1 - 2z 

d 

2 

^2 

- 1 

^ 3-1 " 

3 

"3 

- 1 

Yj_^ - 1 - 4-3 

4 

^4 

- 1 

continue till 

• 



Y^ - 1 - ms 

• 

• 


• 

• 

m 

ffl 

Y 

- 1 

= - 1 

m 


. .Y = mz 

. m 


Exajaple And Experlaental Results: 

To show the dependence of the stress-optic coe- 
fficient CO the wa-velength, three filters with wavelengthr 
5950°A, 552 OOA and were nsed to nieasiire the fringe 

c onst.aiit . The" model, was a ■circ.ptlah .'disc , of ' diaiiete.r (d) ■ 
4.73 cm. with a thiclmess (h) of 0.6^- cm. Figs. 65(a), 

(h) and (c) show the characteristics of the three filters. 
Fig. 66. shows the plots of load in gms. vs. fringes at 
the center of the disc. The geometry of loading is also 
shown in the figure. Table 5.2 gives the values of the 
material fringe constant F and the stress-optic coefficient 
C for each of the filters. 

At the center of the disc (p - q) = 8w>4rt d 
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TABLE 5.2 


Wavelength 

Material fringe 
constant (F) 
gm./cm, fringe 

C 

°A fringe, cm/gm 

595 QOA 

210.1 

28.28 

55200a 

192.0 

28.79 

5oc4°a 

172.0 

29.11 


Though the values of C dcnot differ very much, 
still neglecting them can give erroneous results as will 
he sho-'j^n subsequently. For the determination of the inte- 
gral fringe orders frcmi the measured fractional fringe 
orders, two filters with 5950 °A and 5520 OA wavelengths were 
used. The fractional fringe orders were measured by the 
Tardy method. Table 5 *3 gives the values of the experi- 
mentally dbserved fringe orders (Integral plus fTacticoal 
values) and the theoretically calculated values based m 
the model fringe constant. This table reveals that the 
fractional fringe orders measured by the Tardy method are 
■quite ■ accurate* 



2^2 


5.3 


boad P(gni} at 
end of lever 

5970 '^A filter 

Sxptal. Theor. 

55 20*^ A filter 

Exptal. Theor. 

575 

5.25 

5.23 

5.73 

5.73 

625 

5.70 

5-68 

6.23 

6.23 

675 

6 . 1 if 

6 .14 

6.73 

6.73 

725 

6.61 

6*59 

7.23 

7-22 

775 

7.06 

7.04 

7.73 

7.72 

825 

7.50 

7.50 

S .20 

8.22 

In Fig . 

66 w 

_ % = 39 . 7P 



(art)) 11.2 


= 3.5^7P 

From Table 5*2, we have 

^ = 5950OA, =210.1 gm/cm. fringe 

^ = 5520° A, F 2 = 92. O gffl/cra, fringe 



Tables 5.^ and 
with § and z* , 


= ^ - 1.09^ and z * = ^ 1 ,078 

5:5 give the values of e that one can expect 

These are the ’*2‘ead3r -reckoner '* tables. 


i ab ie 5 5 vhich 


Independence of C cn » ) will give ei-roieous 


(assiiaing 
as shown 


in T ab le ^ ^S • 



^3 


TA3LB 
(Using z) 


in s 

0( Positive ) 

ni 

^2 

G(I‘!egative ) 

ni 

n2 

Z = 1 

0.906 

1 

2 

- O.C94 

1 

,S 

2z =2.1 33 

0.812 

2 

3 

- 0.188 

ni 

2 

3z =3 »k.34i 

CO 

57 

• 

CO 

3 

4 

- 0.282 

3 

3 

4-3=4 .376 

0.624 

4 

5 

- 0.376 

4 

4 

52=5 *4 70 

0.530 

5 

6 

- 0.4 70 

5 

5 

6z = 6,564 

0.436 

6 

7 

- 0.564 

6 

6 

72 = 7.658 

0.342 

7 

8 

- 0.658 

7 

7 

8z = 8.752 

0.248 

8 

9 

- 0.752 

8 

8 

92 = 9.846 

0.154 

9 

10 

- 0.8+6 

9 

9 

10z =10,940 

0.060 

10 

11 

- 0.940 

10 

10 


1 .060 

10 

12 



• m 


etc . 


etc 


etc 


2lfif 


TABLE ^.5 
(using z * ) 


3|C. 

m z 0 (Positive) 

z = 1 .078 0.922 

2z = 2.156 0.84+ 

3z = 3.234 0.766 

4z= 4.312 0.688 

5z = 5.390 0.610 

6z = 6.468 0.532 

7z = 7.546 0.454- 

8z = 8.624 0.376 

9z = 9.702 0.298 

10z =10.780 


”2 G*( Negative) n2 


1 

2 

- 0.073 

1 1 

1 ! 

2 

3 

- 0.156 

2 

2 

.3 

4 

- 0.234 

3 

3 

4 

5 

- 0.312 

4 

4 

5 

6 

- 0.390 

5 

5 

6 

7 

- 0.468 

6 

6 

7 

8 

- 0.546 

7 

7 

8 

9 

- 0.624 

8 

8 

9 

10 

- 0.702 

9 

9 

10 

11 

- 0.780 

10 

10 


0.220 


10 10 




2^5 


TilELB ^.6 


load F 
(gms. ) 

Fr ict . 
■^1 

...Fr. Ord. 

5 

2 

C — 2 ^^ ""®2 

e* 

-52 

Fxptal. 
ni n 2 

fr-omC fromC* 
ni n 2 n-j n 2 

575 

0.25 

0.73 

- 0.46 

- 0.46 

5 

< 

y 

5 

5 

6 6 

625 

0 . 70 

0.23 

+ 0.54 

- 0.53 

5 

6 

5 

6 

6 7 

675 

0.14 

0.73 

- 0,58 

- 0.61 

6 

6 

6 

6 

8 8 

725 

0.61 

0.23 

+ 0.44 

+ 0 ,4 3 

6 

7 

6 

7 

7 8 

775 

0.06 

0.73 

- 0 .66 

- 0 .66 

7 

7 

7 

7 

8 8 

825 

0.50 

0.22 

+ 0.33 

+ 0.33 

7 

8 

7 

8 

9 10 


For the wavelengths ^ and ^2, the fractional 
fringe orders < 5 ^ and 5 ^ are experimentally detenained. 
Using 2 and z* the values of e and e* are then calculated. 
Frcm Tables 5-4 and 5.5, for the values of e and e* , the 
corresponding values of n^ and 02 are noted. These values 
are compar0d;,:With .the ■experimentally .observed va.lues.. As 
pointed out, values based on i.e,, neglecting the depe- 
ndence of 'C: on-. 'X can: lead values, /It is more •' 

appropriate., to use (i.e. v^so, assuming that, the 

values in Table .are c.alculated .more ac.curately,, the ■ 
experimentally.,deteriilned.'.values of e my -.differ - slightly 

as, given;..in,lable,,' 5 '. 6 \:^ ,on the-method employed to 

determine the fractional fringe '...orders..' ■, 
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C-onclusions : 

Using two wavelengths and knowing the material 
or model fringe constants for these' two wavelengths, one 
can determine the integral fringe order existing at any 
point in a plane model by measuring the associated fractional 
fringe orders. It is shown, that the usual assumption of the 
stress-optic coefficient being independent of X may lead to 
erroneous results. For any two given wavelengths and a 
model, a ’ ready -r eckcxier ’ table can be formulated to deter- 
mine the integral fringe orders based on the measured frac- 
tional fringe orders. Alternatively , a nomogram can be 
constructed. 

In this chapter, we have discussed various methods 
of finding the integrated optical effect in a ptiotoelastic 
model under three -dimen si dial state of stress graphically. 

Two new grajtiical methods which are derived without referr- 
ing to Poincare sjSiere representation are presented, (kie of 
these methods has the advantage that caitalns dily straight 
lines and so it is easy to construct. In the secdid section 
we have presented the construction of a large field coHipen- 
sator using 8 retarders of equal retardance. This is esse- 
ntially an extension of the works of Pancharatnam and Tuzi . 
Also presented are a n offlcg ram and a *readr reckciief table * 
with which one can find the integral fringe orders ii sirig the 
fractional fringe orders for two wavelengths. This is essm- 
tially an extension of Pant’s work. Where as with Pant’s 



methcxi or the use of Babinet-Soleil compensator one can find 
the integral fringe order upto few fringes (about 6), with 
the nanogram or the ready -reckoner table one can determine 
the integral fringe order virtually upto any number. 



CHAPTER 6 


REsunrs 


the previ'^us chapters, different experimental 
technlcf.as for the deterainati® of varloue oharacterlstlo 
parameters were dlsouseed. Those were essentially built 
Oh the concept of light ellipse. These experimental tech- 


niques are very general in nature, m the present chapter 
we shall present some experimental results in support of 
these techniques for completeness. In addition to these, 
we have already described some experimental results in tL 
previous chapter, connected iJlth the determlnatlcn of inte- 
gral fringe orders in photoelasticity. 


Test Model : 

Tests were conducted on a rectangular bar under 
combined torsion and tension. Fig. 67 shows the model and 
the light path under oonsideratlon Fig. 68. is the sketch 
of the loading device. Figs. 69. and 70. show the experi- 
m.,ntal set up and the loading arrangements, respectively. 
The model material was paraplex. The experiments were con- 
ducted with a 8" - lens type transmission polarlscope shown 

in Fig. 69. Two models with the following dlmenslais were 

used; '..i-''. 
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a = 1"; b - o.528»; c = 6i" 

( 11 ) a = l"j b = 1 / 2 ". <, ^ 

The light sources used ^.ere nercui-y and sodium vapour 


lai’ips* 

5oo4-Oa 


Three filters of \mvelengths 5520°a and 

were also need with white light. The model is checked 


for residual stresses with a laser (Helium-Neon; 63280 a, 5mw). 


The scattered light fringe pattern is observed ^en (i) the 
laser beam is along the optical path under consideration and 


(ii) the laser beam is perpendicular to this optical path. 

No fringes were observed in case Ci) . The fringe pattern 
observed in ca*se (ii) is shown in Fig. 71 . j^is shows that 
for the optical path under consideration no residual stresses 
exist. Also the residual stresses are normal to the surface 
of the model (i.e. parallel)to the optical path. 

The objects of the various experimental tests are as follows: 
(1) Experimental determination of the characteristic 


parameters for different loadings and checking tl§ 
Talidity of the governing optical equationsj 
Experimental determination of the characteristic 
parameters for different wavelengths of light and 
checking the validity of the governing optical 
equations; 

Tests in support of the existence of an optically 
equivalent system; 

Checking the validity of the assumptions involved 
in the previous scattered light methods^* ^ *^^or 
different loadings. 
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-he details of these tests ^vill be now presented 

(1 ) Experimental Determination of Characteristic Parameters 
Eor various Loadings (Refer Chapter 3) 

ExpBnLuient 3.1 Ppo cg^dnrp* 

Method 1 : Load the modal and keep the first quartar-waTo 
Plate axes at 45° to the polarizer axis so that we have a 
ciroalarly poaarlzed light at entrance to the model. With 
the second quarter-wave plate out of the field, rotate the 
analyser for minimum Intensity and note this direction . 
This direction corresponds to the minor axis of the light 
ellipse. Along the directions D, ±45°, the amplitudes of 
the components of the light vector are equal and these dire- 
ctions correspond to the secondary characteristic directions. 

Now, Introduce the second quarter-wave plate so 
that its axes are ParaUel to n, and + 90 ° directions. 
Rotate the analyser for minimum Intensity which will be 
absolute minimum. Let the minimum Intensity position be Dg. 
Then 2 (Dg - D^) = (90 + S) where S Is the eharaoteristlc 
retardation. Now remove both quarter-wave plates and keep 
analyser axes along one of the secondary characteristic dl- 

^©ctioris* Rots."!}© ib© polupi^pT* fnT* ^ • • 

u le jr-uxurizer lor (absolute) minimum inten- 
sity. The final orientation of the polarizer axis corres- 
ponds to one of the primary characteristic directions. 

Method 2: Let the optical set up consist of the polarizer 
(at arbitrary orientation), the stressed model, the quart er- 
-wave plate (at arbitrary orientation) and the analyser. 
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Rotate the analyser for minimiim intensity thus finding the 
orientation of the major and minor axes of the light ellipse* 
Rotate the analyser through 4-5o in either direction and note 
the intensity with a photometer. All ellipses (of 

different orientations and ellipticities ) will have an inte- 
nsity equal to along directions which make an angle of 
+ 45 ° with their major axis. Now keep the angle between 
the .axes of the second quarter-wave plate and analyser 
equal to 4-5® and rotate the polarizer until the intensity 
from the analyser is equal to the previously found value of 
The final orientation of the axis of the polarizer 
corresponds to one of the primary characteristic directions. 
Now remove the second quarter-wave plate and rotate the 
analyser for absoliite minimm intensity. The final orienta- 
tions of the analyser axes correspond to the secondary cha- 
racteristic directions. The characteristic retardation can 
be found as in the first method or by any other method. 

The characteristic parameters for various loadings 
are presented in Table 6.1 . Also presented in this table 
are the theoretical values. All the theoretical values 
presented in this chapter are found by integrating numeri- 
cally the governing optical equations in terms of Stokes 
parameters. The numerical integration is done using the 

42 

Runge-Kutte-Gill method. The series solution, (Timoshenko ), 
is used to find the values of the shear stress. The compu- 
tatlGns were done using IBM 7C44 computer.* Also, the 
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fractional part of the eharact eristic retardation of the 
Phase difference at the exit along a particular set of axes 
are found directly by experiment, m most of the cases, 
the integral fringe order is one. This can be seen by coun- 
ting the fringe number at the middle point of the model frcm 
the edges where the stress is known. (Along the edges, the 
secondary shear stress is zero). By numerical integration 
of the optical equations also, the same value for the Inte- 
gral fringe order is obtained. 


We shall use the following abbreviations in this 

chapter. 

Th; Theoretical 

Experimental 

Towards right from the vertical (direction of 
tensile loading) when looking towards the source 

tile vertical when looking towards 

the source 
^et; Retardation 


Pri.ch.dir; Primary characteristic direction 
Sec. ch. dir: Secondary characteristic direction 
Tensile load 
'ii* Torque 

Pri.str.dir; Principal-stress direction 
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TABIxE 6,1 

Characteristic Parameters For Various Loadings 
Model dimensions; a = 1"; h = 1/2*’; c = 6-|-” 
Wavelength of| light used: 5893°A 
Material fringe constant: 139 Psi. in/fringe 



Load 

Pri.ch.dir . 

Sec .ch.dir . 

Ch 

. ret. 

Sec . pr , 

T en Sion Torsion Th. 

P T 

lb . lb . in . 

Ex. 

Th. 

Ex. 

Th. 

Ex. 

St r. dir 
at en- 
trance 
Th. 

96 

8 

37°L‘' 

31 °L 

37 OR 

35^- 

-3380 

332 ° 

27°L 

96 

12 

,?45ol 

41 OL 

450 R 

41 OR 

4010 

394 ° 

32 OL 

96 

16 

58°L 

■ 52 °!- 

58'%^V2°R 

1+44° 

445°’ 

T35°i. 

96 

24- 

120 L 8 . 50 L 

120 R 

80R , 

: 470 ° 

45761 

380 L 

90 

15 


50^1 

55°R 

53°R 

436° 

433 ° 

35°L 

100 

18 

71 °L 

63 °L 

71 °R 

63 °R 

00 

0 

472 ° 

36®L 


In Table 6.1, the agreement between theoretical and experi- 
mental values for the characteristic retardation is excelloat, 
the maximnin discrepancy being Jfo . In Chapter 3, we noted 
that in anti-syiiinetric case of stress distribution (like 
the present case), the primary and secondary characteristic 
directions will be inclined equally, but in opposite direc- 
tions to the secondary principal -stress directions at the 
entrance and exit (of the model) respectively. The values 
of the characteristic directions in the table, support 
this statement. The theoretical values for these directions 
computed numerically are synmetrlGral. about the direction of 
tensile loading with 1 to 3 degrees differehce. This may be 
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due to the effect of round off errors involved in the numer- 
ical integration. However in the tables, the average values 
of these are presented. The values of secondary principal- 
-stresses presented are calculated using the value of shear 
stress computed from the series solution. However, for some 
cases, these values are cheered using strain gages cn the 
surface of the model. The theoretically computed values for 
the sec mdaiy principal-stress directions exactly coincided 
with the experimental values. The maximum difference bet- 
ween the theoretical and experimental values of the charac- 
teristic directions is 8®. This can be attributed partly 
due to the errors involved in the numeid.cal integration 
and partly to the possible errors in. the loading alignment. 
However maximum care was taken to minimise the Matter. 

The governing optical equations are thus shown to 
predict the optical phenomena in the presence of the rota- 
tion of the secondary principal-stresses fairly accurately. 
The experimental techniques for the determination of charac- 
teristic iHrameters are successfully employed to determine 
these parameters. From the values obtained, we note that 
difference in orientation of the characteristic and secondary 
principal-stress directions increases as the torsional load 
is inereasedi 

(2) Experimental Determination of Characteristic 
Parameters For Various Wavelengths of Light 
Experimental Procedure; 

The quarter -wave plates of the polarise ope are for 
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the green filter (5^61 °A) attached to the camera of the pola- 
riscope. For othei wavelengths of light these wave plates 
cease to be quarter -wave plates. For green light (5^61 °A) 
the characteristic parameters are found by the methods des- 
cribed earlier. For other wavelengths of light, the two 
quarter -wave plates are removed and the compensator is fixed 
oust after the model. The model is taken out of the field, 
the polarizer and analyser are crossed and the compensator 
axes are rotated for minimum intensity (with some arbitrary, 
but nonzero retardation in the compensator). The axes of 
the compensator are then parallel to the axes of the polari- 
zer, After locating the axes of the compensator thus, its 
axes are kept at a certain angle to the axes of the pola- 
rizer and the retardation in it is varied until the intensity 
is minimum. The corresponding reading of the compensator 
corresponds to a zero for that particular wavelength of light. 
Few other successive readings are found in the same way which 
correspond to zero retardation or a multiple of Ptt retarda- 
tion. Thus the compensator readings are calibrated for 
that particular wavelength of light. Mow the retardation 
in the compensator is made 90 ° • The primary and secondary 
characteristic directions are found using the method 2. The 
polarizer axes are kept at to the primary characteristic 
directions. The ellipticity tan » of the exit light ellipse 
is found by Tardy’s method. The characteristic retardaticn 
8 = 2i» . The values of characteristic parameters for diff- 

erent wavelengths of light in Table 6.2 are found experimen- 
tally in this manner. 
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TABLE 6.2 

Characteristic ih.rameters For Different 
Wavelengths of Light. 

Load; P = I 09 Ib. T = I 9 Ib. in. 

Dimensions; a = 1 ", b = 0.^28”, c = 6-|-'' 


Wavelenght 
of light 
used 

Material Fringe 
Constant F 
?si. in/fringe 

Pri . 

Th. 

ch. dir. Sec. 

Ex. Th. 

ch.dlr 

Ex. 

, Ch. ret 
(fractio- 
nal value ) 
Th. Ex. 

5950 

158 

34r 

36R 


38L 

890 7^0 

5520 

147 

28R 

29 R 

28L 

26 L 

67° 620 

5oo4 

128 

1 L 

0 

1R 

2R 

50 . 5 ° 410 


In Table 6.2 there is excellent agreement between 
the theoretical and experimental values of the characteris- 
tic directions. The discrei^ncy in the theoretical and ex- 
perimental values of the characteristic retardation is 4-^. 

(if the integral fringe order of 1 is added to the fractional 
fringe orders presented in the table. The integral fringe 
order is found experimentally and by numerical integrsjjtion 
to be 1 . When the wavelength of light is changed from 
5004°A to 59^0° Af the orientaticmi of the characteristic direc- 
tions changed by 34°, which is quite a marked variation. 

Thus the data (characteristic parameters) for different 
wavelengths of light can be used to obtain more information 
about the stress distribution along the optical path. The 
values in Table 6.2 aSso shew that the governing optical 
equations dan predict the optical ihenomena sufficiently 
accurately. It may be noted that the wavelength and the 
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material fringe constant donot follow any linear relation- 
ship. The material fringe constant is found experimentally 
within an experimental error + 2 (i.e. about + 1,5^). 

(3) Tests in Support of the Existence of an Optically 
Equivalent System 

It has been proved theoretically, that for a parti- 
cular wavelength of light, any optical path (or a jart of 
it) in a photoelastic model is equivalent to a system with 
a single retarder and a pure rotator. The characteristic 
parameters have also been defined for this optically equi- 
valent system i 

So, the following statements hold true: 

(i) When a certain amount of retardation is added along 
the primary characteristic directions at the entrance 
to the model, it appears as an additional retardation 
along the secondary characteristic directions at the 
exit of the model. 

(ii) Euring this process of addition of retardation, the 
amplitudes of light vector components along the sec <m- 
dary characteristic direeticns will not vary. 

(iii) Using the experimentally determined characteristic 
parameters, one can predict the exit light ellipse 
parameters for any tjrpe of incident light* 

(iv) If the chafact eristic retafdatioh of the optical path 
under con siderati on is nullified with a compensatcsr, 
the optical path behaves as a pure rotator* If the 
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i , 

rotatory power is millified then the optical path 
behaves as a single retarder. 

Conversely, experimental results which confirm these 
statements form supporting evidence for the existence of an 
optically eqiaivalent system. 

For a number of cases (of various loadings) the 
characteristic retardation is nullified with a compensator 
and then the optical path is then found to be acting as a 
Pure rotator. This confirms that the optical path is equi- 
valent to a single retarder and a pure rotator. 

Tables 6,3 and 6.^ show the variation in the inten- 
sity and retardation along the secondary characteristic dire- 
ctions when the retardation along the primary characteristic 
directions is varied. For this the characteristic parameters 
are found for the particular loading. The axes of the first 
quarter -wave plate are kept at to the primary characteris- 
tic directions. If the angle between the axes of the polari- 
zer and first quart er-waTO: plate is e , the retardaticxi along 
the primary chamct eristic directions at the entrance to the 
model is 2 e. The analyser is kept with its axes along the 
secondary characteristic directions. The polarizer is rota- 
ted so that the retardation along the primary characteristic 
directions is varied and the variation of the photometer 
reading is noted. To find the retardation along the secon- 
dary characteristic directions for a particular valpe of 2 ©,; 
the axes of the second quarter-wave plate are kept at 4-5 
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to the secondary characteristic directions so that we obtain 
a line ellipse before the analyser. The analyser is then 
rotated for minimum intensity. If the angle between the 
axes of the analyser and second quarter-wave plate is u , 
the retardation along the secondary characteristic direc- 
tions is 2(j) . Thus the 03 values are f oxond experimentally 
for various values of e . 

In Table 6.5 the parameters of light ellipse at 
the exit for various types of incident light are presented. 
Two sets of values are presented in this table. One set 
corresponds to the values found directly by experiment. The 
second set of values is computed from the characteristic 
parameters which are found experimentally. The light ellipse 
parameters that are tabulated are the phase difference and 
amplitude ratio along , the secondary principal-stress direc- 
tions at the exit of the model. These directicais are com- 
puted using the series solution for the shear stress; To 
find these parameters experimentally, the following procedure 
is adopted. ■ • 

Remove the second quarter -wave plate and keep the 
analyser axes at some inclination to the secondary principal- 
-stress axes at the exit. Introduce a compensator between 
the model and the analyser with the axes of the compensator 
along the secondary principal-stress axes. Vary the retar- 
dation % of the compensator until the intensity of light 
from the analyser becomes a minimum. The phase difference 
along the secondary principdl-st res s axes at the exit is 



264 


then - 5^, Now, we have plane polarized light after the 

compensator, the direction of which can be found by rotating 
the analyser for minimum (zero) intensity. If the angle 
between the axes of the secondary principal-stresses at the 
exit and the analyser is to , the amplitude ratio along the 
secondarjr principal-stress axes is tan oj . The proof for this 
is as follows: 


Lot the Stokes vector before and after the compensa- 
tor with respect to its axes be given by V-j and V 2 which 


are: 











S 2 

and = 

S„ Cos 5 - S_ Sin s 

^ c ^ 




^3 c ^2 5 c 


when V 2 is referred to the axes of the analyser it becomes 


S-j Cos 2 © + Sin 20 (S 2 Cos S S^ Sin 5 ^) 


“nr 

Si 

(S 2 Cos 3^ - S 3 Sin 3q) Cos 20 - S-^ Sin 2e 

E 

S2 

S 3 Cos 3^ + S 2 Sin 

: 

h_ 


where 0 is the angle between the axes of the compen- 
sator and the analyser. 


93, 


= Sin 20 (- S„ Sin 6 - S- Cos 8^ ) = 0 

a C 3 C 

- S„ 


if Sin 20 = 0 or tan <5_ = — ^ = tanC i.e. 

S2 


_ e 


where € is the phase difference along the secondary 
principal-stresses of the -light ellipse from the model. 

Also Sin ^ f O. 
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TABIE 6.^ 

Variation of Intensity and Retardation Along the 
Secondary Characteristic Directions When the R etar - 
datl«on Along the Prima ry Char acter i stic Directions 
is Varied 

Dimensicns; a = 1 ", b = 0 . 528 ", c = 6i" 

Load: P = IO9 lb. T = I9 lb. in. 

Wavelength of light used: 5^61 

Material fringe constant; 1^7 Psi. in. /fringe 

Characteristic Parameters: 



Ex. 

Th. 

Pri, Ch. Dir. 

630L 

6 2^1 

Sec. ch. dir. 

60 ®R 

62 °E 

ch. ret. (fractional value 56 '^ 

67° 


The intensity variation along one of the secondary charac- 
teristic directions when the retardation along the primary 
characteristic directions is varied is 2 divisions on 100 


scale of photometer. 


Ret .along 

Ret . observed difference 

Ret, along 

Ret.obser- 

diff- 

pri.ch.dir . 

along sec. 


pri.ch.dir 

. ved sec. 

ere-- 

(deg.) 

ch.dir . 

(deg. ) 

( deg . ) 

ch.dir. 

nee 


(deg.) 

(deg.) 

( deg . ) 






2 h 

0 

56 

1 8 

too 

162 


20 

7^ 

20 

tao 

18^ 

18 

ifO 


21 

IlfO 

202 

18 

60 

115 

23 

160 

220 

20 

Bo 

138 

tao 

^0 
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TABLE 6M- 

Variation of Intensity and Retardation Along 
the Seccndarv Characteristic Directions When 
the Retardation Along the Primary Characteristic 


Directions is Varied 


Dimensions: a = 1 , b = 0.^28 , c = 6-1- 

Material fringe exist ant; 158 Psi^in. /fringe 
Wavelength of light used; 5893°A 

Load; P = 96 lb. T = 16 lb, in. 
Characteristic Parameters: 

Ex. Th. 

Pri.ch.dir. 38.5°R 38.5^ 

Sec. ch. dir. 36OL 38.5°R 

Ch, ret .(fractional value) -44-° -440 


Ret. along Ret. along diffe- Ret. added Ret. observed diffe- 
pri.ch. sec. ch. dir. rence along pri. along sec. rence 



With analyser axis along one of the secondary prineipal- 
-stress directions at the exit, the intensity reading cn 
100 scale of photometer varied from 21 to 23, when retarda- 
tion along the primary characteristic directions varied* 
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table 6 .? 

faraaeters.of Light mnp o, ^xit. 
t he Optically Ti g ulvalfint Modal 
Load: P = 100 li. J ^ 18 lb. m. 
Wavelength of light used: 5 if 61 °A 
Dimensions: a = l" , b = 0.528" , c = 6i" 

Material fringe constant: Ilf? Psi. in. /fringe 
Character istic 


Pri, ch. dir. 

Sec. ch. dir, 

Ch. ret. 

Sec. pr.str.dir. ) 
at entrance ) 

Sec. pr.str.dir. ) 
at exit . ) 


- 116 ° 

65 °L 

65 °R 

35 °L 

35 °R 


- 100 ° 
61 .50L 

627R 


Ret .at Ret .alcng seFI AmT 

entrance Pr.str.dfr.at 

along sec. exit ^ * 

fdel^T^^* f?' A calculated Ex. 

(.deg.; (deg. ) from ch. par 

(deg.) 


AmpHtude ratio along sec. 
pr.str. dir. at exit 


calculated ftom 
ch.par. 


1 .54 
2.05 
2.41 

3.08 

4.01 

4.01 

2-75 

2.14 

'1 .73 

1 .28 
1.00 
0.78 
0.^ 


1.50 

1-93 

2.^ 

3.34 

4.00 

3.78 

2.97 

2.25 

1.73 

1.35 

1 .07 

0.85 

^.67 
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Prom Tables 6.3 and 6.^, we note that the intensity 
variation along the secondary characteristic directions when 
the retardation along the primary characteristic directions 
is varied is 2 divisions on the 100 scale of photometer. 
Since the ihotometer pointer by itself shows a fluctuation 
of 2 divisions on 100 scale, this can be taken to mean that 
the amplitudes along the secondary characteristic directions 
remain constant. When the retardation along the primary 
characteristic directions is varied in steps of 20^^, the 
retardation along the secondary characteristic directions 
varied in steps of 18° to The error involved in the 

measurement of phase difference is 5° in the present experi- 
mentation. i.e. the experimental error involved can be as 
high as 5° with the present experimental set up. So, what- 
ever retardation is added along the primary characteristic 
direction the same appears as an additional retardation 
along the secondary characteristic directions. In Table 6.5 
the two sets of values of the phase difference C an<i 
amplitude ratio 0 „ are compared* One set of values is 
obtained direcbiy by exfferiment aiftdl the other set of values 
is calculated from the (experimetlt ally determined) charac- 
teristic parameters. The agreement between the two sets 
of values is good. Thus all these experimental resdltai 
confirm the statements made earlier regarding the optically 
equivalent system. These form a supporting evidence for 
the existence of an optically equivalent system f or an^^^^^ 



optical path in a three-dimensional photoelastic system for 
a particular wavelength of light. 

(4) Experimental Check for the Validity of Assumptions 
Involved in the Previous Scattered Light Methods: 

We note that the assumptions involved in the pre- 

1011 12 

vious scattered light methods ’ ’ are satisfied when the 

secondary principal-stress -axes at the entrance and exit 
coincide with the primary and secondary characteristic 
directions . 

Tables 6. 6 - 6.10 show the amplitude ratio e and 
the phase difference € along the secondary principal-stress 
axes at the exit v;hen the retardation is varied alcng the 
secondary principal-stress axes at the entrance. Also shewn 
are the values of e„ and £ which are computed from the 

iX. 

(experimentally determined) characteristic parameters. These 
sets of values are presented for ^ loadings. 

To vary the phase difference along the secondary 
principal-stress axes at the entrance the axes of the first 
quarter -wave plate are kept at 4-5° to these secondary 
principal-stress axes. The polarizer is then rotated to 
vary the retardation along these axes. If the angle between 
the axes of the polarizer and the first quarter-wave plate 
is , the retardation introduced along the secondary prin- 
cipal-stresses at the entrance is 2 C><. The azimuth of the 
light ellipse at the exit is found with the analyser. The 
axes of the second quarter-wave plate are kept parallel to 
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the axes of this light ellipse so that it reduces to a line 
ellipse. The orientation of this line ellipse is foimd with 
the analyser. From this orientation the ellipticity and the 
azimuth of the light ellipse from the model are found, The 
values of e^^. and € are calculated from the larameters of 
this light ellipse. The characteristic parameters are found 
experimentally by the methods described earlier. The values 
of G and e^ for various types of incident light are also 
computed using these parameters. Sodim vapour lamp (589 3°A) 
was used as the source of light. Only the quarter -wave 
plates of the polariscope were used for the measurement of 
ellipticity etc. These introduce only 85° of retardation 
for wavelength 5893°A and not 90°. 

Wavelength of light used: 5893°A 

Material fringe c^onstant; 139 Psi. in. /fringe 



In the tables to(foolow^ the following abbreviations are 

: Retardation added along secondary 
principal-stress axes at entrance* 

: Intensity of light along on§ of the 
secondary principal -stress axes at exit 
as given by the reading on 1000 scale 
of photometer. 

: Azimuth of the out coming light ellipge^^^^^^^^^^^ 
with respect to the vertical (life© of 
tensile loading) 

Q.w.p. : Second quarter-wave plate position 


used: 

Ret . ent . 

Intensity 


Initial azimuth 
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Final azimuth: 


Ph.diff . 
(Ex.) 


Ri. diff. 
(Calc. ) 


Ellipticity 

(Ex) 


Ellipticity 
(Calc . ) 


Ampl. ratio 
(Ex,. ) 


Ampl. ratio 
( Calc . ) 


Pinal azimuth of the outcoming line ellipse 
after the introduction of the second quart er- 
-wave plate. 

Phase difference along the secondary princi- 
pal - stress axes at exit calculated from 
the exit light ellipse parameters which are 
found experimentally. 

Phase difference along the secondary princi- 
pal-stress axes at exit calculated from 
characteristic parameters which are found 
experiment ally. 

Ellipticity of the outcoming light ellipse 
found directly by experiment. 

Ellipticity of the outcoming light ellipse 
calculated from the characteristic parameters 
which are found experimentally. 

Amplitude ratio along the seccndaPy princi- 
pal-stress directions calculated from the 
exit light ellipse parameters which are 
found experimentally. 

Amplitude ratio along the secondsiry princi- 
pal-stress directicais calculated from the 
characteristic parameters which are found 
experimentally. 
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T.^LE 6.6 


P = 90 lb. T = 15 ' Ib^ in. 


Ret . ent . 

Initial azimuth 

1; •W'* p 

Final azimuth 

0 

0 

-d- 

500 L 

T^R 

1+50 

72 Or 

27OR 

46 °R 

60° 

91 °R 

460r 

72 ^^ 

900 

79OL 

569ei 

90Ol 

1 20° 

76 '^L 

59 ®R 

720l 

135 '=' 

71 '='L 

61 fOR 

58ol 

1500 

65°L 

700r 

440L 

1800 

55 '^L 

800r 

23% 


TABIxE 6.7 

P = 96 lb. T =12 

lb . in . 


Ret. ent. 

Initial azimuth. 

Q.w.p, 

Final azimuth 

0 

760 L 

59 ®R 

520L 

30 ° 

67 °L 

68^^ 

35 ^L 

1+50 

55 °L 

SoOr 

1704 

0 

0 


850R 


90P 

93 °R 

ifSOR 

71 

1200 

S'"!' 

480l 

11 OR 

t 35 '^ 

30R 

420 L 

IQOp 

150° 

20R 

lf 4 -OL 


ISO'^ 

50 j^ 

ifO°L 

SPl 



TABLE 6,8 



P = 96 lb. T = 

8 lb. in. 


Ret. ent. 

Initial azimuth 

Q.-w.p, 

Final azimuth 

0° 

35 OR 

4oo.i 

71 °R 

300 

85 OR 

40OR 

95°L 

450 

860 r 

41 Or 

87°L 

600 

83 

38OH 

81 OR 

90° 

8 >+or 

39^1 

630 L 

120° 

120 h 

33^L 

560r 

135° 

9°L 

540L 

^OR 

150^ 

80l 

53 % 

20 Or 

1800 

7°L 

52®L 

40 R 


T.ABLE 6«9 


P = 96 lb. T = 16 lb. in.. 


Ret.ent. Initial azimuth Q.w.p.. Final azimuth 


0 

380 l 

830 L 

IOOR 

300 

17°L 

620 L 


0 

0 

\D 

50L 

50OL 

fOR 

900 

40.R , 

4 iOl 

4OE 

120° 

I6OR 

29OL 

3OR 

135 ° 

I70R 

28OL 

0 

1500 

28OR 

180 L 


1800 

56 °R 

IIOR 

280R 



TABLE 6.1 Q 
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P = 96 lb. T = 24 lb. in. 


rfet • ent . 

Initial azimuth 

Q.-w.p, Pi 

nal azimuth 

0 

40R 


41ol 

120 L 

300 

220 R 


24ol 

50R 

4-^0 

280r 


1-701 

14or 

60O 

330R 


70L 

270R 

90 ° 

48or 


3 °R 

450R 

1200 

660h 


210 a 

730R 

1500 

72OR 


270R 

840r 

I800 

840 L 


51 OR 

650L 





-- _ 



TmiE 6.11 




P = 90 

lb, T 

= 15 lb. in. 


It 0 0 • 0 nt • 

Pn 

Ex. 

. dttff. ’ 
Calc • 

Ellipt icity 

0 

770 

840 

0.65 

• 

0.53 

450 

•660 

520 * 

0.49 

0.49 

60 ° 

399 

• 3 ° 

0.34 

0.38 

90° 

220 

210 

0,19 

0.15 

1200 

80 

140 

+0.07 

+0.07 

135° 

280 

380 

+ 0 .A 

+0.20 

1500 

460 

620 

+ 0»38 

40.31 

1800 

73 ° 

840 

-fO #63 
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TABLE 6.12 

P = 96 lb. T = 12 lb. in. 


Ret. ent.' 

Fh' 

Ife. 

.diff. 
cal* ... 

E 

Ex. 

llipticity 

Pa 1 

0 

48 ^ 

39 ° 

0.45 

OdXC « 

0,29 

300 

650 

710 

0.63 

0.59 

450 

3 oo 

85 ° 

0.78 

0.76 

600 

87° 

820 

0.81 

0 , 84 - 

90 ° 

340 

57 ° 

0.29 

0.53 

1200 

.320 

+290 

0.25 

0.24 

135° 

15° 

+140 • 

0,12 

0.11 

1500 

110 

3 ° 

0.09 

0*02 

180° 

270 

39 ° 

0.23 

0.29 


TM:LS 6.18 

P = 96 lb. T = 8 

lb. in. 


Ret. ent - 

Fh. I 

liiiv • 

liff. 

__ Calc . 

Eliipticity 

0 

37 ° 

30° 

0,25 

GLJUU « 

0.25 

30° 

20 

20 

0.00 

0.02 

450 

22° ' 

I90 

0,14 

0.15 

60*^ 

46 ° 

340 

0.30 

0.29 

900 

7 +° 

63 ° 

0.66 

0.60 

•1 20° 

880 • 

89° 

0.97 

0.93 

135° 

740 

75 ° 

0.^ 

0.73 

1500 

670 

610 

0.53 

0.55 

1800 

31° 

300 

0.19 

0.25 



TaB^6.16 
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Het . ent , 


0 

450 

60^ 

90° 

1200 

135° 

1500 

1800 


^ - 90 lb . 1 lb » in. 


Intensity 

Am pi. ■ 
Ex. 

ratio. 

Calc. 

23 

1 .44 

1 .87 

13 

1.76 

0.99 

11 

1 .28 

0.81 

6 

0.97 

0.53 

4 

0.87 

0.37 

3 

0,75 

0.34 

5 

0.68 

0.37 

8 

0.71 

0.54 

table 6.17 

96 lb. T = 

12 2b, in. 



Rsb * STit « 
0 

300 
^50 
600 
90° 

1200 
1350 
1500 
1800 


Intensity 

8 

10 
12 
1^ 

20 
23 

23 
25 
28 


Ampl. ratio. 

Calc . 


.0.98 

0.86 

0.84 

0.82 

;i'^.35 

0.60 


1 

1 .49 
1 .30 
"1.13 

0*64 


9.71 0.58 

0.67 0.54 

0.78 0*51 
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TABLE 6.18 



P = 96 lb. T 

= 8 Ib.in. 


Ret . ent . 

Int en sit y Ampl . 

Ex. 

fat io 

Calc. 

0 

11 

2.11 

1 .if 8 

300 

11 

2.53 

1 .51 


11 

2.29 

1 .4-5 

600 

13 

2.18 

1.36 

900 

17 

1 .36 

1 .14 

1200 

2h 

0.99 

0.93 

135'' 

27 

0.70 

O.Sf 

1500 

29 

0.62 

0.77 

I 800 

32 

0.4l 

0.67 


TABLE 6.19 



P = 96 lb. T 

= 16 lb. 

in • 

Ret .ent . 

Intensity 

impl. 
_ . ,_£Ss 

ratio 

Calc. 

0 

8 

0.55 

0.57 

300 

14 

0.51 

0,41 

60O 

20 

0.61 

0.39 

900 

25 

0.81 

_ o_. 53 \ 

I200 

28 

/„.-l ',20 

I f 

0.79 

135" 

27 

1,22 

■ 0.96 i"" 

1500 

27 

,1.47 

1 .17 K 

1800 

23 
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TABLE 6.20 




P = 96 lb . T = 24 

lb . in. 


Ret, ent. 

Intensity 

Ampl. ratio 
lx. Calc. 

0 

26 

0.92 

1.53 

300 

32 

1 .54 

2,76 

450 

34 

1 .97 

3.94 

60° 

35 

3.17 

5.23 

900 

33 

10.76 

3.40 

120® 

28 

3.65 

1 .80 

1500 

19 

2.38 

1 .08 

180® 

12 

1.03 

0.65 

In 

Tables 6.6 6.10, the 

experimental values of 

the azimuth 

of the light ellipse at 

the exit 

and its orien- 


tat ion after it is reduced to a line ellipse are presented* 
From these values, the elHpticity, the amplitude ratio and 
phase difference along the secondary principal-»stress direc- 
tions at the exit are calculated. In Tables 6.11 -6.15, 
these values of amplitude ratio and phase difference are 
compared v/ith the values calculated from the characteristic 
parameters (which are found experimentally). In Tables 6.16 
6.20, the intensity along the secondary principal-stress 
directions (as given by the photometer reading) and the amP-^ 
litude ratio along the secondary principal-stress directions 
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at the exit-(ij3§,-calculated from the light ellipse para- 
meters at the exit and (ii) calculated from characteristic 
paramet ers.,,,., are presented. Frcm these tables we note that 
both the sets of values of the ellipticities match very veil, 
whereas for phase difference and amplitude ratio, the match-*, 
ing is not good. The reasons for the discrepancy in the later 
case are! (i) the vave plates used introduce only 85° of re- 
tardation for Wavelength of light 5^893*^1 (ii) errors in 
measurement of the orientation of light ellipse severely 
afj.ect the values obtained for phase difference and amplitude, 
ratio since the constant ellipticity curves .are very steep 
for azimuth values near 0 and 90'^, Fig. 1+4- • (iii) errors 
in measurement of orientation of light ellipse do not have 
a significant effect on the values of ellipticity obtained* 
Also, we note from these tables that the assui&p— 
tictis of the previous scattered light methodl®»^”’»*^^ are not 
valid for the loadings ccosidered. VJhen the retardation 
along the sec. pr. str. directions at the entrance is varied, ♦ 
at . the exit, the intensity along the secondary prinGi|B.l- 
-stress directions does vary and the phase difference 
these directions does not increase by the same amotaat as thp. 
retardation added at the entrance, ; ^ 

In this chapter, we have presented experimental ■ 
evidence for a number of cases in support of the. following 
statements; 



The total light path in a three-dimensional photo- 
elastic model for a particular wavelength of light 
is equivalent to a system with a retarder and a pure 
rotator. The three parameters of this equivalent 
S 3 ’'stem called the characteristic parameters can be 
found by the techniques described earlier. 

The governing optical equations (derived in Chapter 2) 
describe, the optical phenomena in three-dimensional 
cases where rotation of the secondary principal-stress 
axes along the optical path is present. 

The variation of the characteristic parameters with 
wavelength of light e;an be quite marked in some cases. 
In the case considered, the characteristic directions 
varied by 3^^ when the wavelength of light is varied 
by 

The assimiptions involved in the previous scattered 
light methods are not valid when the rotation of 
the secondary axes is large. These assumptions 
lead to large errors in such cases. 



CHAPTER 7 


COHCLUS ION 

This thesis which consists of six chapters in 
addition to the present one deals ’■^ith the problems of , 
stress -induced optical activity and determination of the / 
optically equivalent model in three-dimensional photo- 
elasticity. Chapter 1 , the introduction, describes the 
general problem when the secondary principal-stress axes 
and their magnitude continuously vary along the light path 
in a stressed photoelastic model. It gives a brief review 
of the stress-optic laws in two - and three-dimensions and 
attempts to modify these in the presence of optical activity 
It also discusses briefly some of the scattered light tech- 
niques with specific relevance to the effects of optical 
activity. Finally, the objectives of the present investiga- 
tion are stated along the methods of approach. 

In Chapter 2, the concept of light ellipse is 
introduced. It is shown that studying the properties and 
transformations of light ellipse and then using this infor- 
mation in photoelastic analysis has a definite advantage ^ 
over the usual procedure of resolving the light vector into 
its components trigonometrically along various direct i<ma,s 



Of f ' “ Ihy.lcal asp^t 

tne problem involved readily. The ll^ht .11- 

in te«ns of stov-, 

sector and throughout thio ■ 

satlon the concept of u-t el, • an^esti- 

toen used. vector have 

». .-.e.“"ir::::.’'“” “‘~‘™ « 

tion Of secondary principal-stresi 

that ^->l. , . ^‘^^ess axes is presented to show 

hat the kinematically derivoa ^ ^ 

y at^nvecl equations nrodinf -f-io 

Phenomena sufficient ^ ^ ^ ^ optical 
iciently accurately. ThpoT.^- ^ 

a- then derived kinematically m teims of 

-- out to he a eet of three ^an:!: 

linear different t d co a.. order 

These ■ nations with variable coefficients 

Ihese are very convenient for analirsi. . 

IIsin» ev, uulysis and application. 

mg these equations the analv+inoT 

naj.3rticai solution fnv 

where ratio R of the rat, of r , 

arbltrarv f e- i ® ^ uny 

y unction) to the rate of retardation ?' • 
tant. The case d, ^ is a cons- 

h f. .. r “^ucter and Mlndlln where 

hoth. and r are constants is a particular case Of: this^ 

lence of a^lririp 'V'’' 

^ s3^stGm cc3ntaining a 

e arder and a pure rotator is presented. Use of Poincare' 

P ere representation of polarised foiias for this proof is 
also dencnstrated. The proof is then extended to any length 

path in a three,-dimensicnal^^I^ 

that is, any length of optioal path i„ ^ three -di.ensiona; 
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photoelastic model is again equivalent to a system contain- 
ing a single retarder and a pure rotator. A discussion of 
some special cases (case of a pure rotator, anti-s 3 rmriietric 
case, axi -symmetric case, reversal of optical path, reflec- 
tion polariscope) in which there exists a relation between 
the three characteristic parameters, is presented. 

Two simple and direct methods of determining experi 

mentally the characteristics parameters are presented. It 

is to be noted that there donot exist any such methods to 

find the characteristic parameters directly in a simple •'.■ra.y. 

However Kent and Lawson ^ suggested an ingeneous device to 

measure the ellipticity directy. This has been used by 

1^,15,16 

Robert and Guillemet to find the characteristic para- 

meters. But this method requires the rotation of the analy- 
ser at a particular speed and some additional instrumenta- 
tion. 

In Chapter a review of some of the previous 
scattered light techniques, the assumptions involved and the 
analysis, is made. A critical analysis of these methods is 
made making use of the analytical solution (obtained in 
Chapter 2) for the case where R = is a constant (0' 

can be any arbitrary function) and numerical solutions (of 
the governing optical equations) for some cases. The limi- 
tations of the previous methods are discussed. 

The scattered light techniques essentially involve 
the following steps; (1) measurement of scattered light 
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ntensity along different dlreotlcns taking into account, 
he variation of the lengths of light path Involved, (2) 
fading tte parameters of the light ellipse and the oharact- 
ameters at the point under consideration ( 3 ) eva- 
uating the functions ^ and ^ using these values and the 
g eming o;*ical equations. Some new scattered light tech- 
niques are presented which donot Involve any assumptions or 

he rotation of the analyser at a nar-h • t 

ar a Particular speed. 

Chapter 5’ deal-^ ir-it-v, o 

' 3 mam problems, (i) graphical 

methods in Phot oelasticitv 

, , S.icity, (11) theory of compensation and 

(ill) determination of Integral fringe order In photoelas- 

ticity. Graphical methods are Very convenient to find the 
integrated optical effect in the presence of the rotation 
Of the secondary principal-stress axes, a review of the 
earlier graphical methods Is made. Of these the stereo- 
graphic and paraUel projections are derived from Poincare' 
sphere representation of polarised forms. The methods using 
quaternions needs geometrical construction for each problem. 

In our discussion a new approach to the graphical methods 
is described which doesnot Involve the use of Poincare' 

Ihe light ellipse has essentially four parameters 
the Phase difference, eUthtlclty, amplitude ratio and 
azimuth. These are connected by specific relations which 
can be used to plot two of these fo.ar parameters keeping a 
third parameter constant. Based on these, two graSUoal 
methods are presented. One of these Is veiy useful In find- 
ing the effect of errors In various polariscopes. 
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The second method is ^ery simple to construct since It 
inyolyes only straight lines. Etocedures are described to 
determine the Integrated optical effect and the characte- 
ristic parameters using these graphical methods. 

The second problem deals .rtth the theory of com- 
pensation. nse Of (a) a system of four retarders each hay- 
ing the same retardation ulth a polarizer or an analyser 
and (b) a system of eight retarders each having the same 
retardation, as a large field compensator is described, 
ccnoept of an optically equivalent system is made use 
In deriving this. The angles between the axes of diff- 
erent retarders in either system (a) or in (b) are to be 

adjusted such that the rotatory power of the equivalent 

system for this system of -k ^ 

oytooun oi retarders becomes zero. The 

princip.Ie iryolvad is esseritTai itr an q-w-i- • 

essentially an extension of the work 

uy Pancharatnam ^8 Tuzi^9 , 


The third problem deals with the determination of 
integral .fringe orders in two-dimensional photoelastic 
t^odels under load. The drawbacks of the present methods 
are discussed. Construction of a nomogram and a ready 
reckoner table with which one can determine the integral 
fringe oraers virtually upto any number. This procedure 
doesnot have the limitations of the previous methods. This 
is essentially an extension of Pant 's 80 work in which the 
integral fringe orders are determined from the fractional 
fringe orders for two wavelengths of light when the two 



288 


total fringe orders donot differ by more than unity. 

Experimental results are presented in Chapter 6 
as supporting evidence to the theory developed in the pre- 
vious chapters. The experiments were conducted with a rec- 
tangular model under combined tensicn and torsion in a 
transmission polariscope. These form the supporting evi- 
dence for the following statements; 

(1 ) The governing optical equations in terms of Stokes 

parameters derived in Chapter 2 predict the optical 
phenomena in the presence of rotation of secondary 
principal-stress axes (along the optical path) 
sufficiently accurately. 

(2) The techniques described in Chapter 3 can be 
successfully used to determine the characteristic 
parameters. 

(3) The change in the characteristic directions is 

% 

found to be 3^° when the wavelength is changed 
from ^OOM-^A to 5950 °A (the corresponding change 
in the material fringe constant is from 128 -158 
Psi. in. /fringe). This is quite a significant 
variation and in cases such as this, the data for 
different wavelengths of light can be used to 
obtain more informatico about the state of stress 
along the optical path. 
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Any part of optical path In a three-dimensional 
Photoelastic model is equivalent to a system 
oontamang a retarder and a pure rotator. 

The assumptions on which the existing scattered 

light techniques are based, are not valid In a 
general case. 




APPEIJDIX A 


(i) Determination of the Characteristic Parameters 
When the Direction of Observation is Fixed 

From Eqs.- (^.4-2) 

Cos 2il» + (S^ Cos p^- Sinp^ ) Sin 2^1^ 

If = Cos 26, $2 = Sin 20 CosC , = Sin 26 SinC 

:r Cos 26 Cos 2il» + Sin 2i}» (Sin 2e CosG Cos p-^ 

- Sin 2e Sine Sin P>j ) 

d „ Sin 24» (Sin 2 q Cos Sine + Sin 2e Cose Sin ) 

d C 

The above expression is equal to ?:ero when 
tan€ - - t onp,^ 
or C = ” P '1 

This moans that if the retardation C is varied 
along the primary characteristic direction, the intensity 
along any arbitrary direction at the exit of the optically 
equivalent system becomes extremum when C = - p^ where p-j 
is tha oharaot eristic retardatl® of the eaulvalent optical 
system- Let the angle between the axes of the polarizer - n 
the guartor-wavo plate be « - Mote the intensity I^ . How 
rotate the polarizer through 90° and note the intensity Ig* 
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Change =<, the angle between the axes of the polarizer' and 
quart er -wax’-e plate (keeping, the quart er-\-/ave plate fixed) 
and Tepeat this procedure until In this position, 

the major axis of thelight ellipse at the point under consi- 
deration makes an angle of ^5° ^Arith the direction of observa- 
tion. (i) One possibility is that this light is circularly 
polarized. If it is so, the amplitudes along the secondary 
ch.3-ract eristic directions are equal and the amplitudes along 
the primary characteristic directions are equal. But the 
light ellipse after the quarter -'.rave plate will have its 
major axis along the quarter-wave plate axes so that along 
directions making with these axes, the amplitudes will 
be equal. That is, the quarter-wave plate axes make ^5° mth 
the primary characteristic directions. The phase angle along 
these primary characteristic directions after the quarter- 
-wave plate, is 2=( . So, we have 2c< + p-j = tt/p. All 
these statements will be true only if the light at the point 
under consideration is circularly polarized, (ii) The second 
possibility is ’that the direction of observation is along 
';ne of the secondary characteristic directions, (iii) Neither 
of the cases (i) and (ii) is true. 

To check these, rotate the polarizer and note the 
intensity I. (i) If I doe snct change when the polarizer is 
rot.atod this falls into the second case so that the direc- 
tion ' 3 f observation corresponds to one of the secondary cha-: 
ract eristic directions. The axes of the quarter -wave plate 
are at ^5° to the primary characteristic directions. , 
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The characteristic retardation cannot be found unless the 
direction of observation is changed, (ii) If I changes when 
the polarizer is rotated, find that positi^on of the polari- 
zer for which the condition is satisfied. In this 

position let the angle between the axes of the polarizer 
and the quart er -vrave plate be . Rotate the polarizer for 
minimujii intensity and lot this angle be o( 2 * When the angle 
between the ,axes of the polarizer and quarter-wave plate is 
for obtaining I^ = I^ and if the light is circularly po- 
larized at the point under consideration = V2« 

VJhen the polarizer is r.tated for minimum intensity 2<2 + 

= 0 so that -^ 2 '^ = fj;* Change the orientation of the 
quart or-wato plate until this c-ondition is satisfied. When 
this is satisfied, the axes , of the quarter-wave plate coincide 
with the primary characteristic directions and = tt/r - 2 
= - 2 P^ 2 * To) find the rotat'ory power, nullify the characteris- 
tic retardation with a ctjrapensator and rotate the polarizer 
f':»r minimum intensity. The angle between the axes of the 
polarizer and the direction of observation will give the ro- 
tat'jry power. 
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and = - S2P’Cos2i3 - S-j p’ Sin20 (2.Mf ) 


where 0 is the orientation of the secondary 
principal-stress directions with respect to the fixed axes. 
From Eqs. (2«'+2) and (2.^3) 

S^' = 0 

when = 0 or p' = 0 or Sin20 = 0(i.e. 0=0 or TT/Z) 

Also 


= Cos2oj C os2 ^ 


Similarly, 


when 


s' = 0 

So = 0 orp’ = 0 or Cos20 = 0 (iiOi 0 ='^A 3^A - 


and 

Sg = Coa2w Sin2£>C 

S becomes when o(is replaced by ( '^/+) 

1 

(when w remains constant) 

. SL becomes S ' when \>C is replaced by ( 

• # T 

(when w remains constant) 

For all points where the intensity is extremum along 
the Z-axis (i.e. along the optical path) for a particular 
direction of observation 


dS^ 


= 0 


* I 

If at a point on the optical path in the model, 

oithur P*= 0 or p 0 or both P’ and S3 are zero. ^ To eli- , 

minate the case where S 3 = 0, the incident Ixght is to be 
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modified. If =82 = 0 for any type of incident liglit we 
can conclude that at the point under consideration p’ = 0. 

At this point the directions of the secondary principal-stress 
axes can not he found since no retardation is being added at 
this point. If 8 ^’ 0, but ^ 0, the direction of obser- 

vation is along one of the secondary principal-stress axes. 
(Because 8 ^ ^ 0, P’ ^ 0 otherwise 8 ^' = = 0. But 0 = 0). 

This suggests a method of finding the directions of the secon- 
dary principal-stresses directly from the scattered light 
fringe pattern. 

The scattered light fringe pattern shows the in- 
tensities of scattered light along different directicos 

13 

at different points on the optical lath, 8 winscn et al . 

At each point find the directions of observation ^2* °^3 
etc. along which 8 -j =0. Find the value of along direc- 
tions *^ 2 , “^3 inclined at an angle of 45^ to the 

directions at the corresponding points on the 

optical path. Let the value of’ S^' be nonzero corresponding, 
to the direction say. Then the direction alcng 

one of the secondary principal -stress directions. A laser 
beam of a very small diameter will be an ideal source for 
obtaining the required scattered light fringe pattern. The 
offocts 3f ■varying lengths of (soatt.-rod or transmitted) 
light P'lth can be eliminated in the following way. For this 
take two scattered light fringe BiotograPhs one with right 
and the other with left eircularly polarized light at the 
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entrance. At any point and along a particular direction, 
the intensities in the two patterns are and then the 
required intensity I (in which the effect of varying light 
paths is eliiainatod) in the first pattern is 


I1+I2 


Because , 


I. = K 
^ X 


lo. = K 


1^+1^= K 
2 

= — = 

X K B, +I 2 


After knowing the value of 0 (one of the secondary 
principal -stress orientation) at each point, the value of p’ 
can be calculated by using the relation 

p- 

Sin 20 


The value of S^ at any point can be found in the following 
way. 

At a point on the optical i^th and along a particular set of 
reference axes X and I, 

= a2 - a2 

„ B - A 



where I. and l ar(= -f-vi 

1 J -2 are the intensities of 

when the incident light • ^ ^ scattered light 

J-rgnr is right and inf-f- 

respectively. circularly polarized 

Tind the value of d o '* a 

X direction. This wm • ^ '* 5 ° to the 

« ana s th " ^-Blreotlon. 

^elation 3 -‘■ound from the 

S3 = (1 - - s|)i 
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